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I. 


INTRODUCTION 


A.  BACKGROUND 

As  the  use  of  space  and  satellites  continues  to  grow  in  the  21®^  Century, 
technology  will  yield  smaller  devices  with  reduced  power  requirements.  These 
advancements  are  countered  by  our  growing  demand  for  faster  and  higher-bandwidth 
communications  and  will  predictably  lead  to  larger  satellites.  As  satellites  have  become 
larger  and  more  complex,  historically  the  launch  booster  capability  has  not  been 
appreciably  increased.  This  trend  continues  to  motivate  engineers  to  design  satellites  that 
are  lighter,  more  densely  configured,  and  capable  of  deploying  large  antennas  and 
reflectors  when  on  orbit. 

Another  result  of  our  continued  introduction  of  space  vehicles  and  the  degradation 
of  existing  satellites  is  the  drastic  increase  in  the  population  of  orbital  micrometeorites. 
This  increase,  coupled  with  the  cross-section  provided  by  larger  satellites  implies  that 
orbital  impacts  will  become  much  more  frequent.  Currently  most  damage  assessment  is 
provided  via  performance  degradation  analysis.  This  can  led  to  wasted  manpower  trying 
to  troubleshoot  a  device  that  was  literally  vaporized  at  impact  and  no  longer  exists.  More 
significantly,  the  introduction  of  large,  manned  space  stations  will  require  a  method  of 
damage  assessment  other  than  “go  look  and  see  what  has  been  hit.”  Rapid  and  accurate 
response  to  an  impact  could  mean  the  difference  between  life  and  death  for  the  astronauts 
onboard.  This  thesis  will  propose  a  method  of  rapid  impact  assessment  using  a 
distributed  sensor  system. 

These  larger,  lightweight  satellites  are  also  more  prone  to  low  firequency 


1 


vibrations  than  ever  before.  When  on  orbit,  the  operating  frequencies  of  equipment 
installed  on  the  same  structure  may  induce  performance  degradation  or  premature  failure 
of  other  satellite  components.  In  order  to  protect  sensitive  equipment  from  these 
vibrations,  there  are  two  possible  approaches.  These  options  are  to  dampen  the  vibration 

or  actively  suppress  the  effect  at  a  specific  location.  This  thesis  will  discuss  vibration 
suppression. 

Inherent  in  all  structures  is  a  degree  of  natural  passive  damping.  This  may  be 
enhanced  with  installed  devices,  such  as  visco-elastic  dampers.  These  devices  may  be 
large,  and  thus  not  worth  the  additional  lift  capacity  required  to  bring  them  to  orbit. 
Active  damping  is  another  option  available  for  vibration  reduction.  Implementing  an 
active  damping  system  can  be  difficult  due  to  the  difficulties  in  modehng  the  dynamic 
characteristics  of  the  structure.  However,  if  the  control  law  does  not  depend  on  the 
structure,  then  active  suppression  techniques  may  be  chosen. 

The  Naval  Postgraduate  School  (NPS)  Space  Truss  is  a  small-scale  flexible 
structure  that  has  two  piezoceramic  actuators  installed.  The  piezoceramic  actuators  were 
selected  because  of  their  high  bandwidth,  lightweight,  and  low  power  consumption. 
Previous  experiments  have  shown  that  a  computer-modeling  program  can  predict 
reductions  in  vibration  amplitude  using  a  detailed  dynamic  model  of  the  NPS  Space 
Truss  with  the  active  -damping  simulated  on  the  model.  The  simulated  results  of  the 
controller  developed  by  this  thesis  were  also  evaluated  with  this  model  and  verified  by 
experiments  on  the  NPS  Space  Truss. 
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The  NFS  Space  Trass  is  now  a  test  platform  for  further  experiments  on  active 
vibration  reduction  and  impact  analysis.  The  techniques  learned  using  the  NFS  Space 
Trass  could  be  applied  to  larger  structures,  such  as  the  International  Space  Station. 

B.  OBJECTIVES 

The  objective  of  this  thesis  were  as  follows: 

•  Develop  an  active  vibration  control  system  that  is  capable  of  reducing  single 
frequency  vibrations  at  a  remote  location  and  is  independent  of  structural 
characteristics. 

•  Evaluate  the  controller  using  an  active  FEM  and  compare  the  results  with  the 
experimental  response. 

•  Develop  a  method  to  determine  the  location,  force,  and  orientation  of  an 
impact  using  a  minimal  distributed  sensor  grid. 

C.  SCOPE  OF  THESIS 

The  scope  of  this  thesis  includes: 

•  Creation  of  a  control  algorithm  to  actively  suppress  an  axial  component  of  a 
constant  frequency  disturbance  at  a  random  point  from  a  remote  location. 

•  The  integration  of  the  dSFACE  digital  signal  processing  system  [Ref.  1]  as  the 
heart  of  both  the  signal  analysis  used  for  impact  analysis  and  the  active 
control  system  on  the  NFS  space  trass. 

•  Evaluation  of  the  active  controller  using  an  ANSYS  [Ref  2]  finite  element 
model  of  the  NFS  space  trass. 
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•  Comparison  of  the  FEM  results  with  the  actual  truss  by  experimentation. 

•  Proposal  of  a  method  to  determine  impact  characteristics  using  a  system  of 
distributed  sensors. 

The  use  of  ANSYS,  MATLAB  [Ref.  3],  and  dSPACE  were  instrumental  to  the 
work  contained  herein.  This  thesis  is  not  intended  to  replace  the  use  of  tutorials  available 
in  MATLAB,  ANSYS  or  dSPACE.  On-line  and  printed  guides  and  tutorials  are  available 
to  give  the  layman  an  opportunity  to  learn  the  system  in  order  to  use  the  information 
presented.  A  basic  outline  of  how  to  get  the  programs  and  systems  operating  is  included, 
and  has  sufficient  detail  to  allow  a  user  with  some  experience  to  run  the  desired 
programs.  Where  specific  program  instructions  are  given,  menu  commands  will  be  used 
with  a  “>”  to  indicate  a  sub-menu  selection,  or  the  commands  as  typed. 

D.  METHODOLOGY 

The  research  for  this  thesis  followed  several  different  paths  that  involved  signal 
processing,  active  control  systems,  FEM,  and  the  implementation  of  the  various  software 
packages  used.  The  first  goal  was  to  reduce  the  vibration  of  various  nodal  points  on  the 
NPS  Space  Truss,  with  the  ultimate  objective  of  total  stabilization.  Due  to  the  inherent 
limitations  of  vibration  damping  in  a  structure  with  multiple  transmission  paths,  the  idea 
of  active  suppression  became  the  primary  focus.  A  single  axis  approach  was  chosen  as 
the  first  iteration  on  the  system  for  the  sake  of  testing  and  simplicity. 

The  Adaptive  Multi-Layer  LMS  Controller  was  developed  one  piece  at  a  time. 
After  initial  experimental  testing  proved  that  the  Adaptive  Notch  Canceller  [Ref  4]  could 
be  applied  to  the  NPS  Space  Truss,  it  was  decided  to  include  a  fi’equency  identification 
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routine  that  would  make  the  controller  more  versatile.  After  the  identification  algorithm 
was  applied,  it  was  discovered  that  the  nodal  signal  could  not  be  used  for  identification 
due  to  the  ability  of  the  controller  to  suppress  the  signal  to  the  noise  level.  The 
identification  input  was  then  shifted  to  the  disturbance  generator,  vice  the  nodal  signal. 

Further  testing  quickly  revealed  that  the  LMS  gain  was  a  variable  based  on  the 
frequency  of  the  disturbance.  The  next  modification  added  the  ability  of  the  controller  to 
adjust  the  magnitude  of  the  gain  based  on  system  response.  This  was  only  partially 
successful  due  to  the  phase  shifts  experienced  over  the  frequency  range  of  interest.  The 
ability  of  the  controller  to  change  the  sign  of  the  gain  was  then  added  to  the  system.  The 
logic  for  the  adaptive  gain  adjustment,  translated  into  mathematical  a  format,  came  from 
literally  dozens  of  experimental  runs  attempting  to  determine  a  pattern  to  the  optimal  gain 
values.  The  controller  was  essentially  complete  at  this  point. 

After  initial  testing,  a  second  active  element  was  installed  in  the  NFS  Space  Truss 
in  order  to  evaluate  the  new  controllers  ability  to  adapt  regardless  of  actuator  position  and 
orientation.  This  testing  was  completely  successful  and  no  additional  adjustments  were 
required  based  on  which  element  was  chosen.  The  remaining  experimentation, 
evaluation,  and  tuning  was  dedicated  to  documenting  the  abilities  of  the  system  and 
adjusting  the  various  adaptation  rates  for  maximum  stability. 

The  next  phase  in  the  research  dealt  with  the  modification  of  a  FEM  created  for 
previous  research  [Ref.  5]  to  determine  the  ability  of  a  software  model  to  predict  the  new 
controllers  effectiveness.  With  the  assistance  of  Mr.  Sheldon  Imaoka,  a  Collaborative 
Solutions,  Inc.,  Engineering  Consultant,  the  original  ANSYS  structural  model  was 
modified  to  include  the  new  active  element  and  the  program  to  implement  the  controller 
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was  written.  It  was  determined  that  limiting  the  analytical  controller  implemented  in  the 
FEM  to  include  only  the  phase  adjustment  portion  of  the  controller  would  significantly 
simplify  the  programming  requirements  and  allow  for  a  valid  comparison  of  the  steady 
state  operations.  That  decision  also  dictated  multiple  runs  of  each  experiment  to 
determine  the  optimal  gain  values  required.  This  experience,  coupled  with  the  dozens  of 
experimental  experiments  run  with  the  same  goal,  was  valuable  in  determining 
qualitatively  the  ability  of  the  model  to  duplicate  actual  structural  responses  to  the  control 
system.  The  video  files  that  were  created  by  ANSYS  also  allowed  for  amplification  of 
the  displacements  and  resulted  in  a  greater  understanding  of  the  entire  structures  response 
to  suppressing  a  single  nodes  axial  vibration. 

The  finial  phase  of  the  research  was  dedicated  to  developing  a  method  of 
determining  impact  characteristics  based  only  on  sensory  input.  The  idea  of  a  neural 
network  was  explored  and  discarded  as  overly  complex.  The  concept  that  was  eventually 
decided  on  was  based  on  the  idea  that  the  some  level  of  accuracy  on  impact  location, 
force  and  energy  could  be  obtained  by  comparing  impact  signals  with  known  impact 
responses.  Signal  cross  correlations  were  used  to  compare  signals  and  initial  testing  was 
done.  The  testing  produced  results  that  were  generally  encouraging,  but  were  not 
extensive  enough  for  documentation  purposes.  This  phase  of  the  research  was 
discontinued  due  to  time  limitations  and  is  thus  referred  to  in  the  Future  Research  section. 

E.  THESIS  ORGANIZATION 

Chapter  II  contains  background  information,  theory  and  a  summary  of  the  work 
that  preceded  the  tasks  performed  by  the  author.  Chapter  El  contains  the  details  on  the 
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development  and  implementation  of  an  active  controller  for  the  NPS  space  truss,  using 
MATLAB  /  Simulink,  and  dSPACE.  Chapter  IV  describes  the  finite  element  model 
created  in  ANSYS.  Chapter  V  describes  the  experimental  methodology  used  to  verify  the 
controller  using  the  finite  element  model.  Chapter  VI  is  a  description  of  the  proposed 
method  of  signal  processing  required  for  impact  analysis  and  recommendations  for 
further  study  regarding  the  Adaptive  Multi-Layer  LMS  controller.  Finally,  Chapter  Vn 
provides  conclusions  in  the  areas  researched. 
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n. 


THE  NFS  SPACE  TRUSS 


A.  TRUSS  DESCRIPTION 

1.  Background 

The  NPS  space  truss  obtained  by  the  Spacecraft  Research  and  Design  Center 
(SRDC)  at  NPS  from  the  Naval  Research  Laboratory  (NRL)  and  is  as  exact  replica  of  the 
NRL  space  truss.  The  engineers  at  NRL  obtained  approximately  a  100-fold  reduction  in 
power  amplitude  due  to  the  active  control  system  using  the  controller  [Ref.  6]. 

To  date,  the  NPS  space  truss  has  had  modal  testing  performed,  and  two  active 
control  piezoelectric  struts  have  been  installed  and  tested.  Planned  research  using  the 
NPS  space  truss  includes  the  integration  of  fiber-optic  strain  gages. 

2.  Elements  and  Construction 

The  NPS  space  truss  is  3.67  meters  long,  0.67  meters  tall,  and  is  composed  of  52 
aluminum  nodes  joined  by  161  elements  in  a  cubic  12-bay  structure.  It  is  attached  at  the 
center  to  a  fixed  base  plate  and  mounted  on  a  Newport  Vibration  Control  System  Table. 
The  layout  of  the  unmodified  truss  is  shown  in  Figure  1. 


Figure  1.  The  NPS  Space  Truss  [From  Ref  5] 
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The  node  balls  are  precision  milled  from  7075-T6  Aluminum  and  a  standoff  is 
used  to  provide  an  interface  between  the  elements  and  the  node  ball.  The  elements  are 
5/16-inch  aluminum  alloy  tubing  with  0.035-inch  wall  thickness  and  have  two  lengths,  a 
shorter  length  for  the  battens  and  longerons  (100  total),  and  a  larger  length  for  the 
diagonal  elements  (61).  An  exploded  view  of  the  element  to  node  assembly  is  shown  in 
Figure  2.  A  more  detailed  description  of  the  truss  can  be  found  in  Reference  5  and  the 
assembly  procedure  used  is  Reference  7. 


Figure  2.  Node  Assembly  Details  [From  Ref  8] 

The  isolation  table  uses  a  compressed  Nitrogen  cushion  to  provide  for  high 
frequency  vibration  attenuation  (greater  than  99%  above  12  Hz)  [Ref  9].  Detailed  data 
on  the  NFS  space  truss  are  contained  in  Appendix  A. 

B.  ACTIVE  CONTROL  ELEMENT 
1.  Overview 

The  active  elements  installed  in  the  NFS  space  truss  consist  of  piezoelectric 
actuators  and  co-located  force  transducers.  These  components  are  mounted  to  steel  rods 
in  a  manner  that  is  compatible  with  the  truss  elements.  Together  with  the  nodal  3-axis 
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accelerometers  and  the  computer  data  processing  system,  these  components  comprise  the 
adaptive  suppression  system. 


2.  Piezoelectric  Theory 

Piezoelectric  (also  called  piezoceramic)  materials  elongate  when  a  voltage  is 
applied  in  a  pre-determined  direction.  Conversely,  when  a  deformation  is  applied  to  a 
piezoelectric  material  with  either  an  external  force  or  strain,  an  internal  electric  potential 
will  be  developed  that  can  be  measured.  In  this  manner,  piezoelectric  material  is  used  for 
both  a  sensor  and  as  an  actuator. 

Piezoelectric  materials  are  anisotropic  in  nature  with  orthogonal  directions  related 
to  the  poling  axis  mechanically  by  Poisson’s  ratio.  While  not  truly  linear,  and  containing 
some  hysteresis,  the  piezoelectric  effect  can  be  modeled  as  such,  with  minimal  error.  The 
piezoelectric  electroelastic  relations  can  be  described  with  the  following  equation  pief 
10],  shown  in  stiffiiess  method  form: 


W  Wlj^i 

[^r 


(2.1) 


The  stress  and  strain  vectors  represent  the  three  axial  and  three  shear  stresses  and  strains. 
Rewriting  equation  (2.1)  in  the  applied  force  method,  given  by  [Ref.  1 1]: 


[4 


(2.2) 


These  equations,  when  the  [d]  and  [e]  matrices  are  zero,  reduce  to  Hooke’s  Law 
and  the  dielectric  equation  [Ref  10,  p.  4].  The  relationship  between  the  [d]  and  [e] 
matrices  is  given  by: 


11 


[e]  =  [cnd]  (2.3) 

Where  [c]  is  a  symmetric,  (6  x  6)  matrix  that  contains  all  of  the  linear  elastic  and 
flexural  terms.  The  [e]  and  [d]  matrices  are  (6  x  3)  and  contain  the  piezoelectric  coupling 
terms.  The  [s]  matrix  is  a  diagonal  matrix  with  the  dielectric  permitivities.  The 
individual  constants  in  the  matrices  are  usually  indexed  by  their  coordinate  axes  (1,2  and 
3);  the  3-axis  is  the  poling  axis. 

3.  Application  and  Operating  Characteristics 

Piezoelectric  actuators  are  attractive  due  to  their  high  bandwidth,  compact  size, 
and  low  power  requirements.  Special  care  must  be  used  while  working  with  them 
because  even  though  the  power  required  by  a  piezoelectric  actuator  is  low,  the  voltages 
required  are  high. 

The  basic  schematics  for  a  piezoelectric  actuator  are  presented  in  Figure  3.  The 
actuators  are  mounted  such  that  the  desired  applied  force  is  in  the  direction  of  the  poling 
axis.  In  the  following  figure,  the  poling  axis,  applied  electric  field,  and  the  resulting 
displacement  are  shown  in  the  piezoelectric  crystal. 


Figure  3.  Piezoelectric  Actuator  Configmation 
Piezoelectric  materials  are  limited  to  an  applied  electric  field  of  2kV/mm  [Ref 
12].  To  fimctionally  eliminate  this  problem,  piezoelectric  material  are  stacked  commonly 
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in  the  direction  of  the  poling  axis,  with  interleaving  foil  electrodes  to  provide  a  greater 
cumulative  strain  and  thereby  greater  force  to  the  structure.  The  polarity  of  the 
piezoelectric  wafers  is  inverted  at  each  step  to  simplify  the  electrode  placement  and 
maintain  uniform  global  effect.  This  configuration  for  a  stacked  linear  actuator  is  shown 
in  Figure  4. 


Figure  4.  Stacked  Piezoelectric  Actuator  [From  Ref  13] 

The  active  piezoelectric  element  is  moimted  in  the  structure  such  that  the  poling 
axis  is  directly  in  line  with  the  element  that  it  is  replacing  to  allow  for  the  maximum 
piezoelectric  effect  to  be  felt  by  the  structure.  The  piezoelectric  active  control  element 
selected  was  the  Piezoelectric  Translation  Model  P-848-30,  built  by  Polytek  Physik 
Instrumente  of  Hamburg,  Germany  [Ref  5]  and  shown  in  Figure  5.  The  P-843.30  has  a 
maximum  operating  voltage  of  lOOV  and  is  cylindrical  in  shape,  with  a  14-mm  diameter 
and  73-nim  length.  Other  characteristics  are  listed  in  Appendix  A. 
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Figure  5.  Physik  Instnimente  P-843.30  [From  Ref.  5] 

The  second  connector  shown  in  Figure  5  is  for  an  integrated  metal  foil  strain  gage 
that  was  not  used  for  this  research.  The  P-843.30  has  a  rated  open  loop  travel 
characteristic  of  45  pm/lOOV  +/-  20-percent.  Verification  of  this  was  performed  prior  to 

installation  [Ref.  14,  pp.  25-28].  The  results  fi'om  this  measurement  are  presented  in 
Figiue  6. 


Expansion  Characteristics  of  Model  P-843.30  Piezo 


Volts  (V) 

Figure  6.  Piezo  Model  P-843.30  Expansion  Characteristics  [From  Ref  14] 


Ab  average  of  50  lam/lOOV  was  assumed  for  further  calculations. 


4.  Locating  The  Active  Element 

The  optimal  active  strut  location  was  based  on  the  element  that  experienced  the 
maximum  strain  energy  in  the  first  mode.  A  finite  element  model  of  the  NFS  space  truss 
constructed  in  MATLAB  and  the  strain  energy  was  calculated  based  on  the  deformation 
of  the  truss  [Ref  5].  As  noted  in  the  analysis,  there  is  an  odd  symmetry  to  the  locations 
of  the  maximum  strain  energy  due  to  the  diagonal  elements  in  the  truss. 

Based  on  the  results  of  this  analysis,  Element  101  was  replaced  by  the  first  active 
truss.  The  second  active  element  was  placed  base  on  the  symmetry  of  the  analysis  and 
replaced  Element  107.  Operational  experience  with  the  current  controller  suggests  that 
this  may  not  have  been  the  best  choice.  The  abihty  of  the  user  to  select  the  control  node 
is  limited  by  the  physical  nature  of  the  structure.  If  the  active  elements  were  placed  in 
geometrically  umelated  positions,  the  likelihood  of  both  of  them  being  ineffective  at  a 
give  frequency,  relative  to  the  same  control  node  location,  is  greatly  reduced.  The 
control  element  location  was  not  altered  after  initial  placement  due  to  the  cost  of  the 
piezoelectric  elements  and  their  susceptibihty  to  failure  due  to  applied  torque  diuing  the 
instillation  process. 

The  program  used  to  determine  the  mode  shapes  [Ref  5]  was  verified  and  slightly 
altered  and  is  included  in  Appendix  B.  The  first  ten  natural  frequencies  computed  by  the 
MATLAB  model  are  shown  in  Table  1,  along  with  the  experimental  results  firom  modal 
testing  [Ref  14]. 
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Table  1.  NFS  Space  Trass  Modal  Frequencies  [Ref.  14] 

The  first  four  mode  shapes  for  the  NPS  Space  Trass  were  obtained  by  taking  the 
eigenvector  displacements  and  adding  them  to  their  respective  nodal  coordinates.  The 
mode  shapes  generated  by  this  model  are  included  in  Appendix  C.  It  should  be  noted  that 
these  analysis  were  conducted  on  a  bare  trass  model  (without  the  additional  mass  and 
stiffiiess  of  the  LP ACT  and  the  active  elements).  The  ANSYS  analysis  was  done  on  a 
modified  trass  model  that  includes  all  the  systems  installed  in  the  truss  as  used  in  the 
experimental  system.  The  mode  shapes  from  this  analysis  are  also  included  in  Appendix 
C  and  are  noticeably  different  from  the  bare  trass  model.  It  is  recommended  that  a 
modified  model  be  used  if  further  verification  of  the  optimal  active  element  location  is 
desired. 

5.  Force  Transducer 

Previous  research  using  the  NPS  Space  Trass  used  an  installed  force  transducer 
that  was  collocated  with  the  active  elements.  The  original  PCB  Piezotronics  Model 
208B02  General  Purpose  ICP  Force  Sensor  is  a  still  present  in  the  experimental  setup  and 
a  second  was  installed  with  the  additional  active  element  [Ref.  15].  Their  physical 
properties  were  included  in  the  ANSYS  FEM,  but  neither  output  was  used  in  this 
research.  The  model  208B02’s  operating  characteristics  are  included  in  Appendix  A. 
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6.  Installation  of  the  Active  Element 


The  second  active  element  was  constructed  with  identical  components  and 
procedures  as  the  first  element  [Ref.  14].  The  assembled  active  element  is  shown  in 
Figure  7. 


Figure  7.  Assembled  Active  Control  Element  [From  Ref  14] 


C.  LINEAR  PROOF  MASS  ACTUATOR  (LPACT) 

The  disturbance  force  is  provided  by  a  Linear  Proof  Mass  Actuator  (LPACT) 

installed  at  the  end  of  the  NPS  space  truss.  The  LPACT  is  a  model  CML-030-020-1 
manufactured  by  Planning  Systems,  Incorporated,  from  Melbourne,  FL  and  is  powered 
by  a  separate  amplifier  and  controller  assembly  with  embedded  feedback  electronics. 
The  LPACT  has  mounted  accelerometers  for  use  in  monitoring  or  driving  the  feedback 
loops  that  were  not  used  in  this  research,  but  could  be  used  as  a  disturbance  reference 
signal  for  frequency  identification.  The  LPACT  is  shown  in  Figure  8  and  further  details 
are  included  in  Appendix  A. 
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Figure  8.  LPACT  [From  Ref.  16] 

A  picture  of  the  installed  LPACT  is  provided  as  Figure  9. 


The  LPACT  has  its  own  frequency  characteristics  shown  in  Figure  10.  The 
natural  frequency  at  about  8.5-Hz  is  below  the  range  of  10  to  30  Hz  tested  in  these 
experiments.  The  mass  of  the  LPACT  does  contribute  to  a  shift  in  natural  frequencies  of 
the  truss  and  is  modeled  in  ANSYS.  The  addition  of  a  new  natural  frequency  in  the  8-10 
Hz  range,  that  is  the  product  of  the  LPACT  spring,  is  not  included  in  either  FEM. 
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frequency  (Hz) 


Figure  10.  LPACT  Transfer  Function  (feedback  loops  off)  [From  Ref  16] 

D.  LASER  DIODE  ASSEMBLY 

The  laser  diode  attached  to  the  truss  [Ref  14]  was  relocated  to  node  26  of  the 
truss.  The  laser  was  directed  to  the  far  wall.  The  only  use  of  this  assembly  was  to 
effectively  amplify  the  truss  vibrational  motion  and  thus  serve  as  a  visual  demonstration 
of  the  controller  in  action.  The  mounting  consisted  of  a  thin  aluminum  rod,  with  a  mass 
at  the  end,  housing  the  laser.  The  laser  diode  selected  is  a  1-mW,  635-mn  Model 
PLC6351FW  from  Laseimate  Corporation  of  Walnut,  CA  powered  by  a  Hewlett  Packard 
E3615A  DC  Power  Supply,  operating  at  2-4  Volts  nominal.  A  picture  of  the  laser  diode 
assembly  is  included  as  Figure  1 1  and  further  details  are  can  be  found  in  Reference  14. 
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Figure  1 1.  Laser  Diode  Assembly  [From  Ref.  5] 

E.  PREVIOUS  EXPERIMENTS 

Three  previous  theses  have  been  completed  on  the  NPS  space  truss:  One  by 
CAPT  Brent  K.  Andberg,  USMC,  an  other  jointly  by  LT  Scott  Johnson  and  LT  John 
Vlattas,  USN,  and  the  third  by  LT  Carey  M.  Pantling,  USN. 


1.  Andberg  [Ref.  17] 

In  his  thesis,  CAPT  Andberg  developed  an  FEM  of  the  NPS  space  truss  using  an 
NRL  code  entitled  NRLFEMI.  He  then  performed  modal  testing  on  the  truss  to  confirm 
the  model.  The  experimental  data  was  lacking  in  that  it  failed  to  observe  the  first  mode. 
Finally,  CAPT  Andberg  performed  a  technology  demonstration  of  the  use  of  Fiber-optic 
Bragg  Gratings  (FBG),  used  in  this  example  to  detect  the  motion  of  a  simple  cantilever 

beam.  In  the  future,  FBGs  will  be  installed  on  the  NPS  space  truss  for  shape 
determination. 
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2.  Vlattas  and  Johnson  [Ref.  14] 


This  Master’s  Thesis  focused  on  two  areas,  the  re-performance  of  modal  testing 
with  the  HP  35665A  two-channel  spectrum  analyzer  and  integration  and  testing  of  an 
active  control  device.  They  recommended  that  modal  testing  be  again  done,  to  overcome 
some  of  the  limitations  of  the  dSPACE  system  that  was  installed  at  the  time.  As  will  be 
discussed  in  the  next  chapter,  a  new  dSPACE  system  was  installed  on  the  NPS  space 
trass,  and  would  be  available  for  this  purpose,  if  desired. 

The  active  controls  integration  produced  good  results,  with  a  maximum  reported 
reduction  of  14.817  dB  at  16.85  Hz.  They  reported  using  a  disturbance  amplitude  of  100 
mV  for  the  LPACT  source.  It  was  determine  early  on  through  the  cinrent  course  of  this 
research  that  this  did  not  even  provide  a  sensor  signal  sufficient  to  overcome  system 
noise.  Therefore,  the  14.817-dB  reduction  in  amplitude  is  held  in  question. 

3.  Pantling  [Ref.  5] 

This  Master’s  Thesis  concentrated  on  the  analytical  modeling  of  the  NPS  Space 
Trass  while  active  control  methods  were  being  applied.  Models  were  created  using  both 
MATLAB/Simulink  and  ANSYS  software  packages.  An  BDIFF  control  law  was  used  to 
realize  a  15-20  dB  reduction  in  vibration  at  the  force  transducer  as  observed  using  the 
dSPACE  data  acquisition  and  processing  system.  This  was  compared  to  18-22  dB 
reductions  as  predicted  by  an  ANSYS  actively  controlled  finite  element  model.  This  was 
used  to  demonstrate  the  validity  of  predicting  the  effectiveness  of  a  control  authority 
using  a  FEM. 
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The  use  of  the  IDIFF  controller  did  reduce  the  magnitude  of  the  examined  degree 
of  freedom,  but  could  result  in  large  increases  in  other  DOF’s  that  were  not  used  in  the 
control  law.  He  recommended  the  development  of  other  control  laws  and  the 
implementation  of  additional  active  elements. 
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III. 


CONTROLLER  DESCRIPTION  AND  INTEGRATION 


A.  OVERVIEW 

The  active  vibration  control  research  done  using  the  NPS  Space  Truss  has 
primarily  focused  on  minimizing  the  effect  of  equipment  induced  vibration  on  other 
components  installed  on  the  trass.  Attempts  have  been  made  to  reduce  the  vibration 
transmitted  via  a  single  truss  element  with  the  intent  of  reducing  the  total  energy  of  the 
structure.  The  effectiveness  of  that  approach  is  limited  by  the  percent  of  vibrational 
energy  that  is  transmitted  by  the  control  element.  The  control  authority  of  the  actuator 
becomes  the  limiting  factor  if  energy  is  introduced  into  the  structure  that  cancels  out  the 
disturbance  effect  at  a  specific  node.  As  stated  in  the  OBJECTIVES  of  this  research,  the 
control  system  developed  in  this  thesis  reduces  single  firequency  vibrations  at  a  remote 
location  and  is  independent  of  structural  characteristics.  The  control  law  is  also  flexible 
enough  to  allow  the  user  to  choose  the  node  on  the  trass  that  they  want  to  control,  the 
active  element  they  want  to  use,  and  the  axis  of  concern.  This  allows  for  a  specific  focus 
for  the  purpose  of  the  active  element  and  significantly  improved  results. 

The  last  active  vibration  control  research  on  this  structure  was  done  using  an 
integral  plus  double  integral  force  feedback  (IDIFF)  control  system  that  was  incorporated 
by  Vlattas  and  Johnson  [Ref  14]  and  evaluated  by  Pantling  [Ref  5].  The  effect  of  this 
controller  was  a  reduction  of  15-20  dB  in  the  vibration  force  as  measured  directly 
adjacent  to  the  active  control  element.  This  method  was  replaced  with  a  newly 
developed,  adaptive  multi-layer  LMS  controller  that  utilizes  the  same  hardware 
configuration  as  previous  work,  but  is  able  to  reduce  the  vibration  acceleration  at  a 
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remote  sensor  location.  This  system  consistently  reduces  single  frequency  vibrations  to 
the  level  of  system  noise  at  any  one  of  several  key  nodes,  using  either  of  two  installed 
active  control  elements.  Results  documented  in  this  research  show  reductions  in  excess 
ofSOdB. 

Figure  12  is  an  overview  of  the  entire  system,  showing  the  positions  of  the 
disturbance  source,  nodal  sensors,  and  active  elements.  It  also  includes  the  generalized 
flow  path  of  the  control  signal  and  the  processing  elements,  but  is  not  intended  to  be  an 
all-inclusive  description.  Detailed  descriptions  of  the  mathematics  and  experimental 
application  are  outline  in  the  rest  of  this  section. 

The  adaptive  multi-layer  LMS  controller  was  developed  using 
IVLA.TLAB/SIMULINK  software  and  implemented  using  dSPACE  RTI  on  a  dSPACE 
DS1103  PPC  Controller  Board.  The  DS1103  has  a  Motorola  PowerPC  604e 
microprocessor  for  a  central  processing  unit  (CPU),  and  resides  in  a  triple-wide  ISA  slot 
in  a  host  PC  [Ref.  1,  p.  12].  The  DS  1103  board  uses  connectors  for  an  external 
input/output  (I/O)  box  that  contains  BNC  fittings  and  standard  computer  cable 
connections  for  analog  to  digital  (ADC),  digital  to  analog  (DAC)  and  support  for  other 
cabling  formats  (e.g.  RS-232).  Control  of  the  dSPACE  CPU  and  access  to  its  memoiy 
(128MB)  is  done  with  the  use  of  the  dSPACE  Control  Desk.  Specifics  of  each  sofh^  arc 
component  are  discussed  in  grater  detail  in  the  following  sections. 
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Figure  12.  Control  System  Overview 


B.  ADAPTIVE  MULTI-LAYER  LMS  CONTROLLER 


1.  Adaptive  Phase  Adjustment 

In  1998,  Bertran  and  Montoro  proposed  an  Adaptive  Notch  Canceller  [Ref.  4]. 
Their  work  originated  from  the  need  to  suppress  vibrations  induced  by  rotating 
machinery.  This  controller  needed  only  mild  structural  assumptions  (stable,  linear,  single 
input,  single  output  plant,  with  tonal  disturbances),  which  are  consistent  with  the  NPS 
Space  Truss.  This  control  algorithm  was  used  as  the  foundation  of  the  adaptive  controller 
developed  in  this  thesis.  The  advantages  of  this  algorithm  include  its  computational 
simplicity,  neutral  response  to  uncontrolled  frequencies,  and  the  possibility  to  operate 
multiple  controllers  in  parallel  at  different  frequencies.  A  diagram  block  describing  this 
controller  is  shown  in  Figure  13. 


Figure  13.  Adaptive  Notch  Canceller  [Ref  4] 


Each  of  the  LMS  blocks  implement  an  adaptive  finite  impulse  response  (FIR) 
filter  of  length  one,  using  the  stochastic  gradient  algorithm  known  as  the  Least  Mean- 

Square  (LMS)  algorithm.  The  mathematical  relationships  of  the  signals  shown  in  Figure 
13  are: 
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^  ’  (3.1) 

CO  =  Digital  Frequency  [0,2 /r] 

x„  -  a„  cos[am)  +  b„  sin(^y«) 

n=  Time  Index  (3.2) 

neN 

y„  =  a„acos{(m  +J^)  +  b„asm(on  +  0)  (3.3) 

en=yn+d„  (3.4) 

««+i=^„+A^„cos(ty«)  (3.5) 

^,+i=^„+A^„sin(6w)  (3.6) 

Assuming  that  the  disturbance  is  a  constant  sinusoidal  signal  in  the  form: 

(3.7) 


For  a  convergence  analysis,  assuming  the  control  signal  rewritten  as: 

x„  =  4,e>'' 

where,  (3.8) 

A  =a„+jb„ 

The  corresponding  system  output  is  then: 

(3.9) 

Substituting  equation  (3.9)  into  the  error  equation  (3.4),  you  get: 

=  (3.10) 

Equations  (3.5)  and  (3.6)  can  be  converted  into  complex  form  and  combined: 

4„,=4,+//e„e-^""’  (3.11) 


27 


By  combining  equations  (3.10)  and  (3.11),  the  iterative  equation  for  the 
exponential  gain  term  is  obtained: 

4.i=(l  +  A««^'’)4+</'e'-'“  (3.12) 

Now,  substituting  equation  (3.7)  into  the  above  equation  yields: 

4+1  =  (l  +  Mae^'^)4,  (3.13) 

Because  the  disturbance  and  the  plant  variables  (a  and  P)  are  assumed  constant  at 
a  given  frequency,  the  equation  converges  if: 

|l  +  //^3re-^^|<l  (3.14) 

If  the  proper  value  of  (x  is  chosen,  the  value  of  will  approach  a  constant  as  the 
time  index  (n)  approaches  infinity.  Therefore  equation  (3.11)  can  be  written  as: 

A.,  =  4 +/'«.=■*’ 

Ae,e-^"=4-4., 

jue^  =  0  (as  M  — >  oo) 

lime„=0  (3.15) 

One  of  the  main  advantages  of  this  algorithm  is  computational  efficiency.  It 
requires  only  seven  floating-point  operations  per  iteration  (two  addition,  five  products, 
and  two  trigonometric  evaluations). 
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2.  Adaptive  Gain  Adjustment 


The  convergence  analysis  of  the  adaptive  phase  adjustment  algorithm  requires 
that  the  correct  value  of  the  imit-less  adaptation  constant  (p)  be  determined.  In  order  to 
make  the  controller  more  versatile,  a  method  of  determining  this  value  automatically  was 
developed  and  implemented  as  a  gain  factor  applied  to  the  error  signal. 


Figure  14.  External  Gain  applied  to  LMS  Controller 


By  introducing  the  variable  gain  (G)  to  the  error  signal,  the  system  is  modified  as 
shown  in  Figure  14.  By  defining  an  adaptation  constant  internal  to  the  LMS  algorithm  as 
p',  then  equations  (3.5)  and  (3.6)  become: 

cos  [m)  (3.16) 

Kx  =K+  sin  (fiw)  (3.17) 

If  we  redefine  the  effective  adaptation  constant  (p)  as  the  product  of  an  internal 
adaptation  constant  (p')  and  an  external  gain  (G),  then  the  original  LMS  equations  and 
convergence  criterion  are  xmchanged.  This  configuration  allows  the  stability  and  rate  of 
convergence  of  the  phase  adjustment  algorithm  to  be  controlled  independently  of  p'. 

//  =  //'G  (3.18) 
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The  initial  attempts  to  initiate  a  third  LMS  to  control  the  gain  directly  were  not 
successful  due  to  complex  and  vaiymg  relationships  between  adaptation  rates,  structural 
phase  effects,  and  the  effects  of  various  error  and  gain  magnitudes.  It  was  decided  to 
split  the  adaptive  gain  adjustment  into  two  distinctive  parts.  The  first  part  deals 
specifically  with  the  required  sign  of  the  gain  to  be  applied  to  the  phase  adjust  error 
signal.  The  second  algorithm  adjusts  the  magmtude  of  the  gain  to  improve  convergence 
rate,  thus  provide  a  faster  response. 

The  sign  correction  factor  is  verified  every  second,  based  on  a  comparison  of  an 
instantaneous  error  amplitude  compared  to  the  largest  amplitude  recorded.  If  the  current 
amplitude  is  larger  than  the  previous  maximum,  the  sign  of  G  is  reversed  and  the 
maximum  is  updated.  The  signal  used  to  represent  the  error  amphtude  is  the  raw  error 
signals  standard  deviation,  which  is  directly  proportional  to  the  maximum  amplitude  of  a 

constant  fi-equency  signal.  The  implementation  of  this  logic  statement  is  shown  in  detail 
in  section  [III.C.  1  .e]. 

A  third  LMS  algorithm  was  used  to  adjust  the  magnitude  of  the  gain  (G)  as  shown 
in  Figure  15.  As  previously  noted,  the  adaptation  rates  of  the  various  sections  of  the 
controller  are  critical.  The  initial  magmtude  of  G  must  be  large  enough  to  allou  for  the 
correct  sign  to  be  determined  prior  to  significant  variation  in  the  magnitude.  If  the  gain 
LMS  adjustments  the  magmtude  of  the  gain  through  zero  (firom  positive  to  negative),  the 
effect  of  the  magnitude  adjustment  counters  the  sign  adjustment.  The  two  algorithms 
would  then  become  competitive  as  the  gain  magnitude  oscillates  around  a  zero  mean. 
The  phase  adjustment  error  signal  would  not  necessarily  have  a  zero  mean  due  to 
frequency  variations,  but  the  magnitude  of  the  error  signal  would  be  reduced  to  typical 
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experimental  system  noise  levels,  thus  effectively  stopping  controller  adaptation 
regardless  of  the  of  the  actual  error.  With  the  proper  adaptation  rates  applied  to  the 
system,  the  controller  will  therefore  converge  regardless  of  the  value  of  p.  That  does  not 
assure  controller  effectiveness,  only  convergence,  because  the  gain  corrected  error  (G  en) 
is  artificially  being  forced  to  zero  as  the  gain  LMS  removes  the  system  input.  The  net 
effect  of  this  condition  would  be  to  transform  the  controller  into  another  disturbance 
source  that  may,  or  may  not  suppress  the  original  disturbance. 


Normalized 

i^plirude 


Jbiplitude 

Trend 


Figure  15.  Gain  Amplitude  Adjust  LMS 

In  Figure  21,  the  “Normalized  Amplitude”  is  used  to  establish  an  error  magnitude 
and  is  the  instantaneous  amplitude  used  for  to  determine  the  gain  sign,  divided  by  the 
maximum  value.  As  the  name  suggests,  this  bounds  the  error  magnitude  to  1.  The 
“Amplitude  Trend”  is  set  to  either  +/-  1  based  on  the  approximate  derivative  of  the  error 
amplitude.  This  combination  results  in  a  range  of  errors  of  -1  <  e’n  <1.  This  results  in 
the  following  system  of  equations: 


(3.19) 

|GL=g.(l)  =  #. 

(3.20) 

Sn-^\  Mcain^n 

(3.21) 
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This  results  in  a  maximum  step  size  per  iteration  of  the  gain  to  be  applied  to  the 
controller  phase  adjustment  error  signal  of: 

N«.i-|^LNA(7a«  (3.22) 

This  provides  a  predictable  foundation  upon  which  to  base  the  response  time  of 
the  other  LMS  algorithms  as  described  in  section  [in.B.  4], 

3.  Frequency  Identification 

In  order  to  make  this  controller  as  robust  as  possible,  it  is  necessary  to  identify  the 
disturbance  frequency  and  generate  the  appropriate  sine  and  cosine  functions.  This  also 
supports  the  assumption  in  the  convergence  analysis  that  the  inputs  to  the  phase  adjusting 
LMS  algorithms  are  sine  and  cosine  functions  of  the  actual  disturbance  frequency.  Only 
the  accuracy  of  the  method  of  frequency  identification  is  relevant  to  the  operation  of  the 
controller,  not  the  method.  For  a  proof  of  concept,  the  identification  process  used  here 
assumes  that  there  is  a  single  frequency  disturbance  source  and  updates  the  calculated 
frequency  every  1.5  seconds.  It  is  assumed  that  a  signal  would  be  available  from  the 
disturbance  somce  that  could  be  used  for  identification  purposes.  In  this  case,  the 

disturbance  frequency  identification  is  done  using  the  signal  generator  output  to  the 
LPACT. 

The  disturbance  is  assumed  to  satisfy  equation  (3.23). 

=  Acos{a}n  +  a) 


It  follows  that  dn  also  satisfies  the  linear  difference  equation  (3.24). 


(3.23) 


(3.24) 


with, 

fl,  =  -2cos(ty) 

^2  =  1 

Rewriting  equation  (3.24)  as  a  set  of  linear  equations, 

<=-K-,  ^2]“' 


(3.25) 


The  solution  requires  only  two  samples,  but  in  order  to  account  for  system  noise, 
N  samples  are  taken  and  the  system  is  written  as. 


or 

[d]  =  [D][a\  (3.26) 

Using  the  least  mean  squared  method  to  solve  this  system  yields  the  best  fit  for 
the  values  of  [a], 

[a]  =  ([Z)]^  [/)])■' [Z)nrf]  (3.27) 

Defining  the  following  function: 

1/  \  ,  »  -1  t  -2  +  o.'.z  +  ai 

A(z)  =  \  +  a^z  +a2Z  - - \ - -  (3.28) 

z 

It  can  be  shown  that  an  equivalent  form  of  the  numerator  is. 
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(3.29) 


+  a[z  +  a'  =  (z  -  e^‘^)(z  -  =  z^  -  (6^“^+  e‘^'®)z  + 

flr;  =  -(e^<-+e->) 

fl' =(e-'"'e--'‘^)  =  l 

Using  trigonometric  equalities,  it  is  shown  that  equation  (3.29)  also  satisfies 
equation  (3.24). 


al  -  e  =  -2cos(^y)  = 

fl'=(e^"e-^")  =  l  =  a2 


(3.30) 


Substituting  the  original  values  of  ai  and  a2  into  the  numerator  of  equation  (3.28), 


z  +  a,z  +  =  0 

Solving  for  the  roots  of  the  numerator, 

z^  +a^z  +  a2  =  0 
^_-fli±Vaf-4a^ 


-a,  J-af 

j  = — !-+-2 — L 

f  r\  / 


+  4iz, 


Z2  = 


2  2 


2  2 

Using  the  first  root  to  solve  for  the  fi’equency, 
-«1  ^  +4^2 


2  2 
therfore. 


•  =  4e  +  74 


7m 


Ct)  = 


\tan 


f  -  —{Sample  Frequency)  [Frequency  (Hz)] 


-1 


\^ReJ 


[Digital  Frequency  (Rad)] 


2fr 


(3.31) 


(3.32) 


(3.33) 


The  implantation  of  this  calculation  is  shown  in  section  [m.C.  1  .a] . 
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4.  Component  Integration 


The  interactions  between  the  various  parts  of  an  adaptive  multi-layer  LMS 
controller  are  dictated  by  their  individual  adaptation  rates.  For  the  LMS  algorithm,  the 
update  rate  is  established  by  the  system  sample  frequency  (V^/)  and  the  step  size  as 
determined  by  the  following  equation. 


Aw  =  w„  -  =  Step  Size 

where, 

//  =  LMS  Adaptation  Constant 
=  LMS  Input  Signal 
=  LMS  Error  Signal 


(3.34) 


This  makes  the  each  LMS  algorithm  have  an  adaptation  rate  of, 

D  _/5£l 
^IMS  ~  .  ^ 

At 


(3.35) 


The  exact  value  of  each  step  size  is  a  function  of  the  LMS  error  signal.  As  the 
controller  minimizes  the  error,  the  step  size  is  reduced  and  the  reaction  rate  slows  down 
to  reach  equilibrium  at  zero  error. 

For  the  gain  sign  correction  algorithm,  the  correction  is  applied  at  predetermined 
intervals  (zJ7).  The  size  of  the  correction  is  determined  by  the  magnitude  of  the  initial 
gain  (G„).  This  makes  an  effective  adaptation  rate. 


n 


(3.36) 


The  initial  adaptation  rates  are  estimated  based  on  the  observed  orders  of 
magnitude  of  the  error  signals.  The  orders  of  magnitude  used  in  for  this  research  are 
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listed  below.  The  units  on  the  error  (e)  and  input  (u)  signals  are  irrelevant  as  long  as  the 
cancel  out. 


^Phase  ~  ^ 

/^Gain  =  1^10'^ 

^ Phase  ~  ^ 

G  =  l;cl0' 

^Gain  ~  ^ 

^Phase  =  1^1 0"' 

AT  =  1.5  sec 

s- 

II 

^tiMs  - 1^10'^  sec 

At  LMS  =  1^10  ^  sec 

The  initial  adaptation  rates  are  estimated  as  follows, 


R  _  2G  _  2(1x10^)  _ 


'■Sign 


AT  1.5 


0(100) 


sec 


(3.37) 


^Phase 


^//Gue']  (1x10-')(1x10')(1)(M0-M 


^Gain  ” 


juue 

At 


y  Gain 


MM=o(.)seo- 

1x10^  ^  ^ 


(3.39) 


Both  the  sine  and  cosine  LMS  algorithms  have  the  same  adaptation  rates,  labeled 
here  as  the  “Phase”  rate.  One  of  the  key  relationships  is  the  relative  dominance  of  the 
gain  LMS  algorithm  compared  to  the  phase  LMS  with  respect  to  reaction  time.  This  is 
required  because  the  magmtude  of  G  controls  the  convergence  rate  of  the  phase 
algorithm.  If  the  convergence  rates  were  on  the  same  order  of  magnitude,  the  effect  of 
one  could  conceivable  cancel  out  the  other.  If  the  phase  algorithm  is  made  dominant, 
then  the  benefit  of  the  gain  LMS  would  be  significantly  reduced  or  eliminated.  The  other 
significant  relationship  is  the  relatively  rapid  initial  assessment  of  the  sign  to  be  applied 
to  the  phase  error  signal.  The  order  of  magmtude  difference  in  reaction  rates  is  only  one 
part  of  this  relationship.  The  other  factor  that  must  be  addressed  is  the  time  allowed  for 
the  controller  to  work  between  error  magmtude  assessments  (as  explained  in  section 


36 


[in.B.  2]).  In  this  case,  a  reasonable  reaction  time  was  set  at  T  =  1  sec.  This  time  is  a 
function  of  the  adaptation  rates  of  the  LMS  algorithms.  If  the  controller  does  not  reduce 
the  magnitude  of  the  disturbance  in  the  time  allotted  (1  sec),  then  the  error  signal  sign  is 
reversed.  If  T  is  increased,  the  possibility  of  system  fluctuation  will  be  reduced  at  the 
expense  of  convergence  rate  and  in  extreme  cases  could  significantly  increase  the 
structure  vibration  before  effectively  canceling  it  out. 

The  values  listed  in  this  section  were  used  for  all  the  data  runs  and  dozens  of  other 
combinations  of  frequency,  active  element,  and  control  node  location.  All  results 
demonstrated  stability  and  effectiveness  with  varying  convergence  times  required.  While 
performing  demonstrations  after  the  conclusion  of  this  research,  an  unstable  combination 
was  foimd.  The  combination  of  frequency,  element,  node,  and  axis  are  summarized  at  the 
end  of  the  experimental  results  [V.C.  11].  This  discovery  illustrates  that  the  relationship 
demonstrated  here  is  vestal,  but  not  universal.  Further  experimentation  is  required  to 
either  correct  this  condition  or  fully  document  the  instability. 

C.  SIMULINK/RTI 

1.  Simulink  Controller 

a.  Overview 

The  first  step  in  building  the  controller  was  to  assemble  a  Simulink  model. 
The  Simulink  model  contains  the  processing  connections  necessary,  implemented  with 
block  diagrams,  to  create  the  control  law.  The  complete  model  is  shown  in  Figure  16. 
Individual  functional  blocks  are  broken  down  and  explained  in  the  following  sections. 


37 


LPACT 


Disturbance 

Oenerator 


Disturbar^ce  Signal 


Frequency  Identification 


Reference  Signal 


Freq  (Hz) 


Controller 


W 


Nodal  Signal 


Nodal  Accelerations 


Nodal  Input 
Reset  Signal 


Control  Signal  Output 


Adaptive  Gain  Adjust 


Acf/ve  £}ement 

#f 


Reference  Signal 

Control  Signal  - 

Bror  Signal 

Reset  Signal  -j 

Piezo  #1 


Control  Signal 


Active  Element  Selection 


DAC 


DAC 


Active  Biemeni 
M2 


DAC 


Digits!  Signs} 
Ansiyzer 


Figure  16.  Adaptive  Multi-Layer  LMS  Controller  Overview 
The  nodal  acceleration  selected  as  the  error  signal  for  the  controller  is  also 
sent  directly  to  a  DAC  that  was  used  to  analyze  the  effectiveness  of  the  controller  in  near 
real  time.  This  allowed  for  analysis  independent  of  the  control  law  and  without 
competing  for  system  resources. 

The  dSPACE  system  relies  upon  the  Simulink  model  for  the  selection  of 
its  sampling  frequency.  The  sample  input  and  signal  output  rates  are  both  set  at  IkHz, 
but  the  signal  internal  to  the  controller  is  often  down  sampled  and/or  buffered  as 
explained  in  the  following  sections.  This  results  in  a  multiple  internal  sample  times 
ranging  from  0.001  to  2  seconds.  After  the  Simulink  model  was  complete,  the  command 


‘Tools>  RTW  Build”  created  a  real-time  program  and  object  file  that  was  executed  by 
the  CPU  of  the  dSPACE  system. 


b.  Disturbance  Generator 

The  method  of  generating  of  the  disturbance  signal  was  designed  to  be 
adjustable  from  the  dSPACE  ControlDesk.  Both  the  frequency  and  amplitude  of  the 
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disturbance  can  be  adjusted  using  by  modifying  the  “Disturbance  Freq  (Hz)”  and  the 
“Amplitude”  values.  The  ability  to  turn  off  the  disturbance  was  included  as  a  separate 
gain  factor.  The  default  values  are  shown  in  Figure  17. 


Disturbance  Freq 


(Hz) 


Clock  Hold 


Figure  17.  Disturbance  Signal  Generator 

The  signal  is  generated  by, 

D  =  Asm{2^6}t) 

A  =  Amplitude 

CO  =  Disturbance  Freq  (Hz) 

t  s  Clock 


(3.40) 


Even  though  the  clock  signal  is  digitized  to  a  1  kHz  sample  rate  for 
calculational  purposes,  the  finial  disturbance  signal  requires  it’s  own  1  kHz  hold  in  order 
to  assure  the  correct  synchronization  with  the  other  elements  of  the  controller. 

It  was  also  noted  that  the  clock  signal  would  degrade  over  time,  which 
compromised  the  reference  signal  to  the  phase  adjust  algorithm  and  the  overall 
effectiveness  of  the  controller.  This  was  only  observed  during  extended  experimental 
runs  and  is  a  documented  deficiency  in  the  software  package  used  for  controller 
implementation  [Ref  18,  p.  190]. 


c.  Frequency  Identification 

This  was  initially  constructed  to  use  a  reference  signal  from  a  sensor,  not 
the  digital  signal  from  the  disturbance  generator.  This  requirement  for  additional  signal 
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processing  motivated  both  the  input  filter  and  the  down  sampling.  This  established  the 
highest  effective  frequency  of  the  controller  at  33  Hz.  After  experimental  data  reveled 
the  now  obvious  fact  that  as  the  controller  reduced  the  disturbance  signal  by  30-50  dB, 
the  signal  to  noise  ratio  for  the  frequency  identification  process  increased  to  the  point  of 
ineffectiveness  of  the  algorithm.  This  generated  a  significant  error  in  the  reference  signal 
used  for  the  phase  adjustment  and  effectively  prevented  the  controller  from  converging. 


Although  the  reference  signal  was  switched  to  the  output  of  the  signal 
generator,  the  algorithm  was  not  changed.  The  minimum  algorithm  would  not  require 
either  of  the  filters  or  10  samples  to  determine  the  disturbance  frequency  and  would 
therefore  not  be  limited  to  the  33  Hz  operational  effectiveness  limit.  The  “Contiguous 
Copy”  blocks  are  required  for  the  vector/matrix  manipulation. 

d.  Nodal  Accelerations 

Each  of  the  four  accelerometers  used  on  the  selectable  control  nodes  have 
three-axis  capability.  Each  axis  from  the  four  accelerometers  is  fed  into  a  12-channel 
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Kistler  Piezotron  Coupler,  Model  5124A,  and  then  to  a  Trek  50/750  Voltage  Amplifier. 
Details  of  the  equipment  are  included  in  Appendix  A. 


Nodal  Axis 

Switch  Switch 


Figure  19.  Control  Node  and  Axis  Selection 
After  the  proper  signal  conditioning,  the  signals  are  used  as  the  control 
signal  input.  The  values  of  the  constants  “Node  Selector”  and  “Axis  Selector”  are  both 
controlled  fi'om  the  dSPACE  ControlDesk.  The  “Nodal  Switch”  is  controlled  by  the 
integer  value  of  the  “Node  Selector”  firom  1  to  4  that  corresponds  to  the  desired  node 
(Node  26  =  1,  Node  49  =  2,  Node  18  =  3,  and  Node  41  =4).  The  “Axis  Switch”  is 
similarly  controlled  by  the  value  of  “Axis  Selector”  (X-Axis  =  1,  Y-Axis  =  2,  Z-Axis  =  3, 
Max-Axis  =  4)  but  includes  an  automatic  selection  option.  The  output  fi'om  this  sequence 
of  switches  is  one  of  three  axis  signals  originating  fi'om  one  of  four  nodal  accelerometers. 

It  was  found  through  several  test  runs  that  the  input  contained  a  DC  bias 
that  created  a  DC  bias  in  the  output  and  thus  reduced  the  effective  stroke  of  the  active 
control  element.  This  DC  bias  was  removed  with  high-pass  filter  in  the  signal  path  in 
order  to  maximize  the  effectiveness  of  the  controller.  This  need  is  consistent  with  the 
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results  observed  by  the  researchers  at  NRL,  who  used  a  second-order  Butterworth  filter  to 
remove  the  DC  bias  [Ref.  6,  p.  6].  A  fourth-order  filter  was  used  for  the  NPS  space  truss. 
Higher  order  filters  were  tested  to  determine  their  utility  and  it  was  discovered  that  the 
higher-order  did  not  lead  to  an  appreciable  increase  in  filtering,  but  did  decrease  the 
instantaneous  response  time. 


Figure  20.  Automatic  Maximum  Amphtude  Axis  Selection 
In  general,  minimizing  a  single  axis  disturbance  using  a  random  control 


location  will  result  in  an  increase  in  the  magnitude  of  at  least  one  of  the  other  axis.  This 
motivated  the  rather  simple  option  of  having  an  algorithm  automatically  select  the  largest 
magnitude  axis  and  use  that  signal  as  the  controller  error.  Figure  20  shows  the  method 
used  to  determine  the  axis  with  the  maximum  amplitude  disturbance.  The  standard 
deviation  was  used  because  it  is  proportional  to  the  amplitude  of  a  sinusoidal  signal  and 
the  100  sample  buffer  allowed  for  0.1  second  updates  and  covered  enough  of  a  10  Hz 
disturbance  signal  to  be  effective.  The  “max  sfun”  block  outputs  the  index  of  the  signal 
with  the  largest  magnitude  and  is  used  to  control  the  “Multiport  Switch”  that  determines 
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the  axis  sent  to  the  controller.  This  results  in  an  instantaneous  switching  of  the  control 
signal  as  the  amplitude  of  one  of  the  signals  drops  below  another.  The  controller 
tolerates  this  due  to  the  set  adaptation  rates  of  the  LMS  algorithm.  The  use  of  this  option 
predictably  results  in  a  steady  state  disturbance  of  equal  magnitude  in  two  of  the  control 
axis.  This  option  was  not  used  for  evaluation  of  the  controller  due  to  large  the  effect  on 
controller  performance  due  to  dependence  on  specific  system  response. 

e.  Adaptive  Gain  Adjustment 

The  adaptive  gain  adjust  has  the  net  effect  of  controlling  the  convergence 
rate  of  the  entire  controller.  This  element  uses  two  main  adaptive  elements,  but  both  use 
signals  derived  from  the  magnitude  of  the  error  signal  used  by  the  phase  adjustment 
algorithm.  This  turned  out  to  be  one  of  the  more  complex  sections  of  the  controller, 
using  both  simple  logic  statements  and  LMS  algorithms  to  adjust  the  gain  while  the 
controller  was  attempting  to  reduce  the  error.  An  overview  of  this  function  is  shown  in 
Figure  21. 
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Figure  21.  Adaptive  Gain  Adjustment 
The  “Control  Signal  Output”  is  simply  the  product  of  the  “Nodal  Input” 
and  the  calculated  gain.  If  the  controller  is  reducing  the  magnitude  of  the  error  signal, 
then  the  gain  is  increased,  which  increases  the  convergence  rate  of  the  phase  adjust 
algorithm.  Conversely,  if  the  magmtude  of  the  error  is  increasing,  the  gain  is  reduced, 
thus  slowing  the  controller  convergence  rate.  To  derive  a  signal  proportional  to  the 
disturbance  magnitude,  the  standard  deviation  was  calculated  as  shown  in  Figure  22. 


Figure  22.  Amplitude  Calculation 

The  LMS  algorithm  used  to  adjust  the  magmtude  of  the  gain  requires  an 
error  signal  that  oscillates  around  zero  in  order  to  converge.  This  was  accomplished  by 
using  the  product  of  the  error  magnitude  and  it’s  derivative.  This  produced  a  LMS  step 
size  that  was  positive  and  increasing  in  size  when  the  error  was  increasing,  and  negative 
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and  decreasing  in  size  as  the  error  was  decreasing.  To  force  the  gain  up  as  the  error  was 
reduced,  the  sign  of  the  derivative  was  reversed.  Attempts  to  use  the  amplitude  of  the 
derivative  to  increase  the  convergence  rate  resulted  in  unpredictable  and  thus  undesirable 
results.  To  correct  for  this,  the  derivative  signal  was  used  to  determine  the  sign  of  the 
error  signal  only  and  as  noted  before,  the  signs  were  reversed  to  allow  for  LMS 
convergence. 


Normalization  Trend 


Value 

(Maximum) 

Figure  23.  Amplitude  Derivative 

In  order  to  further  stabilize  the  convergence  rate  of  the  controller,  the 
amplitude  of  the  error  magnitude  was  normalized.  The  normalization  value  is  referred  to 
as  the  “Historical  Max  Amplitude”  and  limits  the  magnitude  of  the  error  signal  used  to 
adjust  the  gain  LMS  algorithm.  In  order  to  allow  for  multiple  experimental  runs  without 
reloading  the  program,  a  reset  was  included.  The  activation  of  this  reset  is  either  manual 
from  the  dSPACE  ControlDesk,  or  automatic  based  on  the  stability  of  the  frequency 
identification  signal.  Exact  operation  will  be  explained  later,  but  the  net  effect  on  the 
“amplitude  Process”  block  is  to  trigger  the  “Hold  Max  Value”  block  continuously,  but 
assign  a  value  of  1x10'^®  in  order  to  assure  the  first  calculated  amplitude  will  be  greater 
and  replace  the  reset  value. 
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Figure  24.  Amplitude  Processing 

When  in  normal  operation,  the  historical  maximum  is  determined  by 
taking  the  difference  of  the  current  maximum  and  the  “Amplitude  Input”  signal.  If  the 
input  signal  is  smaller,  then  nothing  happens.  If  the  input  signal  is  larger  than  the 
maximum,  then  it  is  switched  to  the  input  of  the  sample  and  hold  block  and  the  ‘Tulse” 
signal  is  aligned  to  the  trigger.  This  replaces  the  output  of  the  “Hold  Max  Value”  block 
with  the  input  signal.  The  maximum  value  is  used  directly  as  an  output  from  this  block 
and  to  normalize  the  “Amplitude  Input”  signal,  as  shown  in  Figure  24. 

The  adaptive  gain  adjust  also  contains  a  sign  correction  factor  that  is  based 
on  the  historical  maximum  value  compared  to  the  current  magnitude  of  the  error  and  a 
time  delay.  The  input  to  the  algorithm  is  the  “Current  Amplitude”  and  the  “Historical 
Max  Amplitude”  of  the  controller  disturbance  signal.  The  output  is  limited  to  either  V- 1. 
The  controller  is  allowed  1  sec  to  reduce  the  magnitude  of  error  signal  before  the  sign  is 
reversed.  This  time  delay  is  controlled  by  the  pulse  width  of  the  “Trigger”  to  the  “Hold 
the  Sign”  sample  and  hold  block.  The  sign  is  reversed  by  a  simple  application  of  1'", 
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where  n  is  a  counter  signal  that  is  triggered  by  the  same  pulse  generator  that  triggers  the 
sample  and  hold,  as  shown  in  Figure  25. 


Trlger 
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Figure  25.  Rapid  Sign  Correction 

The  “Adaptive  Gain  Fine  Adjustment”  block  controls  the  magnitude  of  the 
gain  as  seen  in  Figure  26.  The  values  internal  to  the  “Gain  Adaptive  Filter”  are  an  initial 
tap  value  (w)  of  400  and  a  step  size  (p’)  of  1x10'^.  These  values  are  reset  by  the  same 
“Reset  Signal”  used  to  restore  the  initial  values  of  the  other  control  systems,  based  on 
wither  manual  override  or  frequency  identification  fluctuation. 


Reset  Speed 


Reset 
T^p  Value 


Figure  26.  Adaptive  Gain  Fine  Adjustment 
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The  “Gain  Stability”  algorithm  shown  in  Figure  27  was  added  to  reduce 
the  fluctuations  of  the  LMS  algorithm  output  after  the  controller  has  effectively 
converged.  It  was  observed  that  the  gain  magnitude  adjustment  continually  floated  at,  or 
near  a  convergence  value.  Unfortunately,  this  oscillation  was  reflected  in  the  phase 
adjust  LMS  and  resulted  in  the  reduced  effectiveness  of  the  controller.  To  prevent  this,  a 
dead  zone  was  built  into  the  gain  magnitude  adjustment.  The  purpose  of  this  block  is  to 
measure  the  fluctuation  in  the  gain  signal  and  produce  a  pulsed  trigger  if  the  average 
difference  between  four  successive  gain  values  is  less  than  4.  When  the  trigger  is  pulsed, 
the  sample  and  hold  block  labeled  “Gain  Lock”  in  Figure  21  will  lock  in  the  current  gain 
value.  The  “Adjust/Lock  Switch”  is  triggered  by  a  gain  stability  value  of  4  or  less  to 
switch  the  gain  used  to  the  output  of  the  “Gain  Lock”  block. 

f.  Adaptive  Phase  Adjustment 

The  phase  adjustment  of  the  control  signal  is  accomplished  using  a  sine 
and  cosine  generated  based  on  the  output  of  the  frequency  identification  block 
Reference  Signal.”  The  same  signal  also  is  used  to  reset  the  controller  presets  if  the 
disturbance  frequency  varies.  The  overview  is  shown  in  Figure  28. 
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Figure  28.  Controller 

The  “Reset  Sum”  block  allows  for  either  a  manual  or  automatic  reset  of 
the  controller.  All  of  the  reset  switches  associated  with  the  LMS  algorithms  and  the 
disturbance  magnitude  history  have  a  set  point  of  1 .  If  the  reset  signal  is  less  than  the  set 
point,  then  the  appropriate  operating  signal  is  allowed  to  pass.  If  the  reset  signal  exceeds 
that  set  point,  then  the  signals  are  switched  to  the  appropriate  values  to  restore  the  initial 
conditions  to  the  controller. 

The  sine  and  cosine  waves  are  produced  as  described  in  section  [III.B.  1]. 
An  additional  gain  factor  was  included  in  this  section  to  allow  for  the  controller  output  to 
be  zeroed.  Switching  both  the  “Sine  Switch”  and  the  “Cosine  Switch”  to  zero  does  this, 
forcing  both  phase  adjustment  LMS  outputs  to  zero  regardless  of  error  signal  or  gain 


Hold  Cos 

Figure  29.  Sine/Cosine  Generator 
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This  controller  requires  steady  reference  signals,  at  the  correct  frequency, 
in  order  for  the  appropriate  gains  to  be  determined.  If  the  reference  signals  fluctuate,  the 
fluctuation  tends  to  be  amplified  by  the  controller.  To  prevent  this,  the  stability  of  the 
frequency  identification  signal  is  monitored.  If  the  signal  is  fluctuation  too  fast,  the 
controller  resets  and  waits  until  the  identification  becomes  stable.  The  “Frequency 
Stability  LMS  Reset”  block  shown  in  Figure  30  is  tuned  to  allow  only  very  minor 
fluctuation  in  the  reference  signal  prior  to  resetting  the  controller.  It  was  found 
experimentally  that  the  controller  was  capable  of  tracking  with  very  slowly  changing 
disturbance  frequencies,  but  this  option  was  not  explored  by  this  research  and  in  fact,  the 
values  of  the  stability  reset  were  set  to  preclude  most  chance  of  this. 


Figure  3  0.  Frequency  Stability  LMS  Reset 
The  reset  is  based  on  the  difference  between  the  mean  of  the  pre\'ious  1 00 
samples  and  the  current  frequency  identified  by  the  algorithm.  If  the  magnitude  of  this 
difference  is  greater  than  0.001  then  the  controller  is  reset.  This  block  generates  a  signal 
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proportional  to  this  difference  with  a  gain  factor  that  brings  the  output  up  to  the  required 
trip  value  for  the  controller  reset  switches.  The  mechanism  for  the  reset  can  be  seen  in 
Figure  28. 

Figure  31  shows  the  “Adaptive  Phase  Adjustment”  block  that  is  essentially 
a  notch  filter  as  explained  in  section  [III.B.  1].  The  addition  of  the  various  reset 
swatches  allows  for  the  reference  and  error  signals  to  be  changed  in  order  to  force  the 
“taps”  back  to  their  default  value  of  zero.  The  step  size  (p’)  used  for  these  LMS  blocks  is 
1x10"^.  This  section  also  contains  a  provision  to  prevent  an  over  voltage  condition  of  the 
piezoelectric  strut.  The  satiuration  limiter  is  configured  to  insure  the  control  signal  in 
maintained  within  the  operation  limits  of  the  strut.  The  DC  bias  is  included  as  part  of  the 
control  signal  in  order  to  allow  the  maximum  oscillation  of  the  strut. 


Reset 


Figure  31.  Adaptive  Phase  Adjustment 
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g.  Active  Element  Selection 


The  last  segment  of  the  controller,  prior  to  the  digital  to  analog  converter 
that  directs  the  signal  to  the  strut,  is  the  “Active  Element  Selection”  block.  This  block 
was  included  to  allow  for  software  switching  between  active  elements.  The  same  DC 
Bias  is  included  in  this  feature  to  minimize  the  voltage  fluctuation  on  the  piezoelectric 
elements. 


Figure  32.  Active  Element  Selection 


2.  RTILIB 


The  interface  from  Simulink  to  dSPACE  is  found  in  the  Real-Time  Integration 
Simulmk  Library  (RTILIB)  and  was  accomplished  the  same  was  as  previous  control 
system  experiments  using  the  NFS  Space  Truss  [Ref  5].  The  program  is  started  by 
typing  “rtilib”  in  the  MATLAB  command  window,  after  Simulink  is  running  [Ref  18,  p. 
17].  The  RTILIB  contains  the  blocks  required  to  provide  the  interface  between  the  CPU 
and  the  I/O  hardware.  The  RTILIB  interface  is  shown  in  Figure  33. 
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Figure  33.  RTILIB  Interface  [Ref.  18,  p.  18] 


The  blocks  were  selected  and  then  “dragged  and  dropped”  onto  the  Simulink 
model  for  integration.  The  VO  box  connections  are  found  under  the  “MASTER  PPC” 
icon  and  there  are  twenty  available  inputs  and  eight  outputs  available  to  the  user  in 
dSPACE,  accessible  by  three  blocks,  shown  in  Figure  34. 


Figure  34.  MASTER  PPC  FO  Connections  [Ref  18] 


The  ADC  channels  available  are  the  ADC  and  MUX  (Multiplexed)  ADC.  The 
MUX  ADC  is  a  series  of  four  channels,  with  four  multiplexed  inputs  each  (ADCH  1-16). 
They  can  be  used  in  Simulink  as  individual  separated  charmels,  or  as  a  vectored  input. 
The  scalar  ADC  block  is  a  single  input,  corresponding  to  the  last  four  input  channels 
(ACDH  17-20).  There  are  eight  output  channels  accessible  with  the  DAC  blocks  (DACH 


1-8)  [Ref.  18,  p.  24].  The  FO  blocks,  when  connected  to  the  FO  interface  box,  exhibit  a 
10:1  gain  in  value. 

D.  dSPACE  CONTROL  DESK 

ControlDesk  also  allows  for  grouping  relevant  files  under  an  “experiment”  and 
then  accessing  and  manipulating  the  data  loaded  in  the  CPU.  It  also  allows  for  running 
macros  that  automate  the  data  taking  process  using  the  dSPACE  macro  language.  Python 
[Ref  19,  pp.  13-14]. 

Using  ControlDesk  allows  for  the  selection  of  variables  for  real-time  display  in 
on-screen  instruments,  and  for  entire  program  to  be  downloaded  and  executed  by  the 
CPU  [Ref  19,  p.  3 1].  To  create  a  working  experiment,  a  program  must  be  loaded  into  the 
CPU’s  memory,  then  an  instrument  panel  may  be  built  that  will  monitor  and  control  the 
parameters  and  signals  desired.  The  ControlDesk  interface  contains  the  main  window. 
Navigator,  and  Tool  Window  that  can  be  manipulated  to  allow  a  user  to  run  his/her 
experiments.  The  Navigator  provides  access  to  load  files  into  the  CPU;  the  Tool  window 
is  used  to  create  the  interface  using  the  program  variables  and  several  interface  options. 
The  Main  window  is  used  to  monitor  and  control  the  simulation  [Ref  19,  pp.43-50]. 

Instrument  panels  are  built  in  the  Tool  window  using  a  pull-down  menu  interface 
to  D&D  desired  instrument  on  the  main  window.  The  user  can  choose  between  sliders, 
pointers,  oscilloscopes,  knobs,  buttons... etc.  to  create  the  interface.  There  are  various 
t3^es  of  variables  that  can  be  assigned  to  each  of  these  interface  options.  The  variables 
are  grouped  according  to  their  block  location  in  the  Simulink  model  used  to  create  the 
program  loaded  into  the  CPU.  The  relevant  types  to  this  thesis  are  shown  in  Table  2. 
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Type _ Description: _ 

B: _ Block  Outputs 

S: _ Inputs  of  Signal  Sinks 

P:  I  Block  Parameters 

Table  2.  dSPACE  Variable  Types  [Ref.  20] 

The  controller  layout  used  for  the  Adaptive  Multi-layer  LMS  controller  on  the 
NPS  Space  Truss  is  shown  in  Figure  35.  The  large  plot  on  the  top  of  the  Layout  is  an 
oscilloscope  plot  of  the  disturbance  signal.  The  upper  plot  in  the  lower  right-hand 
portion  of  the  layout  is  the  control  signal  to  the  actuator.  The  lower  plot  is  to  monitor  the 
convergence  of  the  LMS  algorithms  for  both  phase  and  gain  adjustment.  The  actuator 
used  for  control  is  designated  in  the  selection  labeled  “Active  Control  Element”  located 
in  the  center  of  the  panel.  The  “Gain  Lock”  LED  is  controlled  by  the  stability  signal 
generated  by  the  “Gain  Stability”  block  (section  III.C.  l.e)  and  turns  from  yellow  to 
green  when  the  gain  stops  adjusting.  The  control  node  and  axis  are  selected  with  buttons 
in  the  “Node  and  Axis  of  Interest”  section  located  just  left  of  center.  The  “Auto  Axis” 
LED  is  controlled  by  the  “Automatic  Maximum  Amplitude  Axis  Selection”  block 
(section  III.C.  l.d)  and  changes  color  depending  on  the  relative  magnitude  of  the 
disturbance  and  indicates  RED  when  the  X-axis  is  the  largest,  WHITE  when  the  Y-Axis 
is  the  largest,  and  BLUE  when  the  Z-Axis  is  the  largest.  The  output  from  the  ‘Trequency 
Identification”  block  (section  III.C.  l.c)  is  displayed  in  by  the  digital  readout  labeled 
“Frequency  Identification”  located  on  the  center,  left-hand  edge  of  the  panel.  The 
“Run/Reset”  LED  is  controlled  by  the  ‘Reset  Sum”  value  (section  III.C.  l.f)  and  turns 
green  when  the  controller  is  active,  and  yellow  when  it  is  in  standby  either  due  to 
frequency  identification  instability,  or  manual  override. 
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The  disturbance  frequency  and  amplitude  can  be  adjusted  from  the  “Disturbance 
Source  Control”  block  located  in  the  lower  left  comer.  The  only  adjustment  required  for 
the  controller  to  function  is  to  push  the  “Run”  button  on  the  “Reset”  switch  located  to  the 
left  of  the  “Run/Reset”  LED.  All  other  controls  change  the  disturbance  properties,  the 
control  element  location,  or  the  location  of  the  node  and/or  the  axis  where  the  vibration  is 
to  be  minimized.  The  steps  required  to  start  and  run  the  controller  are  included  in 
Appendix  D. 

E.  AMPLIFIERS  AND  CONDITIONERS 

As  previously  stated,  the  equipment  used  for  this  research  is  virtually  identical  to 
that  used  by  Pantling  [Ref  5].  A  notable  exception  to  this  is  the  addition  of  another 
active  stmt  and  the  associated  signal  conditioners.  The  signals  generated  by  the  foree 


transducers  are  not  used  for  this  controller,  but  are  still  available  if  required  for  future 
research.  These  signals  are  fed  into  PCB  Piezotronics  Model  484B  signal  conditioner, 
which  has  a  unity  gain  and  a  selectable  bias  (either  6.0  or  11.0  Volts,  DC  or  AC).  The 
four  accelerometers  mounted  on  the  control  nodes  have  three-axis  capability.  Each  axis 
from  the  four  accelerometers  is  fed  into  a  12-channel  Kistler  Piezotron  Coupler,  Model 
5124A.  The  control  signal  is  sent  to  a  Trek  50/750  Voltage  Amplifier.  Details  of  the 
equipment  are  included  in  Appendix  A. 

The  mini-calibration  performed  on  the  Trek  amplifier  by  Pantling  [Ref.  5]  was 
determined  to  be  sufficient  for  this  research.  The  Trek  was  calibrated  for  a  20-time 
voltage  gain. 
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IV.  CONTROLLED  FINITE  ELEMENT  MODEL  OF  THE  NFS  SPACE 

TRUSS 


A.  ANSYS  OVERVIEW 

ANSYS  is  a  finite  element  utility  capable  of  modeling  many  different  physical 
properties,  using  more  than  100  types  of  elements.  Systems  that  involve,  mechanical, 
electromagnetic,  thermal,  piezoelectric  and  electric  characteristics  can  be  modeled.  The 
options  available  with  the  software  package  are  based  on  the  licensing  obtained  from 
ANSYS.  For  the  work  in  this  thesis  ANSYS  version  5.6,  Multiphysics  was  used. 

The  ANSYS  preprocessor  permits  the  creation  and  manipulation  of  a  finite 
element  model  that  can  then  be  meshed  into  a  finite  element  grid  and  acted  upon  with 
loads  and  restraints.  The  solver  sub-program  can  be  used  to  specify  different  solution 
types.  For  this  thesis,  a  modal  analysis  was  performed  on  the  same  model  used  for  the 
transient  analysis  on  the  controlled  model.  The  specifics  of  the  ANSYS  commands 
required  are  detailed  in  Reference  21. 

B.  ANSYS  PARAMETRIC  DESIGN  LANGUAGE  (APDL) 

The  ANSYS  Parametric  Design  Language  (APDL)  is  a  macro  language  that 
enables  a  user  to  run  command  streams  from  an  external  file.  This  allows  for  a  closed 
model  that  performs  an  actively  controlled  data  run,  requiring  no  input  from  the  user  save 
to  establish  initial  parameters  and  start  the  simulation. 

An  APDL  file  is  a  text  file  that  contains  the  commands,  as  they  would  be  t>pcd 
into  the  ANSYS  command  interface,  with  user-determined  variables  [Ref  21,  /INPUT, 

*SET]  and  is  used  in  conjunction  with  the  transient  analysis  capability  to  perform  the 
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control  law  application.  ANSYS  allows  the  user  to  examine  each  step  in  a  transient 
analysis  and  the  APDL  program  can  record  these  parameters  to  output  files  for  analysis  in 
the  time  domain. 

The  model  routine  automatically  time-steps  through  a  transient  analysis,  applying 
a  sinusoidal  force  to  the  mass  representing  the  LPACT.  The  control  node  displacement 
was  also  measured  and  the  velocity  and  acceleration  of  that  node  estimated  based  on 
previous  samples.  The  nodal  acceleration  was  then  used  as  the  error  signal  for  the 
controller  section  of  the  adaptive  multi-layer  LMS  controller.  The  entire  logic  statement 
of  the  control  algorithm  was  not  included  for  this  thesis  for  simplicity.  The  numerical 
control  law  was  applied  and  the  control  voltage  was  calculated  and  applied  to  the  chosen 
piezoelectric  actuator.  These  varying  forces  were  recalculated  and  applied  at  each  step. 

1.  Controller  Application 

An  example  of  the  APDL  files  used  in  this  thesis  is  included  as  Appendix  E  and  is 
based  on  the  finial  APDL  used  by  Pantling  [Ref  5].  There  was  no  need  for  the  digital 
filter  used  in  Reference  5,  so  it  was  removed.  The  adaptive  multi-layer  LMS  control 
algorithm,  as  it  is  applied  to  the  NPS  Space  Truss,  uses  nodal  acceleration  to  generate  a 
disturbance  reference  signal.  This  is  not  a  direct  output  of  the  ANSYS  program,  but  the 
nodal  displacement  is.  After  each  time  step,  the  post-processor  was  used  to  obtain  the 
vector  displacement  at  the  control  node.  Then  the  velocity  is  estimated  using  a  backward 
difference,  2"^  order  derivative  approximation. 
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(4.1) 


f 

"  At 
where, 

=  The  vector  change  in  displacement. 
At  =  The  time  step  used  in  the  analysis. 


The  acceleration  is  calculated  in  the  same  manor,  replacing  the  estimated  velocity 
for  the  measured  displacement  and  using  the  same  time  step. 


(K-K-.) 

At 


where, 

(V„  -  F„_j)  =  The  vector  change  in  velocity. 


(4.2) 


For  the  sake  of  simplicity,  the  gain  portion  of  the  adaptive  multi-layer  LMS 
control  algorithm  was  not  included  in  the  simulation.  This  did  require  multiple  runs  to 
establish  the  optimal  gain  that  should  be  applied  to  the  controller,  but  the  saved 
programming  time.  Further  research  in  this  area  could  incorporate  the  entire  control 
system  into  the  APDL  file. 

The  controller  uses  the  LMS  algorithms  developed  in  Reference  5  and  outlined  in 
section  [III.B.  1],  Adaptive  Phase  Adjustment.  The  basic  equations  are  restated  below: 

CS.=G[a.  cos(2^fnAt)  +  sm{27rfnAt)^ 
where, 

CS„  =  Voltage  applied  to  the  piezoelectric  element. 

G  =  Constant  gain  factor  due  to  hardware  amplification. 

/=  Disturbance  frequency  in  Hz.  (4.3) 

n  =  Iteration  counter. 

At  =  The  time  step  used  in  seconds. 

<2„  s  The  Cosine  weighting  factor. 

=  The  Sine  weighting  factor. 
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The  weighting  factors  are  updated  each  time  step  based  on  the  calculated  nodal 
acceleration  (An)  using  the  following  equations; 

+  M  cos  (4.4) 

and 

sin  (4.5) 

The  adaptation  constant  (p)  was  determined  by  trial  and  error,  running  short 
simulations  and  evaluating  the  results  until  the  initial  trends  were  indicative  of 
convergence.  The  finial  run  was  then  made  using  the  predetermined  value  of  p. 

After  each  time  step  is  solved  for  the  CS  voltage  was  applied  to  the  actuator 
element  in  the  next  time  step.  The  control  signal,  weighting  factors,  nodal  displacement, 
velocity,  and  acceleration  were  all  captured  for  analysis. 

C.  NPS  SPACE  TRUSS  MODEL 

The  FEM  created  in  Reference  5  was  used  for  this  analysis  with  an  additional 
active  strut  assembly  replacing  NPS  Space  Truss  element  107.  First  the  existing  truss 
element  was  deleted.  Then  the  piezoelectric  element  and  force  sensor  were  copied  and 
reflected  across  a  plane  of  symmetry  located  in  the  middle  of  the  truss.  Then  the 
elements  were  rotated  90'  around  the  geometric  center  of  the  span.  Then  a  similar 
procedure  was  followed  for  the  steel  connecting  elements.  This  effectively  made  the 
model  look  correct,  but  the  mathematical  formulation  was  incomplete. 

The  next  step  was  to  merge  the  coincident  nodes  created  by  copying  and  pasting 
the  original  active  element  over  the  previously  existing  nodes  on  the  opposite  side.  This 
procedure  also  had  to  include  the  nodal  junctions  of  the  piezoceramic  elements  and  the 
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steel  support  elements.  After  the  nodes  were  merged,  the  boundary  conditions  had  to  be 
applied  to  the  new  piezoelectric  elements.  Boundary  conditions  do  not  copy  with  the 
elements.  As  with  the  previous  model,  a  grounded  condition  was  used  to  act  as  a 
reference  to  the  applied  voltage  and  then  the  displacement  DOF’s  of  the  solid  elements 
were  coupled  to  the  steel  beam  element  nodes  that  contacted  the  piezoelectric  model. 

Details  of  the  original  model  assembly  can  be  found  in  Reference  5  and  the  actual 
component  specifications  are  found  in  Appendix  A.  A  summary  of  the  properties  used  to 
model  the  piezoelectric  properties  is  found  in  Table  3. 


Property 

Quartz 

PZT 

Stiffiiess 

l.OkN/pm 

Sensitivity 

1 1,420  mV/kN 

- 

Axial  Elastic  Modulus 

lOOGpa 

142GPa 

d33 

87.56e-9  mA^ 

500e-9  m/V 

633 

8756 

7100 

Table  3.  Material  Properties  [From  Ref  5] 


Figure  36  shows  one  of  the  active  struts  as  modeled  in  the  ANSYS  FEM. 


Figure  36.  Active  Strut  Modeled  in  the  NPS  Space  Truss  [From  Ref  5] 
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The  LPACT  model  was  not  modified,  but  the  experimental  results  did  show  that 
the  limits  of  the  control  authority  of  the  installed  elements  could  be  exceeded  by  the 
LPACT  at  certain  fi'equency  and  node  combinations.  This  was  addressed  by  reducing  the 
amplitude  of  the  reference  signal  to  the  LPACT.  The  same  transfer  fimction  derived  in 
Reference  5  was  used  to  model  the  force  applied  by  the  LPACT  in  this  research,  however 
it  was  occasionally  required  to  modify  the  amplitude  by  applying  a  gain  factor  (GF)  in 
the  FEM.  This  resulted  in  the  following  equation  that  was  used  in  the  APDL  program: 

AMP  =  (GF)  (0.2037  *  FREQ^  -  7.071 9  *  FREQ  +  68.3564)  [Ib/amp]  (4.6) 

The  APDL  macro  program,  “Act_truss.inp”  that  was  used  in  the  previous  analysis 
[Ref.  5]  was  modified  to  become  Tiuss_Control.txt”,  which  served  as  the  root  program 
for  this  analysis.  The  original  program  “Act_trass.inp”  is  included  for  reference  in 
Appendix  F  and  the  modified  version  (set  up  to  control  the  y-axis  of  Node  18,  with  19Hz 
disturbance,  using  the  second  active  element)  is  included  in  Appendix  E. 

The  new  version  of  ANSYS  seemed  to  eliminate  the  stability  problems  that  were 
experienced  in  previous  research.  No  stability  problems  were  encoxmtered  in  any  of  the 
simulations  that  were  not  explained  by  incorrect  gains  being  applied  to  the  control 
algorithm. 
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V.  COMPARISON  OF  EXPERIMENTAL  RESULTS  WITH  ANSYS  FINITE 

ELEMENT  MODEL 


A.  METHODOLOGY 

The  purpose  of  these  experiments  is  to  demonstrate  the  performance  of  the 
controller  developed  in  this  thesis  and  to  confirm  the  ability  of  an  active  FEM  to  predict 
the  effectiveness  of  the  algorithm.  The  flexibility  of  the  controller  made  a  full 
demonstration  of  its  abilities  infeasible.  In  order  to  show  its  general  abilities,  nine 
experiments  were  performed  using  a  variety  of  frequency,  control  node,  and  active 
element  combinations.  Eight  of  the  experimental  combinations  were  also  simulated  in 
ANSYS  for  comparison.  Although  the  adaptive  gain  portion  of  the  controller  is  not 
included  in  the  ANSYS  model,  the  correct  gains  for  the  comparisons  were  determined  on 
a  trial  and  error  basis  using  the  results  of  short  simulations  to  modify  and  time  the 
controller  gain  manually.  A  ninth  experiment  was  run  to  compare  the  results  of  the 
adaptive  multi-layer  LMS  controller  with  the  previously  used  IDIFF  controller. 

The  first  four  experiments  were  designed  to  demonstrate  the  versatihty  of  the 
controller  with  respect  to  active  element  location,  disturbance  frequency,  and  control 
location.  The  first  and  second  experiments  control  different  nodes,  axis,  and  fi-equencies 
using  different  active  elements.  The  third  and  fourth  experiments  control  different  axis  of 
the  node  that  is  the  closest  to  the  disturbance  source  with  the  same  disturbance  fi-equency, 
using  different  active  elements.  The  next  four  experiments  were  done  to  demonstrate 
more  challenging  control  situations.  They  use  the  first  three  modal  fi'equencies  predicted 
by  the  FEM,  axis  based  on  the  displacement  shown  by  modal  analysis,  and  a  randomly 
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chosen  active  element.  The  last  experiment  was  done  at  the  frequency,  node  and  axis  at 
which  the  IDIFF  controller  demonstrated  the  best  performance. 

Each  experiment  was  run  multiple  times  with  nearly  identical  results  in  most 
cases.  A  representative  set  of  results  of  each  data  is  included  either  in  the  text,  or  in 
Appendix  G. 

B.  EXPERIMENTAL  CONTROL  SYSTEM 

1.  Experimental  Setup 

The  nodal  accelerations  used  as  the  controller  system  feedback  were  generated  by 
one  of  four  Kistler  8690C10  three-axis  accelerometers.  Two  accelerometers  were  place 
at  the  ends  of  the  truss,  nodes  26  and  41.  The  other  two  were  positioned  near  the 
midpoint  of  the  truss  arms,  nodes  49  and  18.  This  configuration  is  the  same  as  the  one 
used  in  previous  research,  but  the  node  numbering  is  not  consistent.  If  direct  comparison 

of  results  is  desired,  both  the  nodal  numbering  system  and  axis  orientation  should  be 
verified. 

A  Kistler  Piezotron  Signal  Conditioner,  model  5124A  was  used  as  a  signal 
conditioner  for  the  accelerometers.  It  does  not  amplify  the  signal  and  only  served  to 
power  the  accelerometers  and  filter  the  output.  Accelerometer  drift  that  was  experienced 
in  previous  research  [Ref.  5]  was  removed  using  a  high  pass  filter.  This  was  not  required 
in  the  past,  because  the  signals  were  used  only  for  spectral  analysis.  This  controller  uses 
the  signals  as  feedback  and  would  attempt  to  zero  out  the  low  frequency  drift  as  well  as 
the  vibration  and  this  could  have  lead  to  actuator  saturation. 
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The  dSPACE  DS1103  digital  signal  processor  received  ADC  inputs  from  the 
accelerometer  signal  conditioners.  The  control  signal  DAC  output  was  feed  to  a  TREK 
power  supply,  model  50/750  and  then  directly  to  the  piezoelectric  elements.  A  second 
DAC  ou^ut  was  used  to  send  the  nodal  acceleration  signal  to  a  Hewlett  Packard  HP 
35665A  Dynamic  Signal  Analyzer.  This  output  was  not  required  for  the  operation  of  the 
controller  and  only  served  to  provide  a  visual  signal  representative  of  the  performance  of 
the  controller  that  did  not  use  the  same  computational  resources  as  the  control  system. 

The  generation  of  the  excitation  signal  for  the  LPACT  was  incorporated  into  the 
controller  in  order  to  facilitate  a  common  dSPACE  ControlDesk  capable  of  controlling  all 
required  experimental  parameters.  This  signal  used  a  third  DAC  output  from  the  DSl  103 
board. 

Generic,  off  the  shelf  BNC  connector  cables  of  various  lengths  were  used  to  make 
all  connections.  The  experimental  control  system  can  be  seen  in  Figure  37  and  the 
description  of  coimections  is  provided  in  Table  4. 
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Figure  37.  Experimental  Layout 
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Table  ^ 

k  Cable  Coimections 

The  dSPACE  ControlDesk  used  to  run  the  experiments  is  called  “Adaptive 
Conrtoller.cdx”  and  is  located  in  the  “C:\Space_Truss\TrussOr’  directory  of  the  host 
computer.  This  ControlDesk  was  used  to  run  the  majority  of  the  experiments  and  has  an 
additional  layout  that  allows  the  researcher  to  monitor  all  the  nodal  accelerations  at  once. 
A  similar  ControlDesk  called  “Data  Capture.cdx”  is  located  in  the  same  directory  and 
was  used  to  save  the  data  analyzed  in  this  chapter.  A  separate  ControlDesk  was  used  in 
order  to  reduce  the  computational  requirements  of  displaying  all  the  nodal  accelerations 
at  one  time.  The  Data  Capture  version  does  not  include  the  monitoring  layout.  Both 
versions  have  all  the  controls  required  to  vary  the  frequency,  amplitude,  control  node,  and 
active  element  used  for  each  experiment  (Figure  35). 

After  all  the  required  connections  are  made  and  the  equipment  is  powered  up,  the 
a  Newport  Vibration  Control  System  Table  should  be  “floated”  on  the  compressed 
Nitrogen  supports  to  lower  the  noise  level  prior  to  data  collection.  The  ControlDesk  is 
then  loaded  and  the  switch  LPACT  switch  on  the  LPACT  Driver  is  enabled.  The 
remaining  parameters  can  be  manipulated  via  the  ControlDesk  with  general  security  to 
equipment  integrity,  as  the  required  signal  limits  have  been  built  into  the  controls  as 
hardware  limits  and  cannot  be  exceeded  without  modifying  the  ControlDesk  Layout. 
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After  the  required  disturbance  parameters  are  set,  the  control  node  and  axis  are 
chosen  and  the  active  element  selected.  The  “Run/Reset”  button  activates  the  controller 
indicated  by  the  “Run/Reset”  button  turning  GREEN.  Controller  operation  is  monitored 
by  the  gain  level  values,  the  “Gain  Lock”  LED  and  the  “Control”  and  “Disturbance” 
signal  displays.  Variations  of  the  disturbance  frequency  or  magnitude  will  automatically 
reset  the  controller  and  will  be  obvious  by  the  “Run/Reset”  button  returning  to 
YELLOW.  After  the  frequency  identification  stabilizes,  the  controller  will  automatically 
return  to  “Run”  unless  manually  overridden. 

2.  Modal  Verification 

The  system  modal  frequencies  have  been  verified  with  a  Hewlett-Packard 
HP35665A  Dynamic  Signal  Analyzer  (DSA)  [Ref.  5].  The  results  are  shown  below  and 
correspond  well  with  those  predicted  by  the  ANSYS  FEM  modal  analysis. 


Figure  38.  Frequency  Response  from  Random  Noise  [From  Ref.  5] 
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C.  EXPERIMENTAL  PARAMETERS  AND  RESULTS 


1.  Data  Capture 

Data  capture  was  done  using  the  Data  Acquisition  feature  of  dSPACE 
ControIDesk  [Ref  23],  A  layout  was  added  to  the  experiment  to  add  the  “Capture 
Settings”  window  shown  in  Figure  39. 


Figure  39.  dSPACE  ControIDesk  Data  Capture  Window 


The  PPC”  block  defines  the  set  of  variables  that  you  wish  to  capture  and  the 
length  defines  the  time  in  seconds  to  record  the  data.  The  data  is  selected  for  capture  by 
defimng  the  “Capture  Settings”  for  the  desktop  display  of  that  variable  and  then  selecting 
that  defimtion  in  the  before  mentioned  ‘TPC”  block.  The  data  sample  rate  was  down 
sampled  for  capture  to  reduce  file  size.  The  500  Hz  sample  rate  (1  kHz  down  sampled  by 
2)  was  more  than  sufficient  to  provide  adequate  resolution  below  30  Hz. 

The  “Settings”  button  provides  access  to  the  acquisition  options  available.  A 
“Simple  Capture”  with  “Automatic  Storage”  was  used  for  all  data  capture.  These  options 
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generate  a  “*.mat”  file  and  allow  the  file  name  to  be  preset.  The  automatic  storage 
facilitates  setting  up  the  experiment  fi:om  start  to  finish  prior  to  execution.  Data  was 
written  directly  to  the  working  directory  to  be  used  for  storage  and  manipulation  by 
MATLAB. 

The  sample  times  used  for  each  experiment  were  set  to  allow  a  continues  data  set 
to  hold  the  steady  state,  transition,  and  equilibrium  conditions  of  the  controller  at  each  set 
of  experimental  conditions.  The  limiting  time  was  based  on  having  enough  samples 
before  control  and  after  it  converged  to  allow  for  a  meaningfiil  FFT  to  be  performed.  The 
variable  time  was  the  time  required  for  convergence.  The  data  sets  were  typically  180 
seconds,  but  they  were  extended  up  to  360  seconds  if  required  to  show  system  response 
to  the  controller.  The  parameters  stored  were  the  nodal  acceleration,  the  controller  gain, 
the  sine  LMS  gain,  the  cosine  LMS  gain,  and  the  control  signal  sent  to  the  piezoelectric 
actuator.  The  program  “PROCESS l.M,”  included  as  Appendix  H,  contains  the  code  that 
was  used  to  generate  the  plots  for  the  first  experiment.  The  same  code  was  modified  to 
use  the  sample  length  and  file  names  to  analyze  the  other  experiments. 

The  required  adjustments  to  the  disturbance,  control  location,  and  control  element 
were  all  done  from  the  ControlDesk.  All  experiments  used  the  same  control  logic  and  no 
adjustments  were  made  to  any  gain  or  other  controller  parameter.  The  frequency 
response  shown  for  each  experiment  represent  the  power  spectrum  of  the  nodal 
acceleration  in  g’s  that  the  controller  was  attempting  to  reduce. 
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2.  Experiment  1  (19  Hz,  Y-axis  of  Node  18,  Piezo  #2) 

The  first  experiment  performed  as  a  baseline  for  the  controller  performance  and 
not  specifically  to  test  stress  its  capabilities.  The  controller  was  set  to  “Reset”  and  the 
disturbance  was  set  to  its  maximum  amplitude  at  19.0  Hz.  The  system  was  monitored 
with  the  HP  Dynamic  Signal  Analyzer  until  the  power  spectrum  had  reach  equilibrium. 
Then  the  data  capture  was  started  and  allowed  to  run  for  approximately  15-20  seconds  to 
allow  for  enough  data  to  produce  the  uncontrolled  power  spectrum. 

When  the  controller  was  placed  in  run,  the  adaptive  values  began  to  vary  over 
time.  This  response  was  also  recorded  and  is  shown  in  Figure  40.  After  equilibrium  was 
reached,  the  data  captured  continued  from  long  enough  to  demonstrate  stability  and 
record  the  controlled  system  response.  The  time  history  of  the  disturbance  signal  and  the 
control  signal  are  shown  in  Figure  43  and  the  power  spectrum  of  the  frequency  response 
of  the  acceleration  in  g’s  is  shown  in  Figure  42. 


Time  (sec) 

Figure  40.  Exp.  1,  Controller  Gain  Response 
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47.38  dB  Reduction  at  19  Hz 


This  example  showed  a  47.38  dB  reduction  at  the  disturbance  frequency.  The 
limiting  factor  on  the  magnitude  of  the  reduction  was  later  demonstrated  to  be  initial 
amplitude  and  control  authority.  The  adaptive  nature  of  the  controller,  if  convergence  is 
possible,  will  always  take  the  control  signal  down  to  the  noise  level.  This  example  also 
shows  a  marked  drop  in  the  first  harmonic  frequency,  but  this  is  not  a  design  feature  of 
the  controller  and  was  not  expected. 

3.  Experiment  2  (15  Hz,  Z-axis  of  Node  41,  Piezo  #1) 

The  second  set  of  experimental  conditions  changes  the  all  the  base  parameters 
including,  active  element,  disturbance  frequency,  control  node,  and  control  axis.  This 
was  chosen  to  demonstrate  the  ability  of  the  same  control  algorithm  to  be  used  at  a 
different  frequency,  from  a  different  location,  and  controlling  a  different  parameter  with 
no  reduction  in  performance.  The  controller  initially  converges  slightly  faster  with  these 
conditions,  but  the  finial  power  spectrum  still  shows  a  slight  peak  at  15  Hz.  Even  with 
the  remaining  peak,  the  controller  reduced  the  nodal  axial  acceleration  by  42.45  dB.  The 

results  of  this  experiment  are  shown  in  Figure  60,  Figure  61,  and  Figure  62  of  Appendix 
G. 

4.  Experiment  3  (12  Hz,  X-axis  of  Node  26,  Piezo  #2) 

These  conditions  were  chosen  because  the  LPACT  is  at  the  same  end  of  the  truss 
as  Node  26  and  the  disturbance  source  is  applied  in  a  plane  perpendicular  to  the  x-axis. 
This  makes  the  only  sigmficant  contribution  to  acceleration  in  the  x-direction  due  to 
rotation  about  the  base  due  to  the  induced  vibration.  This  signal  is  significantly  smaller 
than  the  other  two  axis,  which  decreased  the  signal  to  noise  ratio  significantly  as  shown 
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in  Figure  64  of  Appendix  G.  This  combination  of  control  node  and  axis  is  also  one  of  the 
worst  cases  for  control  authority  due  to  the  location  and  orientation  of  the  active  element. 
Regardless,  the  controller  still  reduced  the  signal  by  23.07  dB  and  further  reduction  is 
expected  if  longer  time  is  allowed  for  convergence.  The  system  gain  values  were  not 
fully  stabilized  even  after  about  2.5  minutes  of  adaptation.  System  gain  values  and  power 
spectrum  are  included  in  Figure  63  and  Figure  65  of  Appendix  G  respectively. 

5.  Experiment  4  (12  Hz,  Y-axis  of  Node  26,  Piezo  #1) 

This  set  of  conditions  was  chosen  based  on  the  results  of  the  previous  experiment. 
The  axis  of  concern  was  shifted  to  the  Y-axis  in  order  to  provide  more  control  authority 
to  the  active  element  and  reduce  the  signal  to  noise  ratio.  The  shift  firom  active  element 
number  two,  to  active  element  number  one  was  only  to  verify  that  either  position  could 
control  that  nodal  location.  All  other  parameters  were  left  identical.  The  signal  strength 
was  increased  and  but  the  vibration  reduction  was  only  slightly  improved,  showing  a 
reduction  of  25.24  dB.  The  system  response  is  notably  different  for  the  two  cases. 
Comparing  Figure  63  and  Figure  66  of  Appendix  G,  shows  the  difference  that 
disturbance  magnitudes  make  on  the  convergence  rates  of  the  various  adaptive 
components.  In  experiment  3,  die  signal  was  week  and  the  gain  values  for  sine  and 
cosine  were  still  approaching  their  equilibrium  values  after  3  minutes.  In  this 
experiment,  the  gains  initially  overshoot  the  optimal  value  and  approach  them  in  an 
oscillatory  manor.  A  further  reduction  in  both  cases  may  have  been  realized  if  longer 
convergence  time  was  allotted.  The  system  response  and  power  spectrum  are  shown  in 
Figure  67  and  Figure  68  of  Appendix  G. 
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6.  Experiment  5  (1**  Modal  Freq,  Y-axis  of  Node  26,  Piezo  #1) 

The  original  intent  of  the  controller  was  to  adapt  and  control  vibrations  from  10  to 
30  Hz,  based  on  the  assumption  that  the  first  modal  frequency  was  introduced  by  the 
disturbance  source  only.  Upon  review  of  the  mode  shape  predicted  by  the  revised  FEM 
model  and  the  model  frequencies  of  both  the  predicted  and  experimental  results,  it  was 
decided  that  the  first  structural  mode  was  9.4  Hz.  In  fact,  this  frequency  is  exceptionally 
challenging  due  to  the  nonlinearities  introduced  by  the  disturbance  source.  The  LPACT 
frequency  response  curve  [Figure  10]  shows  its  natural  frequency  at  approximately  8.5 
Hz  with  a  significant  increase  in  response  below  10  Hz. 

For  all  of  the  modal  control  experiments,  node  26  was  used  as  the  control  point. 
This  node  is  the  closest  to  the  disturbance  source  and  thus  presented  the  greatest 
challenge  to  the  controller.  A  disturbance  of  9.4265  Hz  was  used  for  this  experiment 
based  on  the  first  ANSYS  modal  frequency  and  the  Y-axis  (horizontal)  because  it  had  the 
largest  signal  to  noise  ratio  prior  to  control.  The  active  element  used  was  still  chosen 
randomly  and  for  this  experiment  Piezo  #1  was  used. 

The  convergence  rate  of  the  controller  was  significantly  slower  for  this 
experiment,  possibly  due  to  the  nonlinearities  and  the  fact  that  the  first  harmonic  (~19 
Hz)  actually  had  about  20  dB  more  power  than  the  disturbance  frequency.  The  controller 
did  manage  to  reduce  the  disturbance  frequency  by  18.80  dB  with  no  noticeable  affect  on 
the  harmonics  [Figure  71  of  Appendix  G]  even  though  the  one  of  the  base  assumptions  of 
the  control  algorithm  was  a  linear  system.  The  controller  gain  and  system  response  plots 
are  Figure  69  and  Figure  70  of  Appendix  G,  respectively.  The  power  spectrum  for  this 
experiment  is  Figure  71  of  Appendix  G. 
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7.  Experiment  6  (2"*'  Modal  Freq,  Y-axis  of  Node  26,  Piezo  #2) 

For  consistency,  the  only  thing  changed  in  this  experimental  condition  was  the 
disturbance  frequency.  The  second  mode  was  predicted  by  ANSYS  to  be  1 1.318  Hz  and 
that  was  used  as  the  disturbance  source.  Both  active  elements  were  used  with  the  same 
results  and  the  results  included  are  those  using  Piezo  #2. 

These  conditions  also  demonstrated  strong  nonlinearities.  It  is  worth  noting  that 
the  nonlinear  nature  of  the  structural  response  was  not  typical  and  these  frequencies 
would  not  have  been  identified  without  the  ANSYS  analysis. 

The  controller  was  able  to  obtain  a  41.54  dB  reduction  with  a  relatively  slow 
convergence  rate.  The  harmonics  in  this  experiment  did  exhibit  some  net  gain,  but 
remained  less  than  the  original  amplitude  of  the  disturbance  frequency.  Again,  the 
harmonic  response  is  not  included  in  the  controller  design  and  a  linear  system  was 
assumed  in  its  development.  Experimental  results  are  shown  in  Figure  72,  Figiue  73,  and 
Figure  74  of  Appendix  G. 

8.  Experiment  7  (2"**  Modal  Freq,  Z-axis  of  Node  26,  Piezo  #2) 

Based  on  the  bare  truss  modal  analysis,  the  second  modal  frequency  should  be  in 
the  x-z  plane  and  have  no  y  component.  This  observation  is  obviously  not  the  case  due  to 
the  applied  force  in  the  experiments  being  in  the  y-z  plane,  but  it  did  motivate  the 
conditions  set  in  this  experiment.  The  conditions  are  the  same  as  those  in  experiment  6, 
except  that  the  Z-axis  was  selected  for  control. 

The  controller  took  a  relatively  long  time  to  converge  in  this  condition  and 
showed  only  a  2.4  dB  reduction  over  the  3-niinute  time  frame  used  in  all  the  previous 
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experiments.  This  situation  was  used  to  illustrate  the  continued  ability  of  the  controller  to 
adapt  and  improve  results  over  time  (until  the  gains  reach  there  equilibrium  values).  The 
data  capture  time  was  thus  increased  to  4-niinutes  and  the  reduction  increased  to  11.44 
dB  as  shown  in  Figure  43. 


Node  26,  Z-Axis,  2nd  Modal  Freq,  Piezo  #2 


Frequency  (Hz) 

Figure  43.  Exp.  1, 4-Minute  Data  Capture  Power  Spectrum 
The  data  capture  time  was  increased  to  5  and  6  minutes  in  the  next  two 
experiment  sets.  Representative  samples  showed  a  reduction  of  19.6  dB  and  21.3  dB 
respectively  and  the  power  spectra  are  shown  in  Figure  44  and  Figure  45. 
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Figure  45.  Exp.  7,  6-Minute  Data  Capture  Power  Spectrum 


Further  increases  in  data  capture  times  were  not  pursued  due  to  the  large  file  sizes 
and  data  manipulation  resources.  It  was  also  discovered  while  attempting  to  explain  the 
slow  convergence  rate  that  the  full  truss  model  included  in  ANSYS  showed  this  mode 
shape  as  a  nearly  linear  displacement  of  the  LPACT  in  the  x-direction.  This  made  the 
actual  truss  response  less  critical  and  the  experimental  set  was  concluded  with  a  set  of  6- 
minute  data  collection  runs.  A  representative  sample  of  the  controller  gain  values  and 
system  response  from  the  finial  iteration  show  that  the  gains  are  still  not  fully  converged 
indicating  an  improved  response  is  possible  given  longer  runs.  These  results  are  shown 
in  Figure  75  and  Figure  76  of  Appendix  G. 

9.  Experiment  8  Modal  Freq,  Z-axis  of  Node  26,  Piezo  #2) 

After  reexamining  the  results  form  the  modal  analysis,  it  was  determined  that  the 
third  modal  frequency  predicted  by  ANSYS  corresponded  closest  to  the  bare  truss  second 
frequency  and  had  the  second  greatest  effect  on  the  overall  motion  of  the  truss.  The 
largest  displacement  in  this  mode  was  in  the  z-direction.  The  predicted  frequency  of 
12.529  Hz  also  agreed  with  the  second  modal  frequency  determined  experimentally 
Figure  38.  The  initial  data  capture  time  was  once  again  3 -minutes,  but  the  data  showed 
that  the  gains  had  not  fully  reached  equilibrium  at  the  end  of  the  run  even  though  a 
reduction  of  28.92  dB  was  realized  at  the  disturbance  frequency.  Consequent  data  runs 
with  a  collection  time  of  4-minutes  showed  that  the  gains  essentially  reached  equilibrium 
with  a  total  reduction  of  56.61  dB.  Representative  results  are  shown  in  Figure  77,  Figure 
78,  and  Figure  79  of  Appendix  G. 
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10.  Experiment  9  (Controller  Comparison) 

The  last  set  of  experimental  conditions  was  chosen  to  reproduce  the  conditions 
that  resulted  in  the  greatest  vibration  reduction  using  the  IDIFF  Controller  on  the  same 
structure.  Although  the  same  active  element  was  used  and  the  same  disturbance 
frequency  was  input,  it  is  important  to  realize  that  the  design  intention  of  the  previous 
controller  was  to  actively  dampen  the  vibration  by  interrupting  the  most  significant 
conduction  path.  The  adaptive  multi-layer  LMS  controller  is  designed  with  the  express 
purpose  to  counter  any  induced  vibration  at  the  control  node  by  introducing  a  similar 
signal  that  is  out  of  phase  with  the  disturbance  at  that  point.  The  controllers  obviously 
have  different  apphcations,  but  are  compared  to  show  the  advantage  of  this  controller  if 
the  desired  result  is  to  reduce  a  single  axis  disturbance  at  a  given  node. 

The  best  results  shown  in  Reference  5  were  noted  on  the  z-axis  of  node  41  with  a 
13.81  Hz  disturbance  source.  The  active  damping  resulted  in  a  reduction  of  19.02  dB  in 
that  specific  node  [Ref  5].  The  same  disturbance  frequency  was  chosen  for  this 
experiment  and  the  control  signal  was  the  acceleration  from  the  z-axis  of  node  41.  The 
same  active  element  was  also  used  for  sake  of  comparison.  The  spectral  analyses  of  the 
old  and  new  controllers  are  shown  in  Figure  46  and  Figure  47  respectively.  It  was  also 
observed  that  the  controller  converged  very  quickly  in  this  case  and  reduced  the 
disturbance  down  to  the  noise  level  for  a  reduction  of  57.26  dB.  Based  on  these  results,  it 
is  obvious  that  controller  location  is  near  optimal  for  that  specific  node  and  axis  which 
also  explains  the  performance  of  the  active  damping  system.  The  gain  and  system 
responses  are  shown  in  Figure  80  and  Figure  81  of  Appendix  G. 
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Figure  46.  ffilFF  Controller  Results  [From  Ref.  5] 


Node  41 ,  Z-Axis,  13.81  Hz,  Piezo  #1 


11.  Unstable  Configuration 


As  noted  in  section  [IILB.  4],  there  was  one  unstable  system  configuration 
revealed  after  the  controller  development.  The  discovery  was  made  during  a 
demonstration  of  the  system  under  random  conditions  chosen  by  the  observer.  Ironically, 
the  observer  varied  only  one  of  three  of  the  demonstration  default  values  choosing  the  y- 
axis  of  node  41,  with  a  14  Hz  disturbance  signal.  The  exact  cause  of  the  instability  was 
not  determined,  but  it  is  suspected  that  it  could  be  removed  with  further  adjustment  of  the 
individual  control  element  adaptation  rates. 

As  noted  in  the  analytical  section  this  algorithm  is  theoretically  capable  of 
guaranteed  convergence,  but  not  guaranteed  effectiveness.  The  addition  of  multiple 
adaptation  rates  was  necessary  to  increase  the  applicability  of  this  controller  to  random 
initial  conditions,  but  sacrificed  an  order  of  magnitude  on  the  potential  convergence  time 
of  the  phase  adjust  LMS  algorithm.  This  effect  is  obvious  by  comparing  the  required 
time  for  convergence  on  the  experimental  and  analytical  models  [Appendix  I]. 
Increasing  the  adaptation  rate  of  the  gain  magnitude  LMS  algorithm  minimized  this 
effect,  but  is  also  the  likely  source  of  this  instability.  Additional  testing  is  required  to 
evaluate,  document,  and  correct  this  situation  if  possible. 

D.  ANSYS  FEM  COMPARISON 

1.  Model  Setup  and  Initiation 

A  modal  analysis  was  done  xising  ANSYS  and  the  natural  firequencies  obtained 
were  significantly  different  fi'om  the  original  bare  truss  modal  firequencies.  This  was 
contributed  to  the  added  mass  of  the  LPACT  and  control  apparatus  (~20%  of  the  bare 
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truss  mass).  The  first  four  natural  firequencies  obtained  fi-om  the  modified  FEM  are 
compared  with  the  experimental  fi-equencies  obtained  fi-om  the  DSA  [Ref.  5]  in  Table  5. 


Mode 

ANSYS 

Frequency 

DSA 

Frequency 

1 

9.43 

8.94 

2 

11.32 

12.19 

3 

12.53 

13.87 

4 

18.74 

17.50 

Table  5.  Natural  Frequencies  of  Truss  With  Active  Components 
Upon  examining  the  mode  shapes  firom  the  postprocessor,  it  was  determined  that 
the  second  natural  fi-equency  was  essentially  due  to  the  2.5  kg  mass  of  the  LPACT 
vibrating  almost  entirely  in  the  x-direction.  If  the  model  is  assumed  to  be  accurate,  then 
this  mode  is  also  nearly  entirely  contained  in  that  single  element  and  thus  not  controllable 
by  the  current  system.  Due  to  the  position  of  the  LPACT  and  the  orientation  of  the  force 
axis,  this  fi-equency  can  be  detected  by  the  system  and  subsequently  controlled,  however 
the  control  authority  of  the  current  configuration  is  minimal  and  thus  leads  directly  to  the 
long  convergence  times  and  reduced  effectiveness  demonstrated  in  the  experimentation 
described  in  sections  [V.C.  7]  and  [V.C.  8]  which  describe  experiments  six  and  seven 
respectively. 

One  factor  that  was  not  included  in  the  ANSYS  model  was  the  active  element 
preloaded  required  on  the  actual  truss.  The  preload  was  done  using  shims  to  give  a 
compressive  bias  for  the  piezoelectric  actuator  and  then  the  control  signal  was  added  to  a 
50V  DC  signal  that  centered  the  piezoelectric  element  in  the  center  of  its  operating  band. 
In  the  ANSYS  model,  this  effect  was  necessary  due  to  the  mathematical  modeling  of  the 
SOLID5  elements,  but  a  limiter  was  included  to  prevent  the  control  signal  firom 
exceeding  +/-  50V. 
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The  simulations  were  run  for  times  ranging  from  5  to  50  seconds  of  simulated 
time.  The  time  was  varied  based  on  capturing  the  optimal  performance  of  the  controller 
based  on  the  convergence  rate  obtained  for  each  set  of  experimental  conditions.  The 
longest  data  run  (50  seconds)  was  extended  to  gather  enough  data  points  for  a  power 
spectrum  analysis  of  the  system  before  and  after  control.  Other  power  spectra  shown 
were  generated  with  fewer  data  points,  but  were  considered  adequate  for  comparison. 
The  time-step  interval  was  set  at  1/20*  the  driving  frequency. 

The  active  control  APDL  macro,  “TRUSS_CONTROL.TXT”  [Appendix  E],  was 
used  to  vary  the  experimental  parameters  to  match  those  used  in  truss  experiments.  The 
time  required  for  each  data  run  varied  from  approximately  15  minutes  to  several  hours 
based  on  the  niunber  of  iterations  required  by  the  combination  of  simulation  time  and 
time-step  required.  The  files  saved  by  the  overarching  APDL  macro  were  ASCII  files 
that  contained  the  control  node  displacement,  velocity,  and  acceleration  as  well  as  the 
actual  control  signal  and  controller  gain  values.  A  MATLAB  algorithm  was  developed  to 
process  and  display  the  data,  included  as  Appendix  J.  A  second  algorithm  was  \Mitten  to 
allow  for  comparison  of  the  experimental  and  FEM  results,  included  as  Appendix  K. 

The  different  nodal  numbering  and  axis  definition  used  in  the  model  had  to  be 
accounted  for  in  the  APDL  file  when  duplicating  the  experimental  conditions.  The 
changes  made  include  switching  the  Y-axis  and  Z-axis  and  verifying  the  node  numbering. 
As  previously  discussed,  the  adaptive  gain  adjustment  was  not  included,  so  an  average  of 
3  abbreviated  data  runs  were  required  to  determine  the  proper  gain  to  be  applied  to  the 
controller  for  each  set  of  experimental  conditions.  The  simulated  convergence  rates  were 
significantly  faster  due  to  the  near  optimal  gains  used  and  the  removal  of  the  layered 
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adaptive  nature  of  the  experimental  control  system.  This  was  expected  and  is 
representative  of  the  possibilities  if  a  simplified  controller  were  to  be  applied  and  tuned  to 
specific  situations.  The  equilibrium  results  in  either  the  multi-layered  controller  or  the 
tuned  controller  are  the  same,  so  the  comparison  is  still  a  valid  verification  of  the  model 
response. 

The  results  of  each  data  run  are  included  in  Appendix  L,  and  comparisons  of  the 
results  are  either  shown  in  the  text  or  in  Appendix  1. 

2.  Simulation  1  (19  Hz,  Y-axis  of  Node  18,  Piezo  #2) 

In  this  data  run,  the  conditions  of  the  first  experiment  were  duplicated  with  a  19 
Hz  disturbance,  controlling  the  y-axis  of  node  18,  using  element  number  two.  To  account 
for  the  variation  in  nodal  numbering  and  axis  definition,  the  APDL  file  was  written  to 
control  the  z-axis  of  node  49,  which  are  the  model  equivalents  to  the  experimental 
parameters.  The  yalue  used  for  the  controller  gain  (p  in  the  LMS  algorithm)  in  this  run 
was  —0.045  and  a  5  second  acquisition  time  was  sufficient  to  determine  system  response. 

Simulation  results  of  control  node  response,  control  signal,  and  adaptive  LMS 
gains  are  shown  in  Figure  98,  Figure  99,  and  Figure  100  of  Appendix  L.  A  comparison 
of  the  nodal  acceleration  is  shown  in  Figure  48. 
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Figure  48.  Case  1,  Control  Node  Acceleration  Comparison 


No  attempt  was  made  to  intentionally  match  the  initial  magnitude  of  the  nodal 
accelerations  felt  by  the  control  node.  This  was  not  deemed  necessary  due  to  the  adaptive 
nature  of  the  controller,  but  it  does  reflect  in  the  magnitude  of  the  resulting  control  signal 
and  the  rate  of  convergence.  While  the  rate  of  convergence  is  typically  a  critical  value, 
the  process  of  manually  determining  the  gain  for  the  simulations  removed  the  effect  of 
magnitude  on  the  convergence  rate.  If  the  entire  control  system  were  to  be  modeled,  the 
signal  strength  would  have  to  be  match  to  predict  the  control  system  response  correctly. 
The  control  signals  generated  by  both  the  experimental  and  simulated  cases  are  shown  in 
Figure  49. 
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Node  18,  Y-Axis.  19  Hz,  Piezo  #2 


Time  (sec) 


Figure  49.  Case  1,  Control  Signal  Comparison 

The  finial  reduction  in  the  nodal  acceleration  predicted  by  the  FEM  was  68.43  dB. 
The  difference  between  this  reduction  and  the  47.38  dB  reduction  realized  in  the  first 
experiment  was  due  to  the  difference  in  noise  level.  No  signal  noise  was  introduced  into 
the  FEM,  so  the  noise  level  was  effectively  due  to  the  mathematical  precision  of  the 
computations.  This  difference  supports  the  assumption  that  the  control  law  is  capable  of 
reducing  a  specific  frequency  disturbance  to  the  noise  level,  regardless  of  initial 
magmtude,  if  the  system  has  sufficient  control  authority. 

3.  Simulation  2  (15  Hz,  Z-axis  of  Node  41,  Piezo  #1) 

The  second  simulation  corresponds  to  the  conditions  of  the  second  experiment, 
using  a  15  Hz  disturbance,  controlling  the  z-axis  of  node  41,  using  element  number  one. 
The  APDL  file  was  written  to  control  the  y-axis  of  node  26,  which  are  the  model 
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equivalents  to  the  experimental  parameters.  The  value  used  for  the  controller  gain  (p  in 
the  LMS  algorithm)  in  this  run  was  slightly  higher  than  the  first  simulation,  at  -0.075  but 
a  5  second  acquisition  time  was  still  sufficient  to  determine  system  response. 

The  simulation  results  are  shown  in  Figme  101,  Figure  102,  and  Figure  103  of 
Appendix  L.  A  comparison  of  these  results  with  the  experimental  data  is  shown  in  Figure 
82  and  Figure  83  of  Appendix  I.  Due  to  the  very  rapid  convergence  of  the  LMS  gains  in 
this  simulation,  the  reduction  in  vibration  was  greatly  enhanced  as  the  controller  was 
allotted  more  time  to  “fine  tune”  the  response  and  drive  the  signal  down  to  the 
mathematical  noise  of  the  computational  algorithm.  The  finial  reduction  shown  in  the 
simulation  was  175.28  dB.  The  results  of  experiment  number  two  showed  a  42.45  dB 
reduction.  Both  are  essentially  the  difference  between  the  original  signal  and  the 
applicable  system  noise. 

4.  Simulation  3  (12  Hz,  X-axis  of  Node  26,  Piezo  #2) 

The  third  simulation  corresponds  to  the  conditions  of  the  third  experiment,  using  a 
12  Hz  disturbance,  controlling  the  x-axis  of  node  26,  using  element  number  two.  The 
APDL  file  was  written  to  control  the  x-axis  of  node  41,  which  are  the  model  equivalents 
to  the  experimental  parameters.  The  value  used  for  the  controller  gain  (p  in  the  LMS 
algorithm)  in  this  run  was  the  largest  of  any  of  the  simulations,  at  0.9,  but  a  5  second 
acquisition  time  was  sufficient  to  determine  system  response.  The  large  gain  required  for 
rapid  convergence  is  an  indication  that  the  signal  strength  for  the  disturbance  was  weak. 
It  was  also  necessary  to  reduce  the  magnitude  of  the  disturbance  signal  by  a  factor  of  in 
order  to  keep  the  disturbance  fi'om  exceeding  the  control  authority  of  the  active  element. 
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This  reduction  kept  the  control  signal  to  within  the  +/-  50V  operational  limit.  The 
amount  of  controller  input  required  to  minimize  this  nodal  acceleration  is  consistent  with 
the  node  shape  and  control  element  orientation,  however  inconsistent  with  the 
experimental  results.  This  is  credited  to  the  orientation  of  the  LPACT  perpendicular  to 
the  control  axis,  thus  minimizing  the  disturbance  energy  at  the  real  control  node  and 
increasing  the  signal  to  noise  ratio  as  noted  in  the  experimental  results  section. 

The  simulation  results  are  shown  in  Figure  104,  Figure  105,  and  Figure  106  of 
Appendix  L.  A  comparison  of  these  results  with  the  experimental  data  is  shown  in  Figure 
84  and  Figure  85  of  Appendix  I.  The  finial  reduction  in  the  simulation  was  73.39  dB  and 
the  experimental  reduction  was  23.07  dB. 

5.  Simulation  4  (12  Hz,  Y-axis  of  Node  26,  Piezo  #1) 

As  before,  this  simulation  mimics  the  conditions  of  the  corresponding  experiment, 
using  a  12  Hz  disturbance,  controlling  the  y-axis  of  node  26,  using  element  number  one. 
The  APDL  file  was  wntten  to  control  the  z-axis  of  node  41,  which  are  the  model 
equivalents  to  the  experimental  parameters.  The  value  used  for  the  controller  gain  (p  in 
the  LMS  algorithm)  in  this  run  was  at  0.005,  which  is  the  smallest  value  used.  The  5- 
second  acquisition  time  was  not  sufficient  in  this  case  to  determine  system  response 
because  of  a  relatively  slow  convergence  rate.  The  simulation  time  was  extended  to  45 
seconds  and  the  controller  gains  were  near,  but  not  fully  converged  to  equilibrium  at  the 
end  of  the  simulation.  The  significantly  reduced  gain  suggested  that  the  magnitude  of  the 
disturbance  in  that  axis  was  larger  than  the  previous  simulation.  That  was  consistent  with 
the  mode  shape  and  the  measurements  taken  on  the  actual  truss.  The  magnitude  of  the 
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disturbance  signal  again  had  to  be  reduced,  this  time  by  a  factor  of  V3  the  original  value 
in  order  to  keep  the  distxirbance  from  exceeding  the  +/-  50V  operational  limit  of  the 
actual  active  element.  Both  the  experimental  and  simulated  LMS  gain  factors 
approached  their  equilibrium  values  in  an  oscillatory  manor  as  seen  in  Figure  87. 

The  simulation  results  are  shown  in  Figure  107,  Figure  108,  and  Figure  109  of 
Appendix  L.  A  comparison  of  these  results  with  the  experimental  data  is  shown  in  Figure 
86  and  Figure  87  of  Appendix  I.  The  finial  reduction  in  the  simulation  was  26.06  dB  and 
the  experimental  reduction  was  25.24  dB. 

6.  Simulation  5  (1**  Modal  Freq,  Y-axis  of  Node  26,  Piezo  #1) 

This  simulation  possible  explains  the  nonlinearities  seen  in  experiment  number 
five.  The  same  conditions  were  input  into  the  model  (9.4265  Hz  disturbance,  controlling 
the  y-axis  of  node  26,  using  Piezo  #1),  but  no  nonlinearities  were  noted  either  before  or 
after  the  controller  was  activated.  It  is  theorized  that  this  is  primarily  due  to  the  FEM 
model  of  the  LPACT  being  a  rigid  part  of  an  element.  The  failure  of  the  harmonics 
frequencies  to  show  up  in  a  model  that  does  not  include  the  LPACT  spring-mass 
interactions  and  the  amplification  of  the  second  harmonic  seen  in  the  initial  power 
spectrum  of  experiment  number  five  bolh  support  this  theory.  As  with  previous 
simulations  the  APDL  file  was  written  to  the  FEM  equivalent  parameters,  controlling  z- 
axis  of  node  41.  The  value  used  for  the  controller  gain  (p  in  the  LMS  algorithm)  in  this 
run  was  at  -0.08  and  the  simulation  time  was  set  to  10  seconds. 

The  simulation  results  are  shown  in  Figure  110,  Figure  111,  and  Figure  112  of 
Appendix  L.  A  comparison  of  these  results  with  the  experimental  data  is  shown  in  Figure 
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88  and  Figure  89  of  Appendix  I.  The  finial  reduction  in  the  simulation  was  54.85  dB  and 
the  experimental  reduction  was  18.80  dB.  The  difference  is  contributed  to  the  absence  of 
harmonics  in  the  model  and  the  excitation  of  the  disturbance  frequency  by  the  first 
harmonic  in  the  experimental  case. 

7.  Simulation  6  (2"**  Modal  Freq,  Y-axis  of  Node  26,  Piezo  #2) 

As  noted  in  the  experimental  version  of  this  simulation,  the  mode  shape  of  this 
modal  frequency  is  nearly  entirely  described  by  displacement  of  the  LPACT  mass  in  the 
x-direction.  The  simulation  was  run  on  the  ANSYS  model  under  the  same  conditions  as 
experiment  six,  to  test  the  model’s  ability  to  predict  the  controller  response  in  a  less  than 
ideal  control  situation.  The  disturbance  frequency  was  set  to  11.318  Hz  and  element 
number  two  was  used  to  control  the  y-axis  of  node  26.  The  APDL  used  the  model 
equivalents  of  the  z-axis  of  node  41.  The  controller  gain  (p  in  the  LMS  algorithm)  was 
determined  to  be  -0.015  and  the  simulation  time  was  extended  to  45  seconds  to  allow  for 
convergence. 

The  model  again  failed  to  predict  the  nonlinearities  seen  in  the  experimental  case, 
but  it  did  exhibit  the  same  slow  convergence  and  reduced  effectiveness.  These  results 
seem  to  indicate  that  the  harmonics  do  not  directly  affect  the  controller  performance  in 
this  case.  That  may  be  true  for  the  previous  experiment  as  well,  which  would  mean  the 
only  sigmficant  error  on  that  evolution  was  the  failure  to  model  the  mass-spring 
relationship  of  the  LPACT.  More  testing  is  required  to  either  confirm^  or  deny  the 
discrepancy. 
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The  simulation  results  are  shown  in  Figure  113,  Figure  114,  and  Figure  115  of 
Appendix  L.  A  comparison  of  these  results  with  the  experimental  data  is  shown  in  Figure 
90  and  Figure  91  of  Appendix  I.  The  finial  reduction  in  the  simulation  was  32.42  dB  and 
the  experimental  reduction  was  41.54  dB. 

8.  Simulation  7  (2"**  Modal  Freq,  Z-axis  of  Node  26,  Piezo  #2) 

The  simulation  of  these  experimental  conditions  showed  the  same  relatively  slow 
convergence  rate  as  seen  in  experiment  seven.  This  slow  convergence  could  not  be 
corrected  by  adjustments  to  the  gain,  even  though  multiple  attempts  were  made.  The 
fastest  convergence  was  an  oscillatory  approach  to  gain  convergence.  This  approach  was 
normally  avoided  due  to  the  obvious  proximity  to  instability.  However,  for  the  FEM 
simulation  it  was  desirable  to  force  the  system  to  equilibrium  and  the  absence  of  noise 
made  the  transition  to  instability  predictable  based  on  the  graphs  of  the  gain  over  time. 
Even  with  the  optimal  gains  and  a  45  second  data  capture  simulation  time,  the  gains  were 
still  not  fully  converged.  The  simulation  was  stopped  at  this  time  because  similarity 
between  system  and  model  response  had  been  verified. 

The  simulation  introduced  an  11.318  Hz  disturbance  and  controlled  the  z-axis  of 
node  26,  using  element  niunber  two.  The  APDL  file  was  written  to  control  the  y-axis  of 
node  41,  which  are  the  model  equivalents  to  the  experimental  parameters.  The  value 
used  for  the  controller  gain  (p  in  the  LMS  algorithm)  in  this  run  was  at  -0.005. 

The  simulation  results  are  shown  in  Figure  116,  Figure  117,  and  Figure  118  of 
Appendix  L.  A  comparison  of  these  results  with  the  experimental  data  is  shown  in  Figure 
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92  and  Figure  93  of  Appendix  1.  The  finial  reduction  in  the  simulation  was  18.29  dB  and 
the  best  experimental  reduction  was  21.26  dB. 

9.  Simulation  8  Modal  Freq,  Z-axis  of  Node  26,  Piezo  #2) 

The  last  set  of  conditions  used  to  test  the  simulation  was  also  used  to  generate 
power  spectra  both  before  and  after  model  control.  The  previous  power  reductions  were 
calculated  based  on  smnlar  spectra  that  were  calculated  with  significantly  reduced 
resolution.  The  spectra  generated  in  this  simulation  are  included  in  this  thesis  and 
required  an  expended  simulation  time  of  50  seconds  for  data  capture.  The  exact  time  for 

the  simulation  is  not  known  because  it  was  run  overnight,  but  it  was  in  excess  of  3  hours 
in  real  time. 

The  conditions  of  the  8*  experiment  were  duplicated  with  a  12.529  Hz 
disturbance,  controlling  the  z-axis  of  node  26,  using  element  nrnnber  two.  To  accoimt  for 
the  variation  in  nodal  numbering  and  axis  definition,  the  APDL  file  was  written  to  control 
the  y-axis  of  node  49,  which  are  the  model  equivalents  to  the  experimental  parameters. 
The  value  used  for  the  controller  gain  (p  in  the  LMS  algorithm)  in  this  run  was  0.005. 
The  magnitude  of  the  disturbance  source  was  reduced  to  25%  of  the  full  strength 
magmtude  in  order  to  keep  the  controller  within  operating  limits.  This  turned  out  to  be 
over  conservative  and  the  control  signal  stayed  below  +/-  20  volts  (+/-  50  V  max). 

The  power  spectrum  of  the  FEM  during  this  simulation  shows  that  it  did  predict 
the  second  tonal  frequency,  but  it  failed  to  show  the  first.  The  difference  in  noise  level 
between  the  two  spectra  is  attributed  to  the  linear  approximation  (reverse  difference 
approximation)  used  in  the  calculation  of  the  nodal  velocities  and  accelerations  based  on 
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FEM  calculated  displacements.  The  plot  shows  the  power  spectrum  of  the  nodal 
acceleration,  which  is  used  as  the  input  to  the  controller  in  both  the  real  and  simulated 
experiments. 


Node  26,  Z-Axis,  3rd  Modal  Freq,  Piezo  #2 


Figure  50.  Simulation  8,  Power  Spectrum  Before  &  After  Control 
The  magnitude  of  the  reduction  exhibited  in  this  simulation  was  due  mostly  to  the 
initial  magnitude  of  the  disturbance  and  the  length  of  time  allowed  for  convergence. 
From  Figure  50,  it  appears  that  another  70  dB  of  reduction  may  have  been  realized 
(assuming  that  the  mathematical  floor  noise  level  had  stabilized)  if  the  controller  was 
allowed  to  continue  adapting.  This  is  not  the  case  because  of  the  nature  of  the  adaptive 
LMS.  The  controller  gains  will  converge,  but  will  never  stop  “hunting”  for  the  correct 
value.  This  introduces  slight  variation  in  the  control  gain  that  will  always  prevent  a 
perfect  cancellation,  if  any  signal  noise  is  present.  The  difference  between  a  reduction  of 
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112  dB  and  1 82  dB  in  the  power  of  an  actual  vibration  with  an  initial  power  on  the  order 
of  20  dB  is  not  considered  significant  for  most  applications.  The  simulated  system 
response  and  LMS  gain  values  for  this  simulation  are  shown  in  Figure  51  and  Figure  52 
respectively. 


Time  (sec) 


Figure  51.  Simulation  8,  System  Response 
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Node  26,  Z-Axis,  3rd  Modal  Freq.  Piezo  #2 
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Figure  52.  Simulation  8,  Controller  Gains 
A  comparison  of  these  results  with  the  experimental  data  is  shown  in  Figure  94, 

Figure  95,  Figure  96,  and  Figure  97  of  Appendix  1.  The  finial  reduction  in  the  simulation 

was  1 12.3  dB  and  the  best  experimental  reduction  was  56.61  dB. 

10.  Comparison  and  Summary  of  Results 

A  summary  of  the  different  cases  is  presented  in  Table  6.  These  figures 
demonstrate  that  the  controller  is  able  to  adapt  to  various  frequencies  and  structural 
configurations.  The  validation  of  the  FEM  should  include  similarity  of  the  system 
responses  shown  in  Appendix  I  and  not  just  the  vibration  reductions  stated  in  this  table. 
The  adaptive  nature  of  this  controller  makes  the  reduction  in  vibration  a  function  of  the 
signal-to-noise  ratio,  controller  gain,  and  the  time  allotted  for  controller  gains  con\  crgc. 
It  may  be  desirable  in  future  research  to  prove  that  the  reductions  (both  analjlical  and 


experimental)  converge  to  near  the  signal-to-noise  ratio  as  the  convergence  time  goes  to 
infinity. 


Case 

Number 

Disturbance 
Frequency 
_ (Hz) 

Control 

Element 

Control 

Node 

(#/axis) 

Analytical 

Gain 

Analytical 

Reduction 

(dB) 

Experimental 

Reduction 

(dB) 

1 

19 

2 

18  A" 

-0.045 

68.43 

47.38 

2 

15 

1 

41/Z 

-0.075 

175.28 

42.45 

3 

12 

2 

26/X 

0.9 

73.39 

23.07 

4 

12 

I  1 

26fY 

0.005 

26.06 

25.24 

5 

9.4265 

1 

^lefY 

-0.08 

54.85 

18.80 

6 

11.318 

2 

im 

-0.015 

32.42 

41.54 

7 

11.318 

2 

26/Z 

-0.005 

18.29 

21.26 

8 

12.529 

2 

26/Z 

0.005 

112.3 

56.61 

9 

13.81 

1 

41/Z 

N/A 

N/A 

57.26 

Table  6.  Summary  of  Results 


The  option  to  tune  to  the  ANSYS  model  of  the  controller  also  demonstrated  the 
controller’s  ability  to  converge  faster  if  the  application  is  specified  and  fi-equency 
specific.  This  was  not  expressly  demonstrated  in  the  experimental  version,  but  was 
observed  in  the  early  stages  of  controller  development.  The  addition  of  the  automatic 


fi'equency  identification  and  adaptive  gain  features  were  designed  to  improve  versatility 


at  the  expense  of  convergence  rate,  without  affecting  the  steady  state  performance  of  the 
controller. 

During  the  development  of  the  Adaptive  Multi-Layer  LMS  Controller,  hundreds 
of  experiments  were  conducted  on  the  NPS  Space  Truss.  Using  this  experience  made  the 
optimal  gains  that  were  manually  determined  in  the  ANSYS  model  a  relatively  simple 
task.  That  was  because  the  empirical  response  of  the  FEM  to  variations  in  disturbance 
fi'equency  and  magmtude  had  the  same  affect  on  the  required  gain  as  the  actual  truss 
system.  To  quantify  the  accuracy  of  its  prediction  of  LMS  gain  values,  further 


98 


experimentation  must  be  performed  which  matches  the  analytical  disturbance  force  and 
uses  a  more  precise  method  of  calculating  nodal  acceleration  in  the  model.  There  would 
also  have  to  be  representative  noise  added  to  the  FEM  signals  used  as  controller  input  and 
output. 
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VI. 


FUTURE  RESEARCH 


A.  IMPACT  ANALYSIS 

1.  Background 

High-speed  motion  photography  can  capture  the  structural  response  of  impacts 
such  as  a  bullet  hitting  a  flack  jacket.  Those  types  of  images  motivated  this  approach  to 
impact  analysis.  They  clearly  capture  a  structoal  response  to  the  impact  besides  the 
axially  transmission  of  energy  through  the  material.  The  difference  between  these  forms 
of  propagation  allows  them  to  be  differentiated  and  utilized  independently. 

An  analysis  of  structural  impacts  using  installed  accelerometers  should  be 
possible  using  cross  correlations  with  known  reference  signals.  The  analysis  output 
would  include  impact  location,  energy,  and  angle.  This  would  provide  the  ability  to 
rapidly  evaluate  the  potential  damage  done  by  unexpected  impacts. 

2.  Analysis  Technique 

The  reference  signals  required  for  this  analysis  are  three  orthogonal  impacts  at 
each  sensor.  Two  variations  of  each  reference  impact  should  saved,  one  normalized  with 
respect  to  the  energy  of  impact  and  the  other  with  respect  to  the  peak  response.  These 
signals  are  representative  of  each  nodes  response  to  a  unit  impact  in  a  specific  direction. 
The  peak-normalized  data  is  to  be  used  to  evaluate  the  location  and  angle  of  impact  by 
comparing  the  characteristic  shape  of  the  responses  and  the  delay  times.  The  impact- 
normalized  signals  are  used  to  determine  the  magnitude  of  the  impact. 


101 


The  analysis  can  take  two  different  approaches  based  on  sample  frequency.  The 
material  response  is  characterized  by  the  transmission  of  relatively  high  frequency 
vibration  through  the  structure.  In  the  NPS  Space  Truss,  this  energy  is  transmitted  axially 
through  every  element  of  the  truss.  This  correctly  implies  multiple  paths  and  high-speed 
propagation.  The  rate  of  propagation  is  estimated  by  using  the  material  properties  of  the 
structure.  For  the  NPS  Space  Truss,  this  is  Aluminum. 


Based  on  this  propagation  rate,  a  delay  time  can  be  estimated  based  on  the  size  of 
the  structure  in  order  to  access  the  required  sample  time  to  capture  the  signal.  For  this 
estimation,  a  straight-line  path  is  used  and  structural  reflections  are  ignored.  The 
maximum  dimension  of  the  truss  is  3.67  meters  and  the  distance  between  nodes  is  0.33 
meters.  The  sensors  are  only  attached  to  nodes,  thus: 


Delay^ 


3.67m 

5,100"/ 


■  0.7  msec 


n  ;  0.33m  ^  ' 

node-node  =  T  ,  /  «  0.06  mSCC 

5,100"/ 

/  sec 

Assuming  that  10  samples  are  required  for  resolution  between  delay  times,  this 
leads  to  required  sample  rates  on  die  order  of: 
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While  obtainable,  these  sample  rates  are  higher  than  those  that  are  widely 
available  on  current  spacecraft.  It  is  therefore  desirable  to  use  the  other  method  of  energy 
transmittance  through  the  structure  in  order  to  reduce  the  required  sample  rates.  The 
natural  frequency  and  mode  shapes  of  the  structure  as  a  whole  dominate  the  structural 
transmittance.  The  propagation  rates  of  this  energy  are  estimated  based  on  the  mode 
shape  of  the  structure.  For  this  order  of  magnitude  estimation,  it  is  assizmed  that  the  NFS 
Space  Truss  responds  as  a  bare  truss  and  the  mode  shapes  predicted  by  both  MATLAB 
and  ANSYS  FEM  analyses  are  used  as  a  basis.  From  those  mode  shapes,  it  is  assumed 
that  the  path  length  of  the  1®‘  and  2”*^  modes  is  approximately  4  times  the  length  of  the 
truss.  Similarly,  the  3'^^  and  4*  modes  have  path  lengths  of  about  14  that  distance  or  2 
times  the  truss  length.  These  estimations  are  based  on  Figure  53  and  Figure  54. 


Mode  1 


‘  •  ■  ■  •  • 


.’T-K 


•  •  . 


Figure  53.  Extrapolation  of  Mode  1 


Figrure  54.  Extrapolation  of  Mode  3 
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The  natural  frequencies  of  the  first  two  modes  are  both  about  10  Hz  and  the 
second  two  modes  are  about  20  Hz.  These  numbers  are  expected  due  to  the  synunetry  of 
^'Shaped,  NPS  Space  Truss.  Using  the  nodal  frequencies  and  the  extrapolated  path 
length  to  estimate  propagation  rates  results  in  the  following  calculations: 

V  =  (Freq)(Wave  Length) 

^-2  ~  (yinode  =  (10  Hz)(4)(3.67  m) » 150  ^  (6.4) 

^3-4  =  {f^2){2L)  =  (20  Hz)(2)(3.67  m) « 150  V 

/  sec 

Based  on  this  much  slower  propagation  rate,  a  new  delay  time  can  be  estimated 

using  the  same  basis  as  before.  The  delays  of  structurally  transmitted  energy  are  as 
follows: 


r>  /  3.67m 

end-end  =  - - 7“  «  25  HlSeC 


Delay^ 


150™/ 

/  sec 

0.33m 


(6.5) 


node-node 


150 


m/ 


■2  msec 


sec 


The  longer  delay  times  reduce  the  required  sample  rate  to: 

„  10  samples 

^mend-end  =  ^  »  400  Hz 

25  msec 

E,  10  samples  ^ 

P^^qnode-node  =  “ ^  «  5  kHZ 
2  msec 


(6.6) 


These  rates  may  be  able  to  be  further  reduced  if  the  required  resolution  can  be 
obtained  with  fewer  than  10  samples. 


3.  Application 

There  are  four,  3-axis  accelerometers  and  two  single  axis  accelerometers  mounted 
on  the  NPS  Space  Truss.  The  single  axis  sensors  are  positioned  at  an  angle  such  that  x,  y. 
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and  z-axis  accelerations  are  all  measured.  The  cross  coirelations  are  done  using  the  peak 
normalized  impact  signals  from  each  of  the  sensors.  This  corresponds  to  the  peak 
normalized  reference  signals  from  known  impacts  in  an  attempt  to  compare  only  the 
shape  of  the  response.  The  magnitude  is  determined  in  a  different  process. 

The  cross  correlations  contain  information  on  both  the  relative  delay  of  the  signals 
and  there  approximation  of  the  reference  signal.  The  delays  are  required  in  order  to 
determine  the  location  of  impact  and  the  quality  of  the  comparison  is  used  to  determine 
the  direction  of  impact.  The  “sharper”  the  cross  correlation,  the  better  the  fit  an  impact 
signal  is  with  the  reference.  Several  methods  were  attempted  to  quantify  this 
“sharpness”,  but  none  provided  consistent  results.  Assuming  that  this  can  be 
accomplished,  the  relative  fit  can  be  used  to  determine  the  approximate  angle  of  impact 
by  compiling  a  unit  vector  of  those  quantities. 

The  delay  times  are  the  key  to  estimating  the  location  of  impact.  It  is  assumed 
that  the  system  has  a  preset  minimiun  sensitivity  level  set  to  trigger  the  data  collection. 
The  first  sensor  to  meet,  or  exceed  this  level  becomes  the  trigger  node.  The  distance 
from  the  impact  to  the  trigger  node  is  unknown.  The  delay  times  determined  from  the 
cross  correlations  from  the  other  sensors  measure  the  time  from  the  trigger  to  the  energy 
reaching  the  relative  sensors. 
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D„=^{x-x„y  +{y-y^f  +  {z-z„y  (6.8) 

The  distance  from  any  node  is  also  equal  to  the  distance  measured  using  the  delay 
time,  plus  the  distance  the  energy  traveled  before  reaching  the  trigger  node: 

A=^r+k|  (6.9) 

The  result  is  “n”  equations  with  the  impact  location  and  trigger  distance  as 
unknowns.  The  trigger  distance  (Dt)  is  a  constant  for  all  these  equations,  so  by 
minimizing  the  deviation  in  Dt,  you  localize  the  impact  location.  The  result  is  a  point 


that  corresponds  to  an  estimate  of  the  impact  location  (x,  y,  z).  The  larger  the  number  of 
sensors  used  for  the  localization,  the  smaller  the  error  (limited  by  the  error  in  measuring 
the  delay  times). 

Dr  =  sj{x-^.f  +  {y-y.f  +{z-z,f  (6.10) 

The  last  characteristic  of  the  impact  is  the  energy.  This  analysis  uses  only  the 
signal  with  the  best  correlation  to  a  known  impact.  A  method  of  evaluating  the  quality  of 
a  cross  correlation  is  suggested  to  quantify  this.  The  closest  signal  is  then  compared  to 
the  energy  normalized  reference  signal.  The  peak  response  of  the  raw  signal  is  divided 
by  the  energy  normalized  peak  response  to  give  an  estimation  of  the  actual  impact 
energy.  The  units  of  the  signal  compared  are  irrelevant  because  they  cancel  out. 


Impact  Energy  = 


Peak  Response 
Ref  Peak(energy  ■' ) 


(6.11) 


4.  Preliminary  Testing 

Initial  data  sets,  collected  using  only  4  of  the  sensors  and  a  sample  rate  of  1  kHz 
were  not  sufficient  to  produce  any  meaningful  results.  That  data  run  did  provide  the 
required  information  to  develop  the  information  discussed  above. 

The  next  data  set  was  taken  at  10  kHz,  in  an  attempt  to  capture  some  of  the 
structural  energy.  All  6  sensors  were  used  for  this  set,  resulting  in  15  channels  of  data 
being  captured  (14  firom  sensors  and  one  fi-om  the  impact  hammer).  The  volume  of  data 
limited  the  sample  time  to  only  0.2  seconds  per  impact.  This  size  can  be  increased,  but  at 
the  time  it  seemed  sufficient.  A  dozen  known  impacts  were  used  to  attempt  to  develop  a 
method  of  weighting  the  cross  correlations  in  such  a  manner  as  to  permit  the  required 
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calculations.  Using  various  methods  of  weighting  the  cross  correlations  of  this  data, 
results  varied  in  accuracy  from  <lm  to  3m  from  impact  location.  Based  on  the  size  of  the 
structure,  this  was  meaningless.  The  impact  angle  and  force  were  more  accurate.  Force 
was  used  instead  of  energy  due  the  ease  of  measuring  force.  Further  data  manipulation 
could  be  used  to  solve  for  impact  energy  after  the  method  is  established.  The  direction 
was  typically  accurate  to  less  than  30  in  spherical  coordinates  and  the  force  was  within 
15%.  The  results  prompted  another  set  of  experiments  in  an  attempt  to  isolate  the 

structural  response  and  clarify  the  resulting  resolution  of  the  cross  correlation  delay 
times. 

The  third  and  final  data  set  was  taken  from  all  6  sensors  at  1  kHz  sample  rate. 
The  sample  size  was  extended  to  capture  4  seconds  of  data  and  a  software  filter  was  built 
into  the  processing  to  allow  low  pass  filtering  at  various  comer  frequencies.  A 
Butterworth  filter  was  used  with  a  comer  frequency  of  30  Hz  in  order  to  capture  the  first 
4  modes  and  filter  out  the  high  frequency  noise.  The  delay  times  derived  from  this  setup 
were  inconsistent  and  the  suspected  causes  are  addressed  in  the  next  section. 

Further  research  into  this  question  is  precluded  due  to  a  lack  of  time,  facilitated  by 
the  early  transfer  of  the  principle  investigator. 

5.  Required  Analysis 

The  quality  of  the  cross  correlation  is  used  to  determine  the  angle  of  impact  and 
could  also  be  used  to  weight  the  accuracy  of  the  time  delay  signals.  The  methods  tried 
thus  far  have  been  inconsistent  with  respect  to  determining  relative  signal  .similarity,  The 
next  method  tried  should  be  an  initial  data  fit  to  a  normal  distribution  curve.  Because  all 
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signals  are  peak  nomalized  prior  to  correlation,  the  standard  deviation  of  the  distribution 
should  be  related  to  the  quality  of  the  signal  match. 

For  these  developmental  experiments,  the  impact  hammer  sensor  was  always  iwed 
to  trigger  the  data  collection  (both  reference  and  impact  signals).  This  should  have 
resulted  in  positive  delay  times  for  each  node,  and  a  zero  delay  if  the  impact  was  located 
at  the  sensor  node.  This  was  not  the  case  and  the  delay  times  were  often  negative.  The 
most  significant  of  the  errors  suspected  in  causing  this  is  the  method  of  collecting  the 
reference  data  sets.  A  sample  rate  equivalent  to  the  impact  sample  rate  was  always  used. 
This  seemed  logical  at  the  time,  but  introduced  an  unexpected  error.  The  time  of  impact 
has  to  be  captured  to  high  accuracy  in  order  to  be  a  useful  tool  in  correlation.  The 
hammer  trigger  was  observed  to  start  consistently  at  the  trigger  level  of  0.05  volts,  but 
often  at  the  down  side  of  the  peak.  This  resulted  in  delays  long  enough  to  cause  negative 
correlations  when  compared  to  actual  impact  signals.  To  correct  this,  the  trigger  level 
should  be  lowered  and  the  sample  time  of  the  reference  sets  increased  to  at  least  10  times 
the  impact  sample  rate  with  a  recommended  rate  of  150  kHz.  The  resulting  data  should 
be  down  sampled  to  correspond  to  the  sample  rates  used  for  analysis,  but  only  after  the 
initial  capture.  This  will  minimize  the  error  introduced  in  determining  the  time  delay. 

Another  problem  discovered  in  this  data  set  is  the  lack  of  resolution  between 
sensors  on  neighboring  nodes.  This  is  lead  to  and  is  consistent  with  the  assumption  that  a 
5  kHz  sample  rate  would  be  required  to  resolve  nodal  differences  using  structural 
response  [Equation  (6.6)]. 

The  ability  of  the  single  axis  accelerometers  to  resolve  the  angular  difference 
between  impacts  is  also  questionable.  The  technique  developed  to  determine  the 
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significance  of  the  cross  correlations  may  prove  this  a  mute  point  and  they  do  serve  as 
useful  data  points  in  determining  both  impact  location  and  energy,  regardless  of  angular 
resolution. 

B.  VIBRATION  SUPPRESSION 

1.  Simultaneous  Suppression  of  Multiple  Frequencies 

The  current  configuration  of  the  NPS  Space  Truss  includes  two  active  elements 
that  have  already  been  proven  to  be  capable  of  controlling  several  different  nodes  and 
axes.  This  significantly  reduces  the  implementation  requirements  of  developing  a 
suppression  system  capable  of  acting  on  two  ifrequencies  at  the  same  time.  Variations  on 
this  theme  include  controlling  either  the  same  node  or  different  nodes  and  axes. 

The  most  useful  application  of  this  process  would  also  include  a  fi-equency 
identification  technique  capable  of  high  resolution  (+/-  0.01  Hz)  and  amplitude 
classification.  This  would  allow  the  most  significant  disturbance  firequencies  to  be 
automatically  identified  and  suppressed. 

2.  Suppression  of  Multiple  Signals 

The  signal  used  to  generate  the  control  system  error  has  to  contain  both  amplitude 
and  phase  information.  There  is  no  requirement  that  the  signal  be  firom  a  single  sensor. 

It  may  be  possible  that  by  manipulating  the  signals  fi-om  various  sensors,  you  could 
develop  a  method  of  reducing  a  combination  of  disturbance  effects. 

There  have  been  rudimentary  attempts  to  mmimize  the  total  system  energy,  with 
no  significant  reduction.  This  outcome  is  predictable  due  to  the  inherent  nature  of 
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suppression  introducing  energy  into  the  system.  The  equilibrium  conditions  were 
consistently  low  amplitude  and  often  resulted  in  structural  effects  at  the  noise  level. 
Combinations  that  included  only  specific  nodes  and/or  axes  were  not  tried  but  should 
produce  measurable  reductions  related  to  the  specific  location  of  the  nodes,  firequency, 
control  element,  and  axis  chosen. 
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VIL 


CONCLUSIONS  AND  RECOMMENDATIONS 


The  Adaptive  Multi-Layer  LMS  Controller  has  been  shown  to  be  both  flexible 
and  capable  of  reducing  a  disturbance  signal  to  the  level  of  system  noise,  but  there  can  be 
significant  drawbacks  to  its  application.  If  the  control  element  is  placed  randomly  in  a 
structure  and  set  to  suppress  an  equally  random  point  and  axis,  the  controller  will  attempt 
to  apply  ANY  force  required  to  drive  the  nodes  response  to  the  disturbance  to  zero  (noise 
level).  The  actual  control  of  the  node  is  due  to  the  structural  amplification  of  the  force 
applied  by  the  control  element  (a  piezoelectric  element  in  this  case).  If  the  controller  has 
to  excite  the  x-axis  to  reduce  the  z-axis,  it  will.  This  can  result  in  significantly  increased 
off  axis  vibration  and  possible  an  increase  in  total  nodal  vibrational  energy. 

The  ANSYS  FEM  has  a  feature  that  allows  visualization  of  the  structural  response 
over  time,  during  actual  controller  application.  The  output  is  a  “.avi”  file  that  plays  like  a 
movie.  By  amplifying  the  magnitudes  of  the  displacements  by  5,000  times,  it  becomes 
clear  that  the  controller  uses  the  structural  response  to  the  control  signal  for  nodal  control. 
For  instance,  if  the  ay-axis  is  being  used  as  a  control  reference,  the  steady  state  condition 
will  show  that  the  controller  has  excited  either  the  x-axis  or  z-axis  to  coimter  the  LPACT 
disturbance.  This  technique  of  system  modeling  can  be  used  to  analyze  the  system 
response  prior  to  implementation  to  avoid  component  damage. 

The  best  application  of  the  suppression  system  would  be  a  single  axis  disturbance 
with  the  active  element  positioned  in  the  system  to  introduce  energy  into  that  specific 
axis  only.  The  next  best  thing  is  to  determine  the  best  physical  location  of  the  control 
element  to  minimize  the  energy  introduced  into  the  system  that  is  orthogonal  to  the  nodal 
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signal  being  minimized.  In  either  case,  the  relative  location  of  the  active  element  must  be 
optimized  with  respect  to  the  control  node  and  axis,  but  no  modification  to  the  control 
law  is  required.  If  the  active  element  location  is  optimized,  then  it  may  be  possible  to 
significantly  simplify  the  control  law. 
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APPENDIX  A.  NPS  SPACE  TRUSS  CHARACTERISTICS  AND  COMPONENTS 


Mass  Properties  of  the  E 

tare  and  Modified  Truss: 

Part  Name 

No.  In 
Bare  Truss 

No.  In 
Mod  Truss 

Component 

Masses  (kg) 

Node  Balls 

52 

52 

0.0663 

3.445 

3.445 

100 

100 

0.0448 

4.475 

4.475 

[51S?HiT5nMI 

61 

58 

0.0522 

3.181 

3.025 

LPACT  strut 

0 

1 

2.2760 

0.000 

2.276 

Active  Strut 

0 

2 

0.2900 

0.000 

0.580 

Screw 

322 

322 

0.0019 

0.607 

0.607 

Total  Mass 
Truss 

11.708 

14.408 

Table  7.  NPS  Space  Truss  Mass  Properties 
[After  Ref.  13  and  the  addition  of  a  second  active  strut.] 


Bare  Truss  Natural  Frequencies 


Mode  Number 

Modal  Testing 
[Ref.  13] 

MATLAB  FEM 

ANSYS 

1 

14.64 

14.13 

14.25 

2 

16.26 

15.44 

15.57 

3 

30.41 

28.72 

28.93 

4 

33.97 

32.04 

32.26 

5 

62.93 

60.23 

60.76 

6 

74.54 

72.24 

72.93 

7 

80.66 

79.71 

81.67 

8 

101.01 

97.41 

96.61 

9 

126.23 

120.21 

115.41 

10 

135.97 

129.68 

122.26 

Table  8.  NPS  Space  Truss  E 

tare  Natural  Frequencies 

Property 

Steel 

Outer  Radius 

3.975  mm 

Inner  Radius 

3.078  mm 

Solid 

Inertia 

■maiiisMi 

Cross-sectional  Area 

1.96856e-5  m^ 

4.96e-5  m^ 

Table  9.  Truss  Element  Properties 
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COMPONENTS 


Piezoelectric  Translation  Model  P-848-30  S/N 


Signal 

[Volts] 

Expansion 

[Microns] 

Contraction 

[Microns] 

Hysteresis 

[Microns] 

Hysteresis 

[Percent] 

0.00 

0.00 

0.00 

0.00 

0.00 

1.00 

3.48 

6.43 

2.94 

5.98 

2.00 

7.58 

12.55 

4.97 

10.12 

3.00 

12.16 

18.38 

6.21 

12.64 

4.00 

17.20 

23.86 

6.66 

13.56 

5.00 

22.53 

29.06 

6.53 

13.28 

6.00 

27.96 

33.90 

5.95 

12.10 

7.00 

33.44 

38.38 

4.94 

10.05 

8.00 

38.83 

42.43 

3.61 

7.34 

9.00 

44.05 

46.07 

2.02 

4.11 

10.00 

49.14 

49.14 

0.00 

0.00 

Table  10.  Expansion  and  Contraction  Data  for  Model  P-843.30  [ 

'From  Ref  12] 

Open  Loop  Travel  (O-IOOV) 

45  |xm  +/-  20% 

Closed  Loop  Travel 

45  pm 

Stiffiiess 

33  N/pm  +/-  20% 

Force  Generation  (Blocked) 

1500  N  +/-20% 

Push/Pull  Force  Capability 

800/300  N 

Torque  Limit  (at  tip) 

350  mNm 

Capacitance 

5.4pF  +/-  20% 

Dynamic  Operating  Current  Coefficient 

15  pA/(Hz-pm) 

Unloaded  Resonant  Frequency 

10  kHz  +/-  20% 

Operating  Temperature  Range 

-20  to  +80  °C 

Mass  (w/o  cables) 

53  g 

Length 

73  mm 

Table  11.  P-843.30  Operating  Characteristics  [Ref.  22] 
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Planning  Systems  Incorporated  LPACT  [Ref.  16]  S/N  CML-030-020-1 


1  Item 

Value 

■■■nnnsHB 

5.5  Ib/amp 

Max.  Current 

1  amp 

Coil  Resistance 

9  ohms 

Flexure  Natural  Frequency  (o„  ) 

8  to  10  Hz 

Flexure  Modal  Damping(^) 

~3  %  (or  critical)  without  force  loop, 
up  to  >100%  with  force  loop  on 

Stroke 

±0.2  inches 

Stroke  at  10  Hz  for  3  lbs.  output 
force 

0.1  inches 

Gravity  Offset  Spring  Rate 

2.4  Ib/in 

Allowable  Strut  Diameter 

1.00010.01” 

LPACT  Envelope 

3.8”  OD  X  4.86”  height  (including  strut 
clamp  and  accelerometers) 

LPACT  Total  Weight 

4.0  lb. 

LPACT  Proof  Mass  Weight 

2.9  LB 

LPACT  Model  (low  frequency) 
(refer  to  Figure  2  for  measured 
FRF  from  current  to  force  of 
LPACT) 

Output  Force{lb) 

Current  Conimand(a»2/?) 

+  2^CO„S  +  G>„ 

Servo  Amp  Model 

Current  (amp)  ^  ^  ^  amp 

Servo  Command  (V)  ^  V 

Force  Loop  Model 
(see  section  3.3  for  definition  of 
terms) 

Servo  Amp  Voltage  Command  (volts) 

Pr  oof  Mass  Accel  (g) 

_  ^pre  ^rt^force® 

(s  + Wpre)(s  + Wrt)^ 

Rate  Loop  Model 
(see  section  3.3  for  definition  of 
terms) 

Servo  Amp  Voltage  Command  (volts) 

Primaiy  Accel(g;) 

_  ^ pre^n^rate^rate^r  ^ 

(■S  +  ^pre)is  +  W^)^(5  + 

Table  12.  LPACT  Characteristics 


Cable  Assembly  cable 

connect  to  LPACT  Electronics 

(all  on  rear  panel) 

connect  to  LPACT 

Component 

Black  Coax 

‘To  Coil’ 

(banana  plug  to  BNC  adapter) 

6”  Blue  Pigtail  from  coil 

(BNC) 

Blue  Coax  marked  with 

Red  Tape 

‘From  Secondary 
Accelerometer’  (BNC) 

Secondary  Accelerometer 
on  Proof  Mass  (microdot) 

Blue  Coax 

‘From  Primary  Accelerometer’ 

(BNC) 

Primary  Accelerometer 
on  Co-Locate  Ring  (microdot) 

Table  13.  LPACT  Electronics  Connectivity  Guidelines  [From  Ref.  16] 


FRONT 


LPACT  1 

LPACT  2 

TP  ENABLE 

0 

DISABLE 

Input 

LP  i^ABLE 

0 

DISABLE 

mput 

ON 

i 

OFF 

Outputs 

®  ®  ® 

Outputs 

0  <0  0 

Current  Primary  Secondary 
Command  Accel  Accel 

Current  Primary  Secondary 
Command  Accel  Accel 

LED 


Figure  56.  LPACT  Control  Electronics  Rear  Panel  [From  Ref.  16] 
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Current 


where  m  =  mass  of  proof  mass  (2.9  lb) 

Figure  57.  LPACT  System  Level  Block  Diagram  [From  Ref.  16] 


PCB  Piezotronics  Model  208B02  General  Purpose  ICP  Force  Sensor  S/N 15021 


Sensitivity 

(Specification) 

50  mV/lb  (1 1240  mV/kN) 

Sensitivity  (Measured) 

50.80  mV/lb  [Ref  24] 

Dynamic  Range 

-10  lb  to  100-lb 

Stiffiiess 

1.0  kN/fim 

Temperature  Range 

-54  to  121  "C 

Sensing  Element 

Quartz 

Table  14.  PCB  Model  208B02  Operating  Characteristics  [Ref  15] 

PCB®  Piezotronics  Type  484B  Signal  Conditioner  S/N  2086 

Notes 
Unity  Gain 

Set  CPLG  to  DC  &  Bias  to  6  V 


Kistler  Instrument  Corp.  Accelerometers: 

(Note:  g  =  9.807  m/s^) 


Serial  Number 

+  x-axis 

+  v-axis 

+  z-axis 

8690C10 

Cl  12398 

495 

490 

494  mV/g 

8690C10 

Cl  12399 

487 

490 

490  mV/g 

8690C10 

Cl  12400 

499 

500 

494  mV/g 

8690C10 

Cl  12401 

497 

491 

505  mV/g 

Kistler  Instrument  Corp.  Signal  Conditioners  (Multi-Channel  Couplers): 


Type  Serial  Number 

5 124A  (twelve  channel)  C74930 
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Trek  Voltage  Amplifier: 


Type  Serial  Number 

Trek  50/750  none 


Notes 

Required  calibration  on  10  September  1998. 
Mini-cal  to  verify  setting  and  linearity 
Performed  by  author  on  2  July,  1999 
Two  channels  that  can  be  used  with  both  active 
struts. 


dSPACE  GmbH,  Germany,  S/N  3192 

Motorola  PowerPC  running  at  300  Hz 
Memory  128MB  RAM 

Input  Channels  20  ADC 

Ouput  Channels  8  D AC  0- 1 0 V 


Hewlett  Packard  HP  33120A  Signal  Generator  S/N 

settings 

Hewlett  Packard  HP  54601A  Digital  Oscilloscope  S/N  3134A02713 
Four  Channels 

Frequency  response  up  to  100  MHz 

Hewlett  Packard  HP  35665A  Dynamic  Signal  Analyzer  S/N 
hookups  and  settings 


SOFTWARE  DOCUMENTATION 


Program 

Version 

MATLAB 

5.2.1430 

dSPACERTIllOS 

dSPACE  MLIB/MTRACE 

3.1.1  ■ 

dSPACE  ControlDesk 

3.0 

ANSYS/Multiphysics 

5.5.2 

Table  15.  Software  Documentation 
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APPENDIX  B.  TRUSS  MODES.M 


% - % 

%  Modified  from  Example  8.9,5, "The  Finite  Element  Method  using  MATLAB"% 


%  Young  W.  Kwon  &  Hyochoong  Bang,  1997  CRC  Press  LLC  % 
%  Modified  by  LT  Carey  M.  Pantling  % 
%  Modified  by  LT  Timothy  A.  Barney  % 
%  Eigenanalysis  of  NPS  Space  Truss  % 
%  Bulids  systems  mass  and  stiffness  matrices  for  integating  with  % 
%  active  controls  and  plots  the  mode  shapes. 

% 


%  Variable  descriptions  % 
%  coord  ~  global  x,y  and  2  coordiates  of  each  node  % 
%  nd  =  nodal  connection  vector  % 
%  k  =  element  stiffness  matrix  % 
%  kk  =  system  stiffness  matrix  % 
%  m  =  element  mass  matrix  % 
%  mm  =  system  mass  matrix  % 
%  ff  -  system  force  vector  % 
%  index  =  vector  containing  system  dofs  associated  with  each  element% 
%  bcdof  ~  vector  containing  dofs  associated  with  boundary  conditions% 
%  bcval  =  vector  containing  boundary  condition  values  associated  % 
%  with  the  dofs  in  ’ bcdof ^  % 
% - - - - - % 


%  NPS  Space  Truss 

clear; 

nel=161; 

nnel=2 ; 

ndof=6; 

nnode=52 ; 

sdof =nnode*ndof ; 


%  number  of  elements 

%  number  of  nodes  per  element 

%  number  of  dofs  per  node 
%  total  number  of  nodes  in  system 

%  total  system  dofs 


load  nps_coord.dat; 
load  nps_node.dat; 
load  nps_bcdof.dat  ; 


%  52  by  3,  52  nodes  x,y,z 

%  161  by  2,  161  elements,  nodel  node2 

%  24  by  2,  24  degrees  constrained 


coord=nps_coord; 
nd=nps_node ; 
bcdof =nps_bcdof ( : ,  1)  ; 
bcval=nps_bcdof ( : , 2) ; 


%  all  units  in  standard  metric  units 

el=70E9;  '%  GPa  elastic  modulus 

area=l. 9686982E-5;  %  m^2  cross-sectional  area 

xi-1 . 242217E-10;  %  m^4  moment  of  inertia  of  cross-section 

rho=2800;  %  kg/m^3  mass  density  per  volume 


ff=zeros (sdof, 1) ; 
kk=zeros (sdof , sdof) ; 
mra=zeros (sdof , sdof ) ; 
index=zeros (nel*ndof , 1) ; 


%  initialization  of  system  force  vector 

%  initialization  of  system  matrix 

%  initialization  of  system  matrix 

%  initialization  of  index  vector 
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oV>  oV> 


for  iel-l:nel  %  loop  for  the  total  number  of  elements 

index=feeldof (nd(iel, :) ,nnel,ndof) ;  %  extract  system  dofs  associated 

%  with  the  element 

%  compute  element  stiffness  and  mass  matrices 
[k,m]  =feframe3 (el, xi, area, rho, coord (nd (iel, 1) , :) , coord (nd (iel, 2) ,:)); 
kk=feasmbll (kk, k, index) ;  %  assemble  element  matrix  into  system  matrix 
mm=feasmbll (mm,m, index) ; 
end 

%  add  the  concentrated  masses  at  the  nodal  translation  dof 
for  iel=l:nnode 

starts (iel-l) *ndof+l; 
fini=start+2 ; 

%  add  lumped  masses  for  the  node  points 

mm (start : f ini, start : fini) =mm (start : f ini, start : f ini) + (0 . 0663) *eye (3) ; 
%  kg  added  to  each  node  as  concentrated  mass 

end 


%  apply  boundary  conditions 

%  remove  for  simulink  analysis,  as  this  prevents  m  matrix  inversion 
[kk,mm,  ff  ]  =feaplyc2  (kk,mm,  f f  ,bcdof  ,bcval)  ; 

[V,  D]  =eig  (kk,  mm)  ; 
f sol=diag (D) ; 

for  index=l ; size (f sol) 

output (index, : ) = [index  f sol (index) ] ; 

end 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  Plots  the  first  4  mode  shapes.  % 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  First  Mode  Shape 
figure (1) ; 
amp=2 . 0 ; 
mode=V( : ,256) ; 
disp=nps_coord; 
for  index=l :nnode 

start= (index-1) *6+1; 

disp (index, 1) =disp (index, 1) +amp*mode (start)  ;  %  x-coord 

disp(index,2)=disp(index,2)+amp*mode(start+l)  ;  %  y- 
disp(index,3)=disp(index,3)+amp*mode(start+2) ;  %  z- 

end 

disp=:real  (disp)  ; 

%  form  the  line  combinations 

plot3 (disp (1:2,1)  , disp (1:2,3)  ,disp(l:2,2)  ,  'b*-»)  ; 
hold  on 

plot3  (nps__coord(l:2, 1)  ,nps_coord  (1 : 2 , 3)  ,nps_coord  (1 :2 , 2)  ,  ’r.  •  )  ; 

view(-30,20) ; 

axis (  [-2  3  -12  -12]) ; 

plot3 (disp (3 :14, 1) , disp (3 : 14 , 3 ) , disp (3 : 14 , 2 ) , 'b*-') ; 
plot 3 (nps_coord (3:14,1) , nps_coord (3:14,3) , nps_coord (3:14,2) , ' r . ' ) 
view(-30,20) ; 
axis ( [-2  3  -12  -12]) ; 


plots (disp (15 :26, 1) , disp (15 :26, 3) ,disp (15 :26, 2) , ’b*- ' ) ; 

plots (nps_coord (15:26,1) , nps_coord (15 : 26, 3) , nps_coord (15 :26,2),’r,’) 

view(-30, 20) ; 

axis ([-2  3-12-1  2]  )  ; 

plots (disp (27:28,1), disp (27:28,3), disp (27 : 28 , 2) , »b*- ’ ) ; 

plots (nps_coord (27:28,1) , nps^coord (27:28,3) , nps_coord (27 :28,2) , 'r. ’) 

view(-30, 20) ; 

axis ( [-2  3-12-12]) ; 

plots (disp (29:40,1) , disp (29:40,3), disp (29 :40, 2) , ’b*- * ) / 

plots (nps_coord(29:40, 1) ,nps_coord(29 :40, 3) , nps_coord (29 :40 , 2) , 'r. M 

view(-30, 20) ; 

axis( [-2  3-12  -1  2] ) ; 

plots (disp (41: 52, 1) , disp (41 : 52 , 3 ) , disp (41 : 52 , 2) , 'b*- ’) ; 

plots (nps_coord (41:52,1) , nps_coord (41 : 52 , 3) , nps_coord (41:52,2) , ' r . ’ ) 

view (-30, 20) ; 

axis ( [-2  3-12-12]) ; 

%  cross  ties 
for  index=l:26 

line (1, : ) =disp (index, : ) ; 

line  (2,  ; )  =:disp  (index+26,  : )  ; 

plots  (line  ( : ,  1)  ,  line  ( :  ,  3)  ,  line  ( : ,  2)  ,  'b*  -  ' )  ; 

end 

%  verticals 

for  index=  [1  2  27  28] 

line (1, : ) =disp (index, : ) ; 
line (2 , : ) =disp (index+7, : ) ; 

plots (line ( : , 1) , line ( : , 3) , line ( : , 2) , ’b* - ’ ) ; 

end 

for  index=  3 : 14 

line (1, : ) =disp (index, : ) ; 

line  (2  ,  : )  =:disp  (index+12  ,  : )  ; 

plots  (line  ( : ,  1)  ,  line  ( : ,  3)  ,  line  ( : ,  2)  ,  'b*-  ' )  ; 

end 

for  index=  29:40 

line (1, : ) =disp (index, : ) ; 

line (2 , : ) =disp (index+12 , : ) ; 

plots (line ( : , 1) , line ( : , 3) , line ( : , 2) , ’ b*- ' ) ; 

end 

hold  off 

title ( ’Mode  1 ’ ) ; 

print  -djpeg  mode4ml.jpg 

%  Second  Mode  Shape 
figure (2) ; 
amp=2 . 0 ; 
inode=V(  :  ,257)  ; 
di  sp=nps_coord  ; 
for  index=l :nnode 

start= (index-1) *6+1; 

disp  (index,  1)  =disp  (index,  1)  +aTnp*Tnode  (start)  ;  %  x- coord 

disp (index, 2) =disp (index, 2) +amp*mode (start+1) ;  %  y- 
disp  (index,  3)  =disp  (index,  3) +aTnp*Tnode  (start+2)  ;  %  z- 

end 

disp=real  (disp)  ; 

%  form  the  line  combinations 
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plots (disp(l:2,l) ,disp(l:2,3) ,disp{l:2,2) , 'b*-') ; 
hold  on 

plots (nps_coord (1:2,1) , nps_coord ( 1 : 2 , S ) ,nps  coord(l:2,2) 

view (-SO, 20) ; 

axis ( [-2  S  -12  -12]) ; 

plots (disp(S:14,l) ,disp(S:14,S) ,disp(S:14,2) , 'b*-') ; 

plots (nps_coord (S : 14 , 1) , nps_coord (S : 14 , S) ,nps_coord (S : 14 , 2) , ' r . ' ) ; 

view (-30, 20) ; 

axis ( [-2  3-12-12]) ; 

plots (disp(15:26,l) , disp (15 : 26, 3) , disp (15 : 26, 2) , 'b*- ') ; 

(^PS_coord (15:26,1) , nps_coord (15 : 26, 3 ) , nps  coord (15 :26,2)  'r.  ') 

view(-30,20) ; 

axis ( [-2  3  -12  -12]) ; 

plots (disp{27:28,l) , disp (27 : 28 , 3 ) , disp (27 : 28 , 2) , 'b*-') ; 

plots (nps_coord (27:28,1) ,nps_coord (27:28,3) , nps  coord (27 :28 , 2) ,  '  r  ' ) 

view(-30,20) ; 

axis (  [-2  3  -12  -12]) ; 

plots (disp(29:40, 1) , disp (29 : 40 , 3 ) ,disp (29 :40 , 2) , 'b*-') ; 

plots (nps_coord (29:40,1) ,nps_coord (29:40,3) , nps  coord (29 :40,2) ,  'r.  ') 

view (-30, 20 ) ; 

axis (  [-2  3  -12  -12]) ; 

plots (disp (41:52,1) ,disp (41:52,3) ,disp (41:52,2) , 'b*- ' ) ; 

P^°t3 (nps_coord (41:52,1) , nps_coord (41:52,3) , nps  coord (41:52,2)  '  r .  ' ) 

view (-30, 20 ) ; 

axis ( [-2  3  -12  -12])  ; 

%  cross  ties 
for  index=l : 2  6 

line (1, : ) =disp (index, : ) ; 
line (2, : ) =disp (index+26, : ) ; 

plots (line(:,l) ,line(:,3) ,line(:,2)  ,  'b*-') ; 

end 

%  verticals 

for  index=  [1  2  27  28] 

line (1, : ) =disp (index, : ) ; 
line (2, : ) =disp (index+7, : ) ; 

plots (line( : , 1) ,line( : ,3) , line( : ,2) , 'b*- ') ; 

end 

for  index=  3 : 14 

line (1, : ) =disp (index, : ) ; 

line (2 , : ) =disp (index+12, : ) ; 

plots (line(: ,1) ,line(: ,3) ,line(: ,2) , 'b*- ') ; 

end 

for  index=  29:40 

line (1,  : ) =disp (index, : ) ; 
line (2 , : ) =disp (index+12, : ) ; 

plots (line ( : , 1) ,line ( : ,3) , line ( : , 2) , 'b*- ' ) ; 

end 

hold  off 

title ( ' Mode  2 ' ) ; 

print  -djpeg  mode4ml.jpg 

%  Third  Mode  Shape 
figure (3) ; 
amp=2 . 0 ; 
mode=V(: ,258)  ; 
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di sp=np s_coor d ; 
for  index=l innode 

start= (index-1) *6+1; 

disp (index, 1) =disp (index, 1) +amp*mode (start) ;  %  x- coord 

disp  (index,  2)  =disp  (index,  2)  +anip*mode  (start+1)  ;  %  y- 
disp (index, 3) =disp (index, 3) +amp*mode (start+2) ;  %  z- 

end 

disp=real (disp) ; 

%  form  the  line  combinations 

plot3 (disp (1 : 2,1), disp (1:2,3) , disp (1 :2 , 2) , 'b*- O ; 
hold  on 

plot3 (nps_coord (1 :2 , 1) , nps_coord (1 : 2 , 3 ) , nps_coord (1:2,2) , ' r . ' ) ; 

view (-3  0, 20)  ; 

axis( [-2  3  -12  -12]) ; 

plot3  (disp (3:14,1)  , disp (3:14,3) , disp (3 : 14 , 2 )  ,  ’ b* - ’ ) ; 

plot3 (nps_coord (3 : 14 , 1) ,nps_coord (3:14,3) ,nps_coord (3 : 14 , 2) ,  *  r . ’ ) ; 

view(-30, 20) ; 

axis ([-2  3-12-1  2]  ) ; 

plot3 (disp (15:26,1) , disp (15: 26, 3) , disp (15 : 26 , 2 ) , ’b*-') ; 

plot3  (nps__coord(15  :26, 1)  , nps_coord  (15  : 26, 3 )  ,  nps_coord  (15  :26, 2)  ,  ’  r .  ' )  ; 

view(-30, 20) ; 

axis( [-2  3  -12  -12]) ; 

plot3 (disp (27:28,1) , disp (27 : 28 , 3 ) , disp (27 : 28 , 2 ) , ’b*-’) ; 

plot 3 (nps_coord (27:28,1) , nps_coord (27:28,3) , nps_coord (27:28,2) , 'r. ') ; 

view(-30, 20) ; 

axis (  [-2  3-12-12]) ; 

plot3 (disp (29:40,1) , disp (2 9 : 40 , 3 ) , disp (29 : 40 , 2 ) , 'b*-’) ; 

plot3  (nps__coord(29:4  0,l)  ,  nps_coord  (2  9  : 40 , 3 )  ,  nps_coord  (29  :40 , 2)  ,  'r.  ’ )  ; 

view(-30,20) ; 

axis (  [-2  3  -12  -12]) ; 

plot3 (disp (41: 52, 1) ,disp(41:52,3) , disp (41 : 52 , 2 ) , ^b*- ») ; 

plot3 (nps_coord (41:52,1) ,nps_coord (41:52,3) , nps_coord (41:52,2) , ' r . ') ; 

view(-30, 20)  ; 

axis  (  [-2  3  -12  -12]) ; 

%  cross  ties 
for  index=l:26 

line  (1,  : )  =disp  (index,  ; )  ; 
line  (2  ,  : )  =:disp  (index+2  6,  : )  ; 

plot3  (line  (  : ,  1)  ,  line  (  :  ,  3)  ,  line  (  :  ,  2)  ,  'b*-  ' )  ; 

end 

%  verticals 

for  index=  [1  2  27  28] 

line (1, : ) =disp (index, : ) ; 
line (2 , : ) =disp (index+7 , : ) ; 

plot3 (line ( : , 1) , line ( : , 3) , line ( : , 2) , ’b*- ’ ) ; 

end 

for  index=  3 : 14 

line (1, : ) =disp (index, : ) ; 
line (2 , : ) =disp (index+12 , : ) ; 

plot3 (line ( : , 1) , line ( : , 3) , line ( : , 2 ) , ’ b* - ’ ) ; 

end 

for  index=  29:40 

line (1, : ) =disp (index, : ) ; 

line (2 , : ) =disp (index+12 , : ) ; 

plot3 (line ( : , 1) , line ( : , 3) , line ( : , 2) , 'b*- ' ) ; 
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end 

hold  off 

title { ' Mode  3 ' ) ; 

print  -djpeg  mode4ml.jpg 

%  Fourth  Mode  Shape 
figure (4) ; 
amp=2 . 0 ; 
mode=V( : ,259) ; 
di sp=nps_coord ; 
for  index=l :nnode 

start= (index-1) *6+1; 

disp (index, 1) =disp (index, 1) +amp*mode (start) ;  %  x-coord 

disp (index, 2) =disp (index, 2) +arap*mode (start+1) ;  %  y- 
disp  (index,  3)  =disp  (index,  3) +anip*mode  (start+2)  ;  %  z- 

end 

disp=real (disp) ; 

%  form  the  line  combinations 

plot3 (disp (1:2,1)  ,disp(l:2,3)  ,disp(l:2,2)  ,  'b*-')  ; 
hold  on 

plot3 (nps_coord (1 :2 , 1) ,nps_coord (1 :2 , 3) , nps  coord(l:2,2) , 'r  ')■ 

view(-30,20) ;  “ 

axis ( [-2  3  -12  -12]) ; 

plot! (disp (3:14,1) ,disp (3:14,3) , disp (3:14,2), 'b*-'); 

plot3 (nps_coord (3 : 14 , 1) , nps_coord (3 : 14 , 3)  , nps  coord(3 : 14 , 2)  'r  ')• 

view(-30, 20) ;  ~  >  ■  , 

axis ( [-2  3-12-12]) ; 

plot! (disp (15 : 26, 1)  ,disp  (15:26,3) , disp (15 : 26 , 2)  ,  'b*- ' )  ; 
plot3 (nps_coord ( 15 :26, 1) , nps_coord (15 :26, 3) ,nps  coord (15 :26,2),'r  ')• 
view(-30,20) ;  “ 

axis ( [-2  3  -12  -12]) ; 

plot3 (disp (27:28,1) , disp (27:28,3) ,disp (27:28,2) , 'b*- ' ) ; 

plot3 (nps_coord (27:28,1) ,nps_coord(27 :28, 3) ,nps  coord (27 : 28 , 2) , 'r . ' ) • 
view(-30,20) ;  “ 

axis ( [-2  3  -12  -12]) ; 

plot3 (disp(29:40, 1) , disp (29 : 40 , 3 ) , disp (29 : 40 , 2 ) , 'b*-') ; 
plot3 (nps_coord(29:40,l) , nps_coord (29 :40 , 3) ,nps  coord (29 : 40, 2)  'r  ’) • 
view(-30,20) ;  “ 

axis ( [-2  3  -12  -12]) ; 

plot3 (disp (41:52,1) , disp (41:52,3) , disp (41:52,2) , ' b*- ' ) ; 
plot3 (nps_coord (41 : 52 , 1) , nps  coord (41 : 52 , 3 ), nps  coord (41: 52  2)  'r  ')• 
view (-30, 20) ;  “  ,  >  ,  ■  >  . 

axis (  [-2  3  -12  -12]) ; 

%  cross  ties 
for  index=l:26 

line (1, : ) =disp (index, : ) ; 

line (2, : ) =disp (index+26, :) ; 

plot3 (line(:,l) ,line(:,3) ,line(:,2)  ,  'b*-') ; 

end 

%  verticals 

for  index=  [1  2  27  28] 

line (1, : ) =disp (index, : ) ; 
line (2, : ) =disp (index+7, : ) ; 

plot3 (line ( : , 1) ,line ( : ,3) , line ( : , 2)  ,  'b*-  ' )  ; 

end 
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for  index=  3 : 14 

line (1, : ) =disp (index, : ) ; 
line (2 , : ) =disp (index+12 , : ) ; 
plots (line ( : , 1) / line ( : , 3) , line ( : , 2) , ’b*" * ) 
end 

for  index=  29:40 

line (1, : ) =disp (index, : ) ; 

line (2 , : ) =disp (index+12 , : ) ; 

plots (line ( : , 1) , line ( : , 3) , line ( : ,  2) , ’b*- ' ) 

end 

hold  off 
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APPENDIX  C.  NPS  SPACE  TRUSS  MODE  SHAPES 


The  first  four  mode  shapes  obtained  with  MATLAB  are  presented  first.  This  model  does 
not  include  the  added  mass  of  the  LPACT  or  the  active  elements  and  results  in  symmetric 
shapes.  The  deformed  shapes  are  shown  with  the  undeformed  node  locations  for 
comparison.  The  first  four  mode  shapes  generated  by  ANSYS  are  presented  second,  with 
both  the  deformed  and  imdeformed  shapes.  The  ANSYS  model  includes  the  LPACT  and 
active  element  masses  and  stiffiiess  characteristics. 


Mode  1  Mode  2 


Mode  3  Mode  4 

Figure  58.  NPS  Space  Truss  Mode  Shapes  with  MATLAB 
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The  mode  shapes  generated  by  ANSYS  are  shown  with  both  the  deformed  and 
imdeformed  states.  The  undeformed  states  appear  as  light  gray  lines  in  the  following 
figures. 


MODE  1 
(9.427  Hz) 
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MODE  2 
(11.318  Hz) 


MODE  3 
(12.529  Hz) 


MODE  4 
(18.744  Hz) 


Figure  59.  NFS  Space  Truss  Mode  Shapes  with  ANSYS 
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APPENDIX  D.  CONTROLLER  OPERATION 


To  start  and  run  the  controller,  the  following  steps  must  be  followed  (all  on  the 
NPS  space  truss  computer  in  the  Smart  Structures  Laboratory): 

•  Initiate  ControIDesk  using  the  ^  icon  either  on  the  desktop  of  the  toolbar,  or 
by  using  the  menu  Start>Programs>dSPACE>ControIDesk 

•  Open  the  Experiment  using  the  command  File  >  Open  Experiment.  Select 
C  :\Space_Truss\TrussO  1  VAdaptive  Controller.cdx. 

•  Load  the  controller  to  the  CPU  by  using  the  ■  icon  on  the  toolbar.  Select 
the  program  C:\Space_Truss\Truss01\Finished.ppc. 

•  Ensure  the  Edit  Mode  is  selected  from  the  toolbar  (verify  the  Edit  mode 

;I  °S  ?BS  °5 

button  depressed)  -l-l  *'  ""  I.  The  three  buttons  are  Edit,  Test  and 
Animate  mode  respectively. 

•  Start  the  CPU  (depress  the  green  triangle  on  the  toolbar). 

•  Select  the  Animate  mode  to  enable  the  display. 

•  The  disturbance  and  controller  may  be  turned  on  and  off  with  the  push  buttons 
as  previously  described. 


TfflS  PAGE  INTENTIONALLY  LEFT  BLANK 
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APPENDIX  E.  TRUSS_CONTROL.TXT 

!  *  Truss  Control  APDL  Program  for  ANSYS  version  5.6  * 

!  *  Constant  Frequency,  19  Hz  Disturbance  from  LPACT  * 

!  *  Applied  to  the  NPS  Space  Truss  Active  Controlled  Model  * 

!  *  Controlling  the  Z-Axis  of  Node  49  [Lab  Node  18,  Y-Axis]  * 

!  *  Using  Piezo  #2  [Voltage  Applied  to  Node  86]  * 

!  *  Written  by  LT  Carey  M.  Pantling  * 

!  *  Modified  by  LT  Timothy  A.  Barney  * 

!  *  With  assistance  by  Sheldon  Imaoka  * 

!  *  Last  Modified  2  March  2001  * 

finish 

!  First  Load  Truss  model  with  mesh  and  EC's 

!  *  Increase  the  Max  Result  Set  Size  * 

/config,nres,  1 0000 

!  *  Define  Variables  * 

*set,LMSADAPT,-0.045 
*set,FREQ,19.0 
*set,PV,86 
! 

*set,CTLNODE,49 


*set,STRTTIME,0.0 

*set,FINITIME,5.0 

*set,DDISPl,0.0 

*set,DDISP2,0.0 

*set,DVELOl,0.0 

*set,DVELO2,0.0 


!  Control  Node  Displacement  Value 
!  LPACT  Disturbance  Frequency 
!  Node  to  apply  Piezo  Voltage 
(Node  67  =  Piezo  #1,  Node  86  =  Piezo  #2) 
!  Set  Control  Node  Sensor  Location 
(Node  41  =  Lab  Node  26) 

(Node  18  =  Lab  Node  49) 

(Node  49  =  Lab  Node  18) 

(Node  26  =  Lab  Node  41) 

!  Start  time  =  0.0  seconds 
!  Finish  Time  in  seconds 

!  Previous  Nodal  Displacement 
!  Ciurent  Nodal  Displacement 
!  Previous  Nodal  Velocity 
!  Current  Nodal  Velocity 
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*set,CS,0.0 

*set,G3,0.0 


!  Control  Signal  to  Piezo,  initial  zero 
!  Control  Signal  Amplification  to  Piezo 


TIMESTEP— l/FREQ/20  !  Defines  the  Timestep 

*set,PI,acos(-l) 

*set,APLNODE,75  !  The  Node  Where  the  LPACT  is  Located 

finish 

!  *  Generate  Disturbance  Signal  at  20  Steps  per  Period  * 

j  ***♦**♦♦♦***♦**♦♦♦**♦**♦♦**♦♦♦*!)=*♦***♦*♦♦*** 


*set,NUMSTEP,nint(FINITIME/TIMESTEP) 

*set,FORCFUN, 

*set,FFP, 

*dim,FORCFUN, array jNUMSTEP  !  Array  for  Forcing  function 
•  Used  as  Input  (indexed) 

♦dim, FFP, table, NUMSTEP  !  Table  for  Plotting 

!  Used  as  Output  (not  indexed) 


-  Calculate  Disturbance  Magnitude  Based  on  FREQ  and  LPACT  Response  - 


♦set, MAG, 0.05437*FREQ**2-1 .8874*FREQ+1 8.2439  !  (N)  force  from  LPACT 


-  Create  a  Table  for  the  Disturbance  Force  - 


♦do,ICOUNT,  1  ,NUMSTEP 

FORCFUN(ICOUNT)=sin(2^PI*FREQ^((ICOUNT-l)^TIMESTEP))^MAG 
FFP(ICOUNT)=sin(2^PPFREQ^((ICOUNT- 1  )*TIMESTEP))^MAG 
♦enddo 


♦set,NODEDISP 

♦dim,NODEDISP,table,NUMSTEP 

♦set,CSCheck 

♦dim,CSCheck,table,NUMSTEP 


-  Included  for  Output  Data  - 


♦set,MYDISP, 


Tracks  the  Control  Node  Displacement 
Creates  an  Array  to  Get  Displacement 
Tracks  the  Control  Output  Signal 
Gets  the  Control  Signal 


!  Tracks  the  Node  Displacement 
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*set,MYVELO, 

*set,MYACEL, 

*set,MY  COSINE, 

*set,MYSINE, 

*(iini,MYDISP, table, NUMSTEP 
*diin,MYVELO,table,NUMSTEP 
*dim,MYACEL,table,NUMSTEP 
*dim,MYCOSINE, table, NUMSTEP 
*dmi,MYSINE,table,NUMSTEP 

!  *  Do  loop  for  Each  Time  Step  * 


!  Tracks  the  Node  Velocity 
!  Tracks  the  Node  Acceleration 
!  Tracks  the  Cosine  Weighting  Factor 
!  Tracks  the  Sine  Weighting  Factor 
!  Gets  the  Node  Displacement 
!  Gets  the  Node  Velocity 
!  Gets  the  Node  Acceleration 
!  Gets  the  Cosine  Weighting  Factor 
!  Gets  the  Sine  Weighting  Factor 


-  Clears  any  Constraints  on  Nodes  67  and  86  - 


/prep7 

ddele,67,volt 

ddele,86,volt 

finish 


-  ANSYS  Model  Solution  - 


*do,ICOUNT,l , NUMSTEP 

*set,CURRTIME,STRTTIME+ICOUNT*TIMESTEP  !  Set  Current  Time 
!  (Actually  the  End  of  Current  Time  Step) 

/solu  !  Enter  the  Solution  Processor 

*if,ICOUNT,eq,  1  ,then 

antype,trans,new  !  Start  New  Transient  Analysis 

*else 

antype,trans,rest  !  Restart  or  Continue  Transient  Analysis 

*endif 

time,CURRTIME  !  Set  Current  Time  for  ANSYS  Solution 

deltim,TIMESTEP  !  Set  Delta-t 

f,APLNODE,fy,FORCFUNaCOUNT)*0.707  !  Apply  Force  to  APLNODE 

Node  (LPACT) 

f,APLNODE,fz,FORCFUN(ICOUNT)*0.707  !  Apply  Force  to  APLNODE 

Node  (LPACT) 

d,PV,volt,CS  !  Apply  CS  Voltage  to  PV  Node 

allsel,all  !  Select  Everything  to  Prepare  for  Solve 

solve 
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finish 


-  Turns  on  Controller  After  Initial  Transient  - 


*if,ICOUNT,ge, 200, then  !  Turns  on  Control  at  Specified  Step 

G3=20.0 
*endif 

*if,ICOUNT,le,200,then 

*set,SINEWT,0.0 

*set,COSINEWT,0.0 

*endif 


-  ANSYS  Post  Processor  - 


/postl  j  Enter  Post  Processor 

*get,NODEDISP(ICOUNT),node,CTLNODE,u,z  !  Get  the  Control  Node  Displacement 


!  *  Control  Law  implementation,  for  next  time  step  * 

j  ♦*****************s|c***********************:jj,(.**,|5,,e,^,^ 


*set,DDISPl,DDISP2 

Displacement 

*get,DDISP2,node,CTLNODE,u,z 
Displacement 
*set,DVEL01,DVEL02 
*set,D  VEL02,(DDISP2-DDISP  1  )/TIMESTEP 
*set,DACEL,(DVEL02-DVEL01)/TIMESTEP 
Acceleration 

MYDISP(ICOUNT)=DDISP2 

MYVEL0(IC0UNT)=DVEL02 

MYACEL(ICOUNT)-DACEL 


!  Saves  Last  Control  Node 

!  Reads  Current  Control  Node 

!  Saves  Last  Control  Node  Velocity 
!  Calculates  Current  Mean  Velocity 
!  Calculates  Current  Mean 

!  Tracks  Displacement  Values 
!  Tracks  Velocity  Values 
!  Tracks  Acceleration  Values 


-  Updates  Control  Signal  - 


SINEWT-SINEWT+LMSADAPT*DACEL*sin(2*PI*FREO*ffICOUNT- 

1)*TIMESTEP))  !  Sine  Weight 
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COSINEWT=COSINEWT+LMSADAPT*DACEL*cos(2*PI*FREQ*((ICOUNT- 
1)*TIMESTEP))  !  Cosine  Weight 

TERMl=SINEWT*sin(2*PI*FREQ*((ICOUNT-l)*TIMESTEP))  !  Sine  Term 
TERM2=COSINEWT*cos(2*PI*FREQ*((ICOUNT-l)*TIMESTEP))  !  Cosine  Term 
CS=G3*(TERM1+TERM2)  !  New  Control  Voltage 

MYCOSINE(ICOUNT)=COSINEWT  !  Tracks  the  Cosine  Weighting  Factor 
MYSINE(ICOUNT)=SINEWT  !  Tracks  the  Sine  Weighting  Factor 

*if,CS,gt,50,then  !  Voltage  Limiter,  Based  on  Actual  Piezo 

CS=50 
*endif 

*if,CS,lt,-50,then 

CS=-50 

*endif 

CSCheck(ICOUNT)=CS  !  Tracks  the  Control  Signal  for  Output 

*enddo  !  Ends  the  Loop,  Ready  for  Next  Time  Step 

!  *  Plot  the  Results  to  the  Screen  * 

I  ****** ***♦♦*♦***♦ 

/erase 

*vplot„NODEDISP(l) 

J  ♦:|c*:j£:|c***************  ************  ********** 

!  *  Generates  Output  Files  for:  * 

!  *  Control  Signal  -  "zCSou^ut.out"  * 

!  *  Control  Node  Displacement  -  "zNdisp.out"  * 

!  *  Control  Node  Velocity  -  "zNvelo.out"  * 

!  *  Control  Node  Acceleration  -  "zNacel.out"  * 

!  *  Cosine  Weighting  Factor  -  "zCosine.out"  * 

!  *  Sine  Weighting  Factor  - "  zSine.out"  * 

I  ****************************************** 

*cfopen,zCSoutput,out, 

*vwrite,CSCheck(l), 

(E10.3) 

*cfclos 

*cfopen,zNdisp,out, 

*vwrite,MYDISP(l), 

(E10.3) 

*cfclos 
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*cfopen,zNvelo,out, 

*  vwrite,MYVELO(  1 ), 
(E10.3) 

*cfclos 

*cfopen,zNacel,out, 

*vwTite,MYACEL(l), 

(E10.3) 

*cfclos 

*cfopen,zCosine,out, 
*vwrite,MYCOSINE(  1 ), 
(El  0.3) 

*cfclos 

*cfopen,zSine,out, 

*vwrite,MYSINE(l), 

(El  0.3) 

*cfclos 


APPENDIX  F.  ACT  TRUSSJNP 


SjC  9^  9{c  9|C  3^  9^  9|C  3^  9^  9jC  9^  9|C  9|C  3|C  9}c  9|!  9^  9|C  9{C  9^  9|C  9|C  9^  9)C  9fc  9fc  9|C  9|C  9{C  9^  9^  9|C  ^C  ^C  9{C 

*  act_truss_3.inp  * 

*  used  in  getting  the  fresh  truss  on  line  * 

*  Truss  Control  APDL  Program  for  ANSYS  version  5.5  * 

*  For  simple  sinusoidal  disturbances  * 

*  Applied  to  the  NPS  Space  Truss  Active  Controlled  Model  * 

*  Written  by  LT  Carey  M.  Pantling  * 

*  Last  Modified  14  Oct  1999  * 


!  First  Load  Truss  model  with  mesh  and  EC's 

!  *  Define  Variables  * 

(  :^:lc4:9ie^^3)csic:jcHe9ieHe4:^9ie^4s^aie9)c 


eplot 


!  gives  something  to  look  at  while  waiting 


!  *set  command  will 
*set,FREQ,  11.75 
*set,CS,0.0 
*set,PV,68 
*set,RS, 

*set,SV,60 

*set,Gl,300 

*set,G2,100 

*set,G3,0.0 

*set,OLDRS,0.0 

*set,OLDINT,0.0 

*set,OLDFINT,0.0 

*set,OLDDBL,0.0 

*set,OLDFDBL,0.0 

*set,STRTTIME,0.0 


define  and  clear  variables 
!  disturbance  frequency 
!  control  signal  to  piezo,  initial  zero 
!  Piezo  voltage  node#  for  application 
!  read  signal  from  piezo,  defined  only 
!  Sensor  voltage  node#  for  detection 
!  Control  gain  constants,  see  below 


ItheoldRS,  IC  0.0 
!  old  integral,  IC  0.0 
!  old  filtered  integral 
!prevDI,IC0.0 
!  old  filtered  DI 
!  Start  time  =  0.0  seconds 


!  set  to  more  reasonable  later 
*set,FINITIME,4.0  !  Finish  Time  in  seconds 

TIMESTEP=l/FREQ/20  !  TIMESTEP 
*set,PI,acos(-l) 

*set,DISPNODE,25  !  node  26  real  truss  for  reading  the  output 

*set,APLNODE,53  !  the  node  where  the  LPACT  is  located 

finish 
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j  **************Hc5tc*ais****5ie*5lc5|e**si«**:jcHe*****!i::|s*s!«***:|cHcH«5le**Jj«*5|c5|c*5(c 

!  *  Set  up  load  disturbance,  at  least  20  steps  per  period  * 


*set,NUMSTEP,nint(FINITIME/TIMESTEP) 

*set,FORCFUN, 

*set,FFP, 

*dim,FORCFUN,array,NUMSTEP  !  array  for  forcing  function 

*dim,FFP,table,NUMSTEP  !  table  for  plotting 

!  arrays  for  input  (indexed),  tables  for  output(non  indexed) 

!  create  a  magnitude  based  upon  the  frequency,  from  the  LPACT  chart 

♦set, MAG,0.05437*FREQ**2-1.8874*FREQ+1 8.2439  !  (>0  force  from  LPACT 

!  Make  a  table  for  the  force  disturbance 
*do,ICOUNT,  LNUMSTEP 

FORCFUN(ICOUNT)=sin(2*PI*FREQ*((ICOUNT-l)*TIMESTEP))*MAG 

FFP(ICOlJNT)=sin(2*PI*FREQ*((ICOUNT-l)*TIMESTEP))*MAG 

*enddo 


*set,NODE26Y 

*dim,NODE26Y,table,NUMSTEP 

*set,CSCheck 

*dim,CSCheck,table,NUMSTEP 

*set,RSCheck 

*dim,RSCheck,table,NUMSTEP 


creates  an  array  to  get  displacement 
tracks  the  control  output  signal 
gets  the  control  signal 
tracks  the  read  signal 
copies  the  sensed  signal 


!  these  may  be  deleted,  and  not  be  required  to  be  copied,  included  for  error  checking 
!  *set,INTcheck 

!  *dim,INTcheck,table,NUMSTEP 
!  *set,DBLcheck 

!  *dim,DBLcheck,table,NUMSTEP 
!  *set,FINTchk 

!  *dim,FINTchk,table,NUMSTEP 
!  *set,FDBLchk 

!  *dim,FDBLchk,table,NUMSTEP 

j  sH**************************************** 

!  *  Do  loop  for  loading  at  each  time  step  * 

I  ♦♦******************iit*******************,^ 


*do,ICOUNT,  1  ,NUMSTEP 
!  Set  current  time  (actually,  end  of  current  step) 
*set,CURRTIME,STRTTIME+ICOUNT*TIMESTEP 
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!  enter  the  solution  processor 
!  Start  new  transient  analysis 
!  Restart  or  continue  transient  analysis 


/solu 

*if,ICOUNT,eq,  1  ,then 
antype,trans,new 
*else 

antype,trans,rest 
*endif 

time,CURRTIME  !  set  current  time  for  ANSYS  solution 

deltim,TIMESTEP  !  set  delta-t 

f,APLNODE,fy,FORCFUN(ICOUNT)*0.707  !  apply  force  to  APLNODE  node 
(LPACT) 

f,APLNODE,fz,FORCFUN(ICOUNT)*0.707  !  apply  force  to  APLNODE  node 
(LPACT) 


d,PV,volt,CS 

allsel,all 

solve 

finish 

*if,ICOUNT,gt, 200, then 
G3=20.0 
*endif 


!  apply  CS  voltage  to  PV  node 
!  select  everything  to  prepare  for  solve 


!  turns  on  control  at  specified  step 


/postl 

*get,RS,node,SV,volt 

RSCheck(ICOUNT)=RS 


!  enter  post  processor 
!  get  sensor  voltage 
!  saves  cvurent  RS  in  table 


*get,NODE26Y(ICOUNT),node,DISPNODE,u,y  !  get  the  tip  displacement 
finish 


!  *  Control  Law  implementation,  for  next  time  step  * 


*set,INT,OLDINT+(RS+OLDRS)/2*TIMESTEP  !  approx  the  integral 

*set,FINT,0.95*OLDFINT+INT-OLDINT  !  digital  high  pass  filter 

*set,DBLINT,OLDDBL+(FINT+OLDFINT)/2*TIMESTEP  !  second  integral 
*set,FDBL,0.95*OLDFDBL+DBLINT-OLDDBL  !  digital  high  pass  filter! 

!  INTcheck(ICOlJNT)=INT  !  copies  the  terms  for  examination 

!  FINTchk(ICOUNT)=FINT 
!  DBLcheck(ICOUNT)=DBLINT 
I  FDBLchk(ICOUNT)=FDBL 

OLDRS=RS  !  cycle  the  terms  to  old  values 

OLDINT=INT 

OLDFINT=FINT 

OLDDBL=DBLINT 

OLDFDBL=FDBL 

CS=(G1*FINT+G2*FDBL)*G3  I  CS=RS*(Gl/s+g2/s^2)*G3 
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!  voltage  limiter 


*if,CS,gt,60,then 
CS=60 
*endif 

*if,CS,lt,-60,then 
CS=-60 
*endif 

CSCheck(ICOUNT)=CS  !  save  for  output 

•'  ends  the  loop,  ready  for  next  time  step 

!  *  Plot  the  results  to  screen  * 

j  ******s|«*5|c*s!c*%**^H«****sic* 

/erase 

*vplot„NODE26Y(l) 


[  ********************,(:***4:;jt*^.;^^„;,^,H,|,,U^^j^j^j|.j|jjj.jjjj|^^^ 

!  *  Write  RSCheck  results  to  text  file  "ATR_e#_r#.out"  * 

!  *  and  NODE26Y  results  to  ATN.  * 

!  *  and  CSCheck  results  to  ATS.  * 

J  *********5|c*5lc***4:5ic5iC5!c**5!c**9jc**5ic3|c3ic:|c**Hc*3!e5|c*5|c:!cjjc3lesJ:3(c** 

*cfopen,AR_e  l_rl  ,out, 

*  vwrite,CSCheck(  1 ), 

(El  0.3) 

*cfclos 

*cfopen,AN_el_rl  ,out, 
vwrite,NODE26Y(l), 

(E10.3) 

*cfclos 

*cfopen,AS_el_rl  ,out, 
vwrite,NODE26Y(l), 

(E10.3) 

*cfclos 
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APPENDIX  G.  EXPERIMENTAL  RESULTS 


Node  41 ,  Z-Axis,  1 5  Hz,  Piezo  #1 
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Figure  60.  Exp.2,  Controller  Gain  Response 


Node  41 ,  Z-Axis,  15  Hz,  Piezo  #1 


Time  (sec) 

Figure  61.  Exp.  2,  System  Response 
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Time  (sec) 

Figure  64.  Exp.  3,  System  Response 


Node  26.  X-Axis,  12  Hz,  Piezo  #2 
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Figure  65.  Exp.  3,  Power  Spectrum 


147 


10  15  20  25  30  35  40 

Frequency  (Hz) 


0.15 


Node  26,  Y-Axis,  2nd  Modal  Freq.  Piezo  #2 
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Figure  72.  Exp.  6,  Controller  Gain  Response 


Node  26,  Y-Axis,  2nd  Modal  Freq,  Piezo  #2 
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Figure  73.  Exp.  6,  System  Response 
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Frequency  (Hz) 

Figure  74.  Exp.  6,  Power  Spectrum 

Node  26.  Z-Axis,  2nd  Modal  Freq,  Piezo  #2 


Time  (sec) 


Figure  75.  Exp.  7,  Controller  Gain  Response 
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Figure  80.  Exp.  9,  Controller  Gain  Response 
Node41.Z-Axis,  13.81  Hz,  Piezo  #1 
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Figure  81.  Exp.  9,  System  Response 
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APPENDIX  H.  PROCESS1.M 


% - ^ 

%  Plots  the  ControlDesk  capture  files  from  dSPACE  Experimental  runs  % 
%  using  the  NPS  Space  Truss.  % 
%  Written  by  LT  Timothy  A.  Barney  % 
% 


%  Variable  descriptions 

%  zCosine  =  Weighting  value  of  the  Cosine  LMS 
%  zSine  =  Weighting  value  of  the  Sine  LMS 
%  zGain  =  Adaptive  Gain  value 

%  zCSoutput  =  Controler  Output  Signal  to  Piezo 
%  zNacel  “  Nodal  Acceleration  Being  Controlled 


% 

% 

% 

% 

% 


%  DEMO  #1 

%  Controlling  the  Y-axis  of  Node  18,  with  a  19  Hz  disturbance,  using 
Piezo  #2 


clear  all; 
close  all; 


load  demol.mat 
freq=19;  % 

stop=180;  % 

dt=l/500;  % 

sdt=l/freq/20;  % 

stime= [sdt ; sdt : stop] ; % 
time=trace_x;  % 


Simulation  Disturbance  Frequency  (Hz) 
Simulation  Stop  Time  (sec) 

Downsampled  from  1  kHz  durring  capture 
Time  Step  Used  in  Simualtion 
Simulation  Time  Vector 
Experimental  Time  Vector 


zCosine=trace_y (4, : ) ;  %  Weighting  value  of  the  Cosine  LMS 

zSine=trace_y (5, : ) ;  %  Weighting  value  of  the  Sine  LMS 

zCSoutput =trace_y (1, :) ;  %  Controler  Output  Signal  to  Piezo 

zNacel=trace_y (2 , : ) ;  %  Nodal  Acceleration  Being  Controlled 

zGain=trace_y (3, ;) ;  %  Loads  Displacement  at  the  Control  Node 

zNacel=zNacel./0.05;  %  Converts  Acceleration  Signal  into  g's 

zCSoutput=zCSoutput . *200; %  Controler  Output  Signal  to  Piezo 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%  Plots  the  control  siganl  and  the  response.  % 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


figure (1) ; 
subplot (2 , 1, 1) 
hold  on 

plot (time , zNacel , *  r ’ ) 

title ('Node  18,  Y-Axis,  19  Hz,  Piezo  #2’) 

ylabel ( 'Acceleration  (g) ' ) 

xlabel ( 'Time  (sec) ' ) 

subplot (2,1,2) 

plot (time, zCSoutput, 'b') 

xlabel ( ' Time  (sec) ' ) 
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dp  dp 


ylabel (' Control  Signal  (Volts)') 
hold  off 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%^^^%^^^ 
%  Shows  the  convergence  of  the  LMS  weights.  % 


figure (2) ; 
subplot (3,1,1) 
hold  on 

plot (time, zCosine, 'g') 

title ('Node  18,  Y-Axis,  19  Hz,  Piezo  #2') 

ylabel ( ’ Cosine  Gain ' ) 

xlabel('Time  (sec)') 

subplot (3,1,2) 

plot (time, zSine, 'm' ) 

xlabeK'Time  (sec)') 

ylabel ('Sine  Gain') 

subplot (3,1,3) 

plot (time,zGain, 'b' ) 

xlabel('Time  (sec)') 

ylabel ( ' Gain ' ) 

hold  off 


%  Determines  Frequency  Response  Before  and  After  Control  % 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%^%^^^^^^^^^^^^^ 


tl=5*dt; 

t2=16.8; 

t3= (stop-t2) +6*dt; 
rl=floor (tl/dt) : floor (t2/dt) ; 
r2=floor (t3/dt) ; length (zNacel) ; 
vl=zNacel (rl) ; 
v2=zNacel (r2) ; 

fl=linspace(0,l/dt,length(rl) ) ; 
f 2=linspace ( 0 , 1/dt , length (r2 ) ) ; 
before=f ft (vl ' . *hanning (length (vl) ) ) ; 
af ter=f f t (v2 ' . *hanning (length (v2 ) ) ) ; 
avl=20*logl0 (abs (before) ) ; 
av2=20*logl0 (abs (after) ) ; 
rrl=find (f l<=50S:f  l>0)  ; 
rr2=find(f2<=50&f2>0) ; 
strt=fl (rrl (1) ) ; 
fnsh=fl (rrl (length (rrl) ) ) ; 
big=ceil (max (max (avl) ,max(av2) ) ) ; 
ltl=floor (mean(avl) -3*std(avl) ) ; 


%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%^^^%^^^^^^^^^^^ 
%  Plots  Frequency  Response  Before  and  After  Control  % 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%^^^%%^^^^^^^^ 


figure (3) 

siibplot  (2,1,1) 

plot (fl (rrl) ,avl (rrl) ) 
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axis ( [strt  fnsh  Itl  big] ) 

title (’Node  18,  Y-Axis,  19  Hz,  Piezo  #2’) 

ylabel ( 'Uncontrolled  (dB  [g] )  ' ) 

xlabel ( ' Frequency  (Hz) ’ ) 

subplot (2,1,2) 

plot (£2 (rr2) ,av2 (rr2) ) 

axis ( [strt  fnsh  Itl  big] ) 

ylabel ( ' Controlled  (dB  [g] ) ' ) 

xlabel ( ' Frequency  (Hz) ' ) 

text (5 , 20 , ’ 47 . 38  dB  Reduction  at  19  Hz') 

range=f ind (f l<=freq+l  &  fl>=freq-l); 
peak=avl (range) ; 

figure (4) 

plot (fl (range) ,avl  (range) , 'r: ' , f2 (range) ,av2 (range) , 'k- ' ) 
axis ( [fl (range (1) )  fl (range (length (range) ) )  Itl  big]) 
title ('Node  18,  Y-Axis,  19  Hz,  Piezo  #2') 
ylabel ( ' Power  Spectrum  (dB  [g]  )  ' ) 
xlabel ( ' Frequency  (Hz) ' ) 

text (18 . 2 , 30 , ' 47 . 38  dB  Reduction  at  19  Hz') 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%  Calculate  the  dB  reduction  at  the  peak  undamped  frequency  % 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


[m, i] =max (peak) ; 
reduction=m-av2 (range (i) ) 
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APPENDIX  1.  RESULT  COMPARISON 


Time  (sec) 


Figure  82.  Case  2,  Control  Node  Acceleration  Comparison 


Node  41 .  Z-Axis,  1 5  Hz,  Piezo  #1 


70 


^  3Q  I - 1 _ J - ^ _ I _ I _ I _ I _ \ _ I _ I 

0  20  40  60  80  100  120  140  160  180 

Time  (sec) 


Figure  83.  Case  2,  Control  Signal  Comparison 
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Figure  84.  Case  3,  Control  Node  Acceleratio 
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Figure  85.  Case  3,  Control  Signal  Comp'"- 
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control  (V)  Experimental  Control  (V)  Acceleration  (g) 
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Figure  87.  Case  4,  Control  Signal  Comparison 
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Node  26.  Y-Axis,  1  st  Modal  Freq.  Piezo  #1 
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Figure  88.  Case  5,  Control  Node  Acceleration  Comparison 
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Figure  89.  Case  5,  Control  Signal  Comparison 
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ANSYS  Acceleration  (g)  Experimental  Acceleration  (g) 


ANSYS  Control  (V)  Experimental  Control  (V)  . Acceleration  (g)  Experimental  Acceleration  (g) 


Node  26.  Z-Axis,  2nd  Modal  Freq,  Piezo  #2 
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Figure  92.  Case  7,  Control  Node  Acceleration  Comparison 
Node  26,  Z-Axis,  2nd  Modai  Freq,  Piezo  #2 
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Figure  93.  Case  7,  Control  Signal  Comparison 
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ANSYS  Control  (V)  Experimental  Control  (V)  ANSYS  Acceleration  (g)  Experimental  Acceleration  (g) 
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Uncontrolled,  Node  26,  Z-Axis,  3rd  Modal  Freq,  Piezo  #2 


Frequency  (Hz) 


Figure  96.  Case  8,  No  Control  Power  Spectrum  Comparison 
Controlled,  Node  26,  Z-Axis,  3rd  Modal  Freq,  Piezo  #2 


Frequency  (Hz) 


Figure  97 .  Case  8,  Controlled  Power  Spectrum  Comparison 
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APPENDIX  J.  OUTPUT.M 


% - 

%  Plots  the  ANSYS  output  files  generated  by  simulation  runs 
%  using  a  model  of  the  NPS  Space  Truss. 

%  Written  by  LT  Timothy  A.  Barney 
% 

%  Variable  descriptions 

%  zCosine  =  Weighting  value  of  the  Cosine  LMS 
%  zSine  =  Weighting  value  of  the  Sine  LMS 
%  zCSoutput  =  Controler  Output  Signal  to  Piezo 
%  zNacel  =  Nodal  Acceleration  Being  Controlled 
%  zNdisp  =  Displacement  at  the  Control  Node 

%  zNvelo  =  Velocity  at  the  Control  Node 

% - 


% 


% 


% 

% 

% 

% 


%  DEMO  #1 

%  Controlling  the  Y-axis  of  Node  18,  with  a  19  Hz  disturbance,  using 
Piezo  #2 


clear  all; 
close  all; 


f req=19; 
stop=5; 
dt=l/freq/20; 
time= [dt : dt : stop] ; 


%  Simulation  Disturbance  Frequency  (Hz) 
%  Simulation  Stop  Time  (sec) 

%  Time  Step  Used  in  Simualtion 
%  Time  Vector 


load  zCosine.out; 
load  zSine.out; 
load  zCSoutput .but; 
load  zNacel . out ; 
load  zNdisp.out; 
load  zNvelo.out; 
zNacel=zNacel . /9 . 8 ; 


%  Loads  Weighting  value  of  the  Cosine  LMS 
%  Loads  Weighting  value  of  the  Sine  LMS 
%  Loads  Controler  Output  Signal  to  Piezo 
%  Loads  Nodal  Acceleration  Being  Controlled 
%  Loads  Displacement  at  the  Control  Node 
%  loads  Velocity  at  the  Control  Node 
%  Converts  to  g's 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%  Plots  the  control  siganl  and  the  response.  % 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


figured)  ; 
siibplot  (2,1,1) 
hold  on 

plot (time , zNacel , *  r ’ ) 

title ('Node  18,  Y-Axis,  19  Hz,  Piezo  #2') 

ylabel ( 'Acceleration  (g) ' ) 

xlabel ( 'Time  (sec) ' ) 

subplot (2,1,2) 

plot (time, zCSoutput , 'b') 

xlabel ( 'Time  (sec) ' ) 

ylabel (' Control  Signal  (Volts)') 

hold  off 
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%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%^|.%^^;^^|.^ 
%  Shows  the  convergence  of  the  LMS  weights .  % 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%^%%^^^^^ 


figure (2) ; 
subplot (2, 1,1) 
hold  on 

plot (time, zCosine, 'g') 

title ('Node  18,  Y-Axis,  19  Hz,  Piezo  #2') 

ylabel ( 'Cosine  Gain') 

xlabeK'Time  (sec)') 

subplot (2,1,2) 

plot (time, zSine, 'm') 

xlabel('Time  (sec)') 

ylabel ('Sine  Gain') 

hold  off 


^^^^^%%^%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%^^^;^^^ 
%  Shows  the  control  node  displacement  and  velocity.  % 

%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%^^|.^^^^^^^^^^^^^^^ 


figure (3) ; 
subplot (2,1,1) 
hold  on 

plot (time, zNdisp, 'k') 

title ('Node  18,  Y-Axis,  19  Hz,  Piezo  #2') 

ylabel ( 'Displacement  (m) ' ) 

xlabeK'Time  (sec)') 

subplot (2,1,2) 

plot (time, zNvelo, ' c ' ) 

xlabel('Time  (sec)') 

ylabel ( 'Velocity  (m/s)') 

hold  off 
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APPENDIX  K.  COMPARE.M 


% - ^ 

%  Plots  the  experimental  and  FEM  results  for  comparison.  % 
%  Written  by  LT  Timothy  A.  Barney  % 
%  % 
%  Variable  descriptions  % 
%  eCSoutput  =  Experimental  Controler  Output  Signal  to  Piezo  % 
%  eNacel  =  Experimental  Nodal  Acceleration  Being  Controlled  % 
%  zCSoutput  =  ANSYS  Simulation  Controler  Output  Signal  to  Piezo  % 
%  aNacel  =  ANSYS  Simulation  Nodal  Acceleration  Being  Controlled  % 


%  DEMO  #1 

%  Controlling  the  Y-axis  of  Node  18,  with  a  19  Hz  disturbance,  using 
Piezo  #2 

clear  all; 
close  all; 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  Loads  the  Experimental  Data  % 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


eval ( [ ' cd  Experimental ’ ] ) 
load  demol.mat 


freq=19; 

stop=180; 

dt=l/500; 

time=:trace_x; 


%  Disturbance  Frequency  (Hz) 

%  Experiment  Stop  Time  (sec) 

%  Downsampled  from  1  kHz  during  capture 
%  Experimental  Time  Vector 


eCSoutput=trace_y (1, :)  ; 
eNacel=trace_y (2, :) ; 


%  Controller  Output  Signal  to  Piezo 
%  Nodal  Acceleration  Being  Controlled 


eNacel=eNacel . /O . 05;  %  Converts  Acceleration  Signal  into  g’s 

eCSoutput =eCSoutput . *200; %  Converts  control  signal  to  Voltage 


%%%%%%%%%%%%%%%%%%%%%%%% 
%  Loads  the  ANSYS  Data  % 
%%%%%%%%%%%%%%%%%%%%%%%% 


eval ( [ ’ cd  . . ' ] ) 
eval ( [ *  cd  ANSYS ’ ] ) 
load  zCSoutput .out; 
load  zNacel.out; 


%  Loads  Controller  Output  Signal  to  Piezo 
%  Loads  Nodal  Acceleration  Being  Controlled 


SStOp=5;  % 
sdt=l/freq/20;  % 
stime= [sdt :sdt :sstop] ;% 


Simulation  Stop  Time  (sec) 
Time  Step  Used  in  Simulation 
Simulation  Time  Vector 


aNacel=zNacel . /9 . 8 ;  %  Converts  to  g's 
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%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  Plots  the  Accelerations  % 
%%%%%%%%%%%%%%%%%%%%%%%%%%% 

figure (1) ; 
subplot (2, 1,1) 
hold  on 

plot (time, eNacel, 'r') 

titleCNode  18,  Y-Axis,  19  Hz,  Piezo  #2') 

ylabel ( ' Experimental  Acceleration  (g) ' ) 

xlabelCTime  (sec)') 

subplot (2, 1,2) 

plot ( stime , aNacel , ' b ' ) 

xlabel('Time  (sec)') 

ylabel ( 'ANSYS  Acceleration  (g)') 

hold  off 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  Plots  the  Control  Signals  % 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

figure (2) ; 
subplot (2,1, 1) 
hold  on 

plot (time, eCSoutput, 'r') 

titleCNode  18,  Y-Axis,  19  Hz,  Piezo  #2') 

ylabel ( ' Experimental  Control  (V) ' ) 

xlabel ( ' Time  (sec) ' ) 

subplot (2,1,2) 

plot (stime, zCSoutput, 'b') 

xlabel ('Time  (sec)') 

ylabel (' ANSYS  Control  (V)') 

hold  off 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%|.%%^%^^ 
%  Determines  Frequency  Response  Before  and  After  Control  % 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%^%^^^% 


etl=5*dt; 

et2=16.8; 

et3= (stop-et2) +6*dt; 
stl=5*sdt; 

St2=0.53; 

st3= (sstop-st2) +6*sdt; 


%  Experimental  FFT 
erl=floor (etl/dt) : floor (et2/dt) ; 
er2=floor (et3/dt) : length (eNacel) ; 
evl=eNacel (erl)  ; 
ev2=eNacel (er2) ; 

efl=linspace (0, 1/dt, length (erl) ) ; 
ef 2=linspace ( 0 , 1/dt , length (er2 ) ) ; 
ebefore=f f t (evl ' . ‘banning (length (evl) ) ) ; 
eafter=f ft (ev2 ' . ‘banning (length (ev2) ) ) ; 
eavl=20‘logl0 (abs (ebefore) ) ; 
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eav2=20*logl0 (abs (eafter) ) ; 

errl=f ind (ef l<=50&:ef  1>0)  ; 

err2=f  ind  (ef2<=50&:ef  2>0)  ; 

estrt=ef 1 (errl (1) ) ; 

ef nsh=ef 1 (errl (length (errl) ) ) ; 

ebig-ceil(max(max(eavl) ,max(eav2) ) ) ; 

eltl=floor (mean(eavl) -3*std(eavl) ) ; 

%  ANSYS  FFT 

srl=floor (stl/sdt) ; floor (st2/sdt) ; 
sr2=floor (st3/sdt) : length (aNacel) ; 
svl=aNacel (srl) ; 
sv2=aNacel (sr2) ; 

sfl=linspace (0, 1/sdt, length(srl) ) ; 

sf2=linspace{0,l/sdt, length (sr2) ) ; 

sbefore=fft(svl.*hanning (length (svl) ) ) ; 

safter=fft (sv2.*hanning (length (sv2) ) ) ; 

savl=20*logl0 (abs (sbefore) ) ; 

sav2=:2  0*logl0  (abs  (safter)  )  ; 

srrl=f ind (sf l<=50&sf 1>0) ; 

srr2=f ind (sf 2<=50&sf 2>0) ; 

sstrt=sf 1 (srrl (1) ) ; 

sfnsh=sfl (srrl (length (srrl) ) ) ; 

sbig=ceil  (Tnax(max(savl)  ,max (sav2) ) )  ; 

sltl=f  loor  (mean  (savl)  -•3*std  (savl) )  ; 

slt2=f loor (mean (sav2) -3*std (sav2) ) ; 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  Plots  Frequency  Response  Before  and  After  Control  % 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  Before  Control 

figure (3) 

subplot (2,1,1) 

plot  (ef 1 (errl) , eavl (errl) ) 

axis ( [estrt  efnsh  eltl  ebig] ) 

title (* Uncontrolled,  Node  18,  Y-Axis,  19  Hz,  Piezo  #2') 

ylabel ( ' Experimental  (dB  [g]  )  ' ) 

xlabel ( ’ Frequency  (Hz) ’ ) 

subplot (2,1,2) 

plot (sf 1 (srrl) , savl (srrl) ) 

axis ( [sstrt  sfnsh  sltl  sbig] ) 

ylabel { ’ANSYS  (dB  [g] ) ' ) 

xlabel { ’ Frequency  (Hz ) ’ ) 

%  After  Control 

figure (4) 

subplot (2,1,1) 

plot (ef2 (err2) ,eav2 (err2) ) 

axis ( [estrt  efnsh  eltl  ebig] ) 

title (’ Controlled,  Node  18,  Y-Axis,  19  Hz,  Piezo  #2’) 

ylabel ( ’ Experimental  (dB  [g] ) ’ ) 

xlabel ( ’ Frequency  (Hz) ' ) 

subplot (2,1,2) 

plot (sf2 (srr2) , sav2 (srr2) ) 
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axis ( [sstrt  sfnsh  slt2  sbig] ) 
ylabel ( 'ANSYS  (dB  [g] ) ) ' ) 
xlabel ( ' Frequency  (Hz) ' ) 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%^%^^^^^^^^^^^ 
%  Calculate  the  dB  reduction  at  the  peak  undamped  frequency  % 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%^^^^^^^^^ 


[em, ei] =max (eavl) ; 
ereduction=em-eav2 (ei) 


[sm, si] =max (savl) ; 
sreduction=sm-sav2 (si) 


eval ( [ ' cd  . . ' ] ) 
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APPENDIX  L.  ANSYS  RESULTS 
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Figure  98.  Simulation  1,  System  Response 


Node  18,  Y-Axis,  19  Hz,  Piezo  #2 
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Figure  102.  Simulation  2,  Controller  Gains 
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Figure  103.  Simulation  2,  Control  Node  Response 
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Velocity  (m/s)  Displacement  (m) 


Node  26,  X-Axis,  12  Hz,  Piezo  #2 
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Figure  108.  Simulation  4,  Controller  Gains 
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Figure  109.  Simulation  4,  Control  Node  Response 
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Figure  110.  Simulation  5,  System  Response 


Node  26,  Y-Axis,  1st  Modal  Freq,  Piezo  #1 


Figure  111.  Simulation  5,  Controller  Gains 
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Figure  112.  Simulation  5,  Control  Node  Response 
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Figure  1 13.  Simulation  6,  System  Response 
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Node  26,  Y-Axis,  2nd  Modal  Freq,  Piezo  #2 
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Figure  118.  Simulation  7,  Control  Node  Response 
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